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La nomenclatura usada en su mayor parte se adapta a la nomen- 


clatura aprobada por 1а Sociedad de Ingenieros de Petróleo de la AI 
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YACIMIENTOS DE PETROLEO 0-GAS CON EMPUJE HIDRAULICO 


A 
El término Ме que aparece en las ecuaciones de Balance de Ma- 


teriales!, representa la cantidad de agua que entra a un yacimiento de 
petróleo = pass proveniente de formaciones adyacentes saturadas de agua 
(acuíferos), como congecuencia de la disminución de presión en el con- 
tacto agua-petróleo, debido a la producción de fluidos del yacimiento. 
— 
La magnitud del término We está relacionada con la actividad del acuí- 
fero y forma parte de la energía total disponible para producir el pe- 
tróleo. Al término We se le denomina Intrusión de Agua y los yacimien 
tos que lo presentan se denominan Yacimientos con Empuje Hidráulico o 


Hidrostático- 


El presente trabajo comienza con una sección general sobre ya 
cimientos con enpuje hidráulico. Presenta luego la deducción y solu- 
ción de la ecuación de flujo transitorio © ecuación de difusividad apli 
cada a intrusión de agua, tanto para sistemas lineales como radiales. 
Continua con una discusión de los principales métodos para análisis de 
yacimientos con empuje hidráulico, referente a la evaluación de la in- 
trusión de agua, petróleo originalmente en el yacimiento, etc. Se dis 


cuten luego los métodos de predicción de este tipo de yacimientos. 
GENERALIDADES 


Empuje hidráulico es generalmente el mecanismo de recupera - 
ción primaria más eficiente en yacimientos de petróleo: Al mismo tiem 
po, algunos de los yacimientos más productores del mundo han sido pro- 
ducidos bajo este mecanismo. Entre muchos, pueden: mencionarse el cam- 


po East Texas, Conroe, Yates, Thompson y Friendwood en los Estados Uni 


dos; el campo Kirkut en Irak y los campos de Arabia Saudita. 


Acuíferos. E کو پو کا‎ ii. dese 


Los yacimientos con empuje hidráulico están total o parcial- 


mente comunicados con formaciones adyacentes totalmente saturadas de 
| 


agua, que reciben el nombre de acuíferos (Latin: aqua = agua y ferres 


i 

| 

\ llevar). Tales acuíferos pueden estar situados en la periferia de la 
formación productora del yacimiento y en este caso, se denoninan acuí 
feros laterales o de flanco, o pueden estar situados por debajo de la 

+ formación productora у еп tal caso se denominan acuíferos de fondo. 
Los mecanismos correspondientes 86 denominan empuje hidráulico late - 
ral y empuje hidráulico de fondo respectivamente. En los acuíferos la 
terales sólo una parte del área de la zona de petróleo está en contac 
to con el acuífero; еп cambio, en un acuífero de fondo, éste se encuen 
tra en contacto con la totalidad del área productora. La parte donde 
se une el acuífero y la zona petrolífera se denomina simplemente, con 


tacto agua-petróleo. 


En un acuífero lateral, el agua avanza por los lados o costa 

Н dos del yacimiento y generalmente su migración se dirige a las partes 
| altas de la estructura. Los pozos en la parte baja de la estructura son 
los primeros invadidos por el agua y la producción de agua aumenta rá 
pidamente en tales pozos. En cambio en los acuíferos de fondo, el a- 
gua se mueve en dirección vertical; por tanto, en un yacimiento cuya 
formación productora sea horizontal, los pozos serán inundados por a- 
gua a un mismo tiempo siempre y cuando se produzcan bajo condiciones 
Similares. Estos tipos de yacimientos con acuíferos de fondo han si- 


do estudiados por Muskat?, Hutchinson y Kenp?, Henley y colaborado - 


res*, Elkins?, Coats, etc. 


> 
! La forma de un acuífero puede ser muy variado, dependiendo 
у 


del tipo de estructuras y entrampamiento del sistema yacimiento acuífe- 
ro. Aparte de los acuíferos de fondo mencionados antes, muchos siste- 
mas yacimiento-acuífero pueden aproximarse bien sea a sistemas de tipo 
radial o de tipo lineal, para los cuales existen soluciones analíticas 
7,8,9, рага el cálculo de la intrusión. El tipo radial se encuentra 
más frecuentemente en estructuras tales como anticlinales y domos y el 


tipo lineal generalmente son resultado de entrampamiento por falla. La 


‚ Pig.1 ilustra un tipo común de acuífero radial, donde la intrusión ocu 


rre a través de toda la circunferencia. Ejemplos de este tipo de соп- 
portamiento pueden consultarse en la literatura! ®. Puede ocurrir que 
el flujo sea radial pero que el flujo ocurra sólo en una parte de 1а 


circunferencia; en tal caso, la intrusión será proporcional a la cir- 


cunferencia efectiva en que ocurra la intrusión. 


1. El tamaño del acuífero también puede variar mucho. Muskat 11 
considera que si el volumen poroso del acuífero es alrededor de 1009 
veces o más el volumen poroso del yacimiento, el acuífero puede consi- 
derarse de extensión infinita respecto al yacimiento, y usarse en esta 
forma en lo que respecta al comportamiento de flujo no-continuo. Tal 
tamaño corresponde, para acuíferos de forma radial, a un radio 461 а- 
cuífero de 32 veces el radio del campo de petróleo, es decir, тер = 32. 
Acuíferos де tamaño menor al mencionado se consideran de tamaño finito 
y es necesario conocer los límites del mismo con el propósito de eva - 
luar su radio o longitud respecto al radio o longitud del yacimiento. 

Sin embargo, la ecuación de difusividad, como:se verá más adelante, ha 


sido resuelta para muchos valores de Pep" 


Mecanismo del Empuje Hidráulico. , 


El principal factor que impulsa el agua a fluir a un yacimien 
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Fig. 1 - Ilustración de un sistemo yacimiento- acuífero tipo radial 


Presión, p — 


Fig. 2 ~ Distribución de presión con tiempo en un acuífero de extensión 


infinita. to < t1« 12 «13, 


to de petróleo es la expansión del agua del ácuífero. A medida que la 
presión disminuye en un yacimiento debido a su propia producción de hi 
drocarburos, se crea una presión diferencial entre el yacimiento y ela 
cuífero.. Cuando una presión diferencial (disturbio de presión) llega 
al contacto agua-petróleo, el agua del acuífero cerca al yacimiento со. 
mienza a expandirse a la presión menor a medida que el tiempc avanza, 
el disturbio viaja а través del acuífero y más agua se expande a la pre 
sión menor, aumentando al mismo tiempo la intrusión de agua al yacimien 
to. Esto se ilustra gráficamente en la Fig.2, donde Ri es el radio 101- 
cial exterior del yacimiento (o el radio inicial interior del acuífero) 
Auves de iniciar la producción del yacimiento, la presión esta represen 
сада por la linea to. Para tiempos posteriores tq, t5, tg c etc. las 
curvas de distribución de presión indican un aumento continuo del radic 
exterior de lates curvas. Si el acuífero es lo suficientemente grande. 
tal que el diferencial de presión no alcanza un límite exterior del a- 
сої его durante el periodo de interés, entonces se habla de un acuífer 
infinito. El tipo de flujo durante el cual el disturbio viaja a travé: 
del acuífero es estrictamente no-continuo. Finalmente el disturbio al 
canza el límite exterior del acuífero y a partir de este instante, еї 
flujo se convierte en semi-continuo (seudo-continuo o cuasi-continuo) 
simplemente continuo. La Fig.3 ilustra la distribución de presión par 
este caso, a diferentes tiempos, donde г, es el radio exterior del acu 


fero. Estos tipos de acuíferos se denominan limitados o finitos. 


Pueden exister acuíferos con otras formas diferentes de ener 
gíc al mecanismo descrito, considerado el más común. Algunas veces 1 
formación del acuífero зе extiende y aflora en superficie, donde una 
fuente de agua va reemplazando al agua que entra al yacimiento. Est 


: + : ^. д 
tipo de acuífero se denomina acuífero abierto. En este caso, existe 1 


PETROLEO»: 


Ri Radio, г — fe 


Fig. 3— Distribución de presión con tiempo en un ccuifero de extension limitada 
(finito). to «ti < t2 < t3 «ta«ts. 


COMUNICA A FUENTE EXTERIOR QUI 
SUMINISTRA AGUA AL ACUIFERO. 


PETROLE 


Presión, p — 


Ri Radio, г —* re 


Fig. 4- Distribución de presion con tiempo para un acuífero con suministro exte. 
rior де адио 10 <ti<t2<t3<te<ts. 


acuífero limitado (el radio exterior se extiende a la superficie) y el 

suministro de agua en superficie permite mantener una presión constan- 

te en el límite exterior del acuífero. Se habla entonces de flujo con 

tinuo o flujo arteciano. La distribución de presión en este caso se 

ilustra en la Fig.. La existencia de este tipo de empuje hidráulico 
12 


se ha descrito en la literatura . Generalmente, en este tipo de acuí 


ferp, ocurre un contacto agua-petróleo inclinado. 


Otro mecanismo que puede contribuir al suministro de agua al 
yacimiento es el debido al hundimiento de los estratos superpuestos. 
A medida que la producción aumenta, la superficie terrestre se hunde. 
Esto trae como concacuencia varias cosas:.compactación de lutitas veci 
nas al yacimiento y expulsión de agua de las mismas hacia el yacimien- 
to; compactación tanto de la formación productora (1o que también con- 
tribuye al aumento de la recuperación) y compactación del acuífero, lo 
que contribuye a aumentar el suministro de agua al yacimiento. Este 
tipo de hundimiento es muy común en varias partes del mundo, antea о- 
tras, еп la zona de Lagunillas, en. la Costa Bolivar del Lago de Mara - 


caíbo, Venezuela y la zona de Long Beach, Los Angeles, California, E.U. 


Е1 factor efectivo de expansión, o más comunmente conocido со 
mo el factor efectivo de compresibilidad del agua del acuífero, esta 
formado por el factor de expansión mismo del agua del acuífero y por el 
factor de compresibilidad de la roca. Aunque este último término, en 
valor absoluto es generalmente pequeño, puede tomar un porcentaje con- 
siderable en el valor efectivo, debido al mayor volumen de rocas y ё8- 
pecialmente en áreas donde ocurre hundimiento. Además, en casos en que 

13 


ocurra solubilidad de gas en el agua!?, ésta puede ser otra causa que 


contribuya a aumentar el factor efectivo de expansión, especialmente en 


zonas cercanas al yacimiento si han ocurrido acumulaciones de gas 11- 
bre durante la formación del acuífero. La obtención del factor de com 


presibilidad de la roca fué discutido en otro trabajol. 


Reconocimiento de un Empuje Hidráulico 


Antes de iniciar la explotación y durante la primera fase de 


pa мэн Е 5 EE PAPIRU HE ——À 


producción de un campo, no es posible asegurar la existencia efectiva 
о, no es posible asegurar са билана 21222120 


de un empuje hidráulico. Debe tenerse presente que los principales 
factores que conducen a un buen enpuje hidráulico son: alta permeabi 
lidad, que permita una fácil entrada del agua a la formazión y el ta- 
maño y continuidad del espesor del acuífero con distancia. Por ejem- 
plo, un acuífero delgado, de baja permeabilidad y que presenta lenti- 
cularidad es muy posible que no presente, O presente muy poco empuje 

hidrávlico. Debido a que la gran mayoría de las propiedades, tamafio 
y formación del аси его son desconocidas en el momento de descubrir 
un campo, es por consiguiente necesario esperar hasta tener disponible 


cierta historia de producción, con el fin de determinar su existencia 


y analizar su comportamiento. 


Existen varios indicativos con el fin de reconocer si ocurre 
empuje hidráulico. Comenzando con la perforación, puede decirse que 
si se localizan contactos agua-petróleo, (de registros, pruebas, etc),, 
es una buena indicación que puede ocurrir intrusión de agua una нэ? 
iniciada la producción. También, si la presión inicial del campo es- 
tá por encima del punto de burbujeo (indicación que no existe capa de 
gas inicial), lógicamente debe existir una fuente de energía que con- 


tribuyó al mantenimiento de tal presión en el campo. Tal fuente de e 


nergía puede ser un acuífero. 


A medida que la producción avanza, se dispone de mejor infor 


mación para reconocer si existe o no un acuífero activo, es decir, que 
ocurra intrusión de agua. El principal indicativo es el comportamien-. 
to de presión. Cuando la presión de un campo (o yacimiento) tiende 
a permanecer constante о disminuye muy poco, es indicación de un empu- 
je hidráulico activo. Generalmente, al.comienzo de la producción se 
presenta una caída de presión, ya que como antes se mencionó, es nece- 
sario crear una presión diferencial en el.contacto agua-petróleo con 
el fin de que comience la intrusión y ósto lógicamente, tarda cierto 
tiempo, que depende tanto del.tamafio del yacimiento como de la rata de 
producción. Durante este tiempo, al.no.entrar ningün fluido extra al 
yacimiento, la caída será mayor. .En.general, puede ‘decirse, que en el 
caso de acuíferos infinitos, la rata de disminución de presión es menor 
a medida que avanza el tiempo; en cambio,.en acuíferos limitados, la 
rata de disminución de presión puede ser constante con tiempo. - Tambien 
debe tenerse en cuenta que la caída de presión se hace menor por debajo 


del punto de burbujeo debido a la liberación de gas. 


A medida que se vaya disponiendo de historia de producción se 
recomienda calcular regularmente el petróleo originalmente en el yaci - 
miento por balance de materiales!, considerando que no existe intru- 
sión de agua, es decir, We = 0. Si el valor del petróleo originalmen - 
te en el yacimiento aumenta consecutivamente a medida que 1а produc - 
ción avanza, es también indicación que ocurre intrusión de agua ya que 
la ecuación de balance de materiales se ha considerado. sin el término 
de intrusión, lo que disminuye el valor calculado. Un gráfico del pe- 
tróleo calculado originalmente en el yacimiento, como función de tiem- 
po (o producción acumulada), permite obtener fácilmente el valor verda 
dero del petróleo originalmente en el yacimiento, simplemente leyendo 


el valor extrapolado, cuando el tiempo (о la producción acumulada) sea 
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cero. Una vez obtenido N, se regresa a la ecuación de balance de mate 


riales y se calcula la intrusión de agua, We. 


Finalmente, es indicación de un empuje hidráulico: (a) el va- 
lor relativamente constante del índice de productividad considerado en 
la base al campo como un total; (b) si la presión no disminuye por de- 
bajo del punto de burbujeo, la razón gas-petróleo de producción (ins - 
tantánea) permanece constante; (с) la producción de agua aumenta, espe 
cialmente si existen pozos cerca al contacto agua-petróleo. El estu - 
dio de estos pozos en particular, contribuye notoriamente al reconoci- 


miento de intrusión. 


En acuíferos de fondo, la producción de agua aumenta principal 
mente debido a la formación de conos de agua alrededor de los pozos $31 
la producción disminuye debido a prorrateo u otro factor, la presión de 


fondo en los pozos aumenta, ya que los efectos del acuífero continuan 


y tienden а restablecer un nuevo equilibrio en el yacimiento. 


Debe tenerse presente que mientras más rápido se determine la 
intrusión de agua a un campo, puede planearse mejor su desarrollo, ya 
que la localización de los pozos y el espaciamiento puede variar consi 


derablemente si existe o no un acuífero activo. 


Recuperación. 


El factor de recuperación, definiáo como el petróleo recupera 
do por асге-ріе de yacimiento o como porcentaje respecto al petróleo 
originalmente en el yacimiento, es un término afectado por muchas varia 
bles que dependen de las características de los fluidos y del yacimier 


to, lo mismo del arreglo de los pozos, etc. 


- 11 - 


Teniendo sólo consideraciones teóricas es bastante difícil in 
cluir en una fórmula todas las variables que influyen en el factordere 
cuperación. Se han estudiado por separado algunos factores principal- 
mente en procesos de invasión (inundación) con agua con fines de recu- 
peración secundaria. También se han hecho algunas investigaciones de 
la recuperación usando una combinación de datos teóricos y de laborato 
rio, pero los esfuerzos зе han dirigido principalmente a estudios de 
correlaciones y estadísticos usando datos de diferentes campos cuyo 
principal mecanismo de producción ha sido empuje hidráulico. En 1945, 
Craze y Buckley + publicaron un estudic estadístico, considerado hoy co 


mo clásico, sobre el efecto del espaciamiento de pozos en la recupera- 
ción en campos de petróieo, usando datos de 103 yacimientos. Posterior 
mente en 1955, Guthrie y Greenberger!? usaron la misma información de 
Craze y Buckley para ilustrer el uso de análisis de correlación en la 
obtención de la mejor relación entre la eficiencia de recuperación en 
yacimientos con empuje hidráulico y varias combinaciones de parámetros, 
considerados como los principales que contribuyen a la eficiencia de 
recuperación. Más recientemente, en Octubre de 1967, el Instituto Ame 
ricano del Petróleo, АРТ!6 (American Petroleum Institute) publicó el 
API BUL. D-1% sobre un estudio estadístico de la recuperaión, tanto de 
campos cuyo principai mecanismo de producción es empuje por gas en sò- 
lución como empuje hidráulico. Es un estudio hecho por el Sub-conitá 

Sobre Eficiencia de Recuperación del API, iniciado en 1956, y donde par 


ticipó prácticamente toda la Industria Norte Americana, Compañias Con- 


sultoras y algunas Universidades» 


El trabajo comprendió una evaluación y estudio de 312 yaci- 
mientos de Estados Unidos y otras partes del mundo. De estos yacimien 


tos, fueron seleccionados 72 que representaron una información acepta- 


E 


ble para análisis de empuje hidráulico en arenas y areniscas. Á traves 
de estudios estadísticos de análisis de regresión y diferentes modelos 
y usando computadores digitales para derivar correlaciones, fué posi - 
ble agrupar diferentes parámetros en una ecuación con el fin de deter- 
minar el factor de recuperación en tales yacimientos de empuje hidráu- 
lico. Es conveniente observar que el mencionado trabajo también inclu 


ye un estudio semejante para yacimientos con empuje por gas en solu- 


Para yacimientos con empuje hidráulico, el factor de recupera 
ción encontrado, en barriles normales de petróleo por acre-pie de for- 
mación productora, es, 


ets) 11588221 kuyi 9,0770 


4258|————— ——— 


B. мэн 
ВАР = A A ч шз ее (1) 
0,1903 Pi 0,2159 
(5,) 
Я Ра 


El mismo factor, expresado como un porcentaje respecto al petróleo o- 
riginalmente en el yacimiento, es, 


$(1-5,) |9398221 Kiwi 0,0770 


54,898 | ——— na 
Boi Moi 
RE -  _ _  — —__—__mm_—__ 00-00 (2)- 
0,1903[ pj 0,2159 
(5,9 : 
Pa 


donde рд es la presión de abandono. Los demás términos corresponden 
a la nomenclatura normal de la AIME!?. El subíndice i corresponde a 


condiciones íniciales. Debe tenerse cuidado en el uso de estas ecua 


2013: 


ciones, ya que ellas se aplican dentro determinado intervalo de validez 
para cada una de las variables que aparecen en las ecuaciones. Para va 
lores muy altos o muy bajos de las variables se recomienda consultar el 
estudio originali, A veces es más conveniente escribir estas ecuacio- 


nes en forma logarítmica, 


` д61-5,) kuyi 
log(BAF) = log(4259) + 1,0422 log|------^-- + 0,0770 log|------ E 
Boi Moi 
Pi 
- 0,1903 log(Sy) - 0,2159 log|---- (3) 
Pa 
| Ко: 
log(RE) = 1о6(58,898) + 0, 0422 log|- + 0,0770 log|------ 
Moi 
Pi 
- 0,1903 log(S,) - 0,2159 log|*--- (У) 
Ра 


El estudio en mención determinó un factor de recuperación pro 
medio, RE, de 51,1 por ciento para yacimientos con empuje hidráulico 
en areniscas y un valor de 43,6 por ciento para yacimientos con empuje 


hidráulico en calizas y dolomitas. 


La Ec.(1), o su equivalente la Ec.(3), se ha resuelto a través 
de un nomograma con el fin de facilitar los cálculos. Tal поподгапа a- 
parece en la siguiente página. El ejemplo incluido allí, ilustra cla - 
ramente el procedimiento seguido con el fin de obtener el factor de re- 


evperación, ВАР, en barriles normales por acre-pie de yacimiento. So- 


FACTOR DE 
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RECUPERACION, ВАР BARRILES NORMALES POR ACRE-PIE. 


EJEMPLO 


DATOS: 


Un yacimiento con empuje hidráulico tiene los siguien- 


tes datos: 


Porosigad, 4, 
Saturación одџа connata, Sw, 


Factor volumétrico iniciat Petróleo, Boj 1,289 BY/BN 


Permeabilidad absoluta, К, 
Viscosidad del petróleo a Pi, оа 
Viscosidad del одод а pi, Awi 
Presión inicial, pi, 

Presión da abandono, Pa, 


PROCEDIMIENTO: 


Se calcula los factores; 411 5w) 


K wi 


Hoi 


Entrando o la derecha соп 0,115 


hasta la linea de К 7/0 0,01; lusgo horizontolm. 


recuperación, BAP= 223 barriles por 


(77,58x0,115)= 25 porciento del patrdiso originaimente an al yacimiento. 


Nomoqromo Para Dater: 


4000 Ipco 


ITI s ону 


„10.06250,8) 


sigue horizontalmente hosta la пас de Sw= O, 


16,53 % 
20,00 5. 


0,0625 md. 
50 ср. 
0,8 ср. 


200 Ipca 


(0,1853)(0,8) 


7,289 


= 0,01 


5 


2; Se sigue verticalmente 
Біде = 20 y verticalmente ss les lo 
e un factor ae recuperación de 223/ 


acre pie. Este valor equivi 


minar el Foctor бе Racuperoción en Yacimientos con Empuje Hidráulico. 
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lo es necesario conocer los factores 4(1-5,)/8сү:» би» Крез! Ноћ Y Di/Pa: 
Siguiendo el camino indicado por las flechas, se lee directamente el 


valor de BAP. 


Para determinar el factor de recuperación, RE, como porcenta- 
je del petróleo originalmente en el yacimiento, basta dividir el resul 
tado obtenido del nomograma,BAP, por el factor 77,5849 (1-S4)/Bgis Pero 
ya que el término $(1-S,)/Boj ha tenido que calcularse para entrar al 


nomograma, es muy sencillo la obtención de RE. Analíticamente, 


RE 


Esta relación puede establecerse fácilmente comparando las Ec.s(1) y 
(2). El ejemplo incluido en el nomograma, ilustra la obtención de RE 


-а partir de BAP. 


Rata de Producción para Máxima Eficiencia. 


Este término se conoce comunmente en la literatura como MER 
(del ingles Maximum Efficient Rate). Es un término muy disctutidoi?;1! 
y debe tenerse en cuenta especialmente en yacimientos con empuje hi- 


dráulico. 


Es lógico suponer que en un yacimientos donde existe una intr! 


sión activa de agua, habrá una rata de producción de petróleo ideal, e! 
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tal forma que ocurra el máximo desplazamiento de petróleo del área ba- 
rrida por el agua. Por tanto, puede decirse que un yacimiento con em- 
puje hidráulico es sensitivo a la rata de producción de petróleo. Va- 
rios factores deben tenerse presente: (1) la intrusión de agua puede 
ser pequeña comparada con la producción de petróleo, permitiendo que о 
curra segregación, (2) la eficiencia de desplazamiento debido al agua 
de intrusión disminuye a medida que la rata de producción” aumenta, (3) 
a altas ratas de producción de petróleo, las fuerzas capilares pueden 

llegar a ser insignificantes, en tal forma que el desplazamiento capi- 
lar de petróleo de zonas de baja permeabilidad sea prácticamente nulo, 


(+) a medida que la rata de producción aumenta, еі contacto agua - ре ~ 


tróleo se distorciona, bien sea por conificación o inclinación. 


En general, la rata de producción para máxima eficiencia, pue 
de definirse como la máxima cantidad de petróleo por unidad de tiempo 
que puede permitirse producir de un yacimiento, sin pérdida significan 
te en la producción total. Tal rata de producción debe ser lo suficien 


temente baja, en tal forma que: 


(1) En su mayoría el agua de intrusión reemplace el 
petróleo producido, es decir, impida que el empuje hidráulico disminuya 


o desaparezca su acción. 


(2) Ocurra una alta eficiencia de desplazamiento en 


la zona invadida por agua de intrusión. 


(3) Las fuerzas capilares dezplacen el petróleo de 


zonas de baja permeabilidad. 


(^) El contacto'agua-petróleo:avance lo más hóri- 


zontalmente posible. 


DEDUCCION DE LA ECUACION DE DIFUSIVIDAD. 


En esta sección se hará la deducción de la ecuación de difusi 
vidad para 11441405 que fluyen a través de un medio poroso tanto en el 
sistema de coordenadas cartesianas con en.el sistema de corrdenadas ra 


diales. 


La ecuación de difusividad que describe matemáticamente elflu 
jo de una sola fase a través de un medio poroso, es la combinación de 
tres ecuaciones: (a) la Ecuación de Continuidad, que no es más que un 
balance de masa o una forma diferencial de la Ley de Conservación de 
Masa, (b) la Ley de Darcy y (c) la Ecuación de Estado, que describe el 


tipo de fluido en movimiento. 


Coordenadas Cartesianas. 


Considérese un elemento de un medio pororso homogóneo, talco 
mo el ilustrado en la Fig.5, donde ocurre flujo en tres direccíones, 
x, y, 2. Durante un instante cualesquiera de observación, At, deter- 
minado flujo entra al sistema a través de las secciones (áreas) (ћућа), 
(ухћа) y (Ахлу) y durante el mismo instante de observación, At, deter 
minado flujo sale a través de estas mismas áreas en las secciones de 
salida localizadas à distancias Ax, ^y Y Az de las secciones de entra 


da. 


Si м; (і = х, у, z) es la velocidad volumétrica de flujo (та 
ta de flujo. рог área unitaria, (L3/T)LÊ = L/T) en la divección iy Р 
la densidad del fluido, M/L?, el flujo de masa por área unitaria en la 
dirección i será puj = (M/T)/L?. Aplicando la Ley de Conservación de 


Masa a este elemento, puede escribirse para el tiempo de observación, 


At, 


[pue + 20021] 


PUx 


е 


dimensiones pora 
d «en coordenadas cortesionos. 


Fig. 5- Elemento de volumen de un medio poroso en tres 


deducir la ecuación de continuida 


[pur + А Gur] 


t+ ûr 


Fig. 6 - Elemento de volumen de un medio poroso para deducir la ecuación 


de continuidod en un sistema rodial horizental, 
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Cantidad de masa que Cantidad de masa que 
entra al elemento dul - sale del elemento du 
rante el tiempo At. rante el tiempo At. 


Cantidad neta de masa 


introducida (о retira Acumulación de masa en 
+ |да) al elemento por = lel elemento durante el 
fuentes o (sumideros) tiempo At. 


durante el tiempo At. 


Aplicando este ley a cada de las direcciones х, у, Z, la masa total 
que entra al elemento será el flujo de masa por área unitaria, pui, 
multiplicada por el área correspondiente, (AzAy), (А24х) o (Ахлу) por 
el instante de observación At. Aplicando lo anterior en la dirección 


x, la cantidad de masa que entra al elemento será, 
(рух) CAyaz)CAt) 


y la cantidad de masa que sale durante el mismo instante de observa- 


ción At, será, 
(ou, * Арц, CAyAz)CAt) 


Asumiendo que no entra ni sale fluido del elemento a través de fuentes 
o sumideros, la cantidad neta de masa acumulada en el elemento en la 
dirección x será la cantidad de masa que entra menos la cantidad de mas: 
que sale en el instante de observación At, o sea, 
A(puy) 
- [ACpuyJCAyaz)(5t)] = - Axbybz| —— —— 
Ах 
Haciendo balances de masa similares para las direcciones y, 2, la acu- 


mulación total en el elemento será la suma de tales acumulaciones, 


= BXxAyAz|————— + + |m At (5) 


' 
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Por otro lado, la acumulación de masa en el elemento durante 
el instante de observación At, puede obtenerse en forma diferente. A 


un tiempo t cualesquiera, la masa de fluido existente en el elemento es 
(ахдудлэ (40031, 
у а un tiempo t*át, la masa de fluido existente en el elemento será, 


(AxAyAz 1002], үс 


Por tanto, la acumulación de masa durante el instante At, será la can- 


tidad de masa al tiempo (t*At) menos la cantidad de masa al tiempo t, 


Џ(ахлућа)(%)(р)] - [Caxayaz) (0202), 


Asumiendo que las dimensiones del elemento no varian durante el At, la 


expresión anterior puede escribirse, 


Marte PERE - (pp) . (6) 


Ч 


Las Ecs.(5) у (6) individualmente representan la acumulación de masa 


en el elemento durante el instante At. Igualándolas y simplificando, 


А(роу) 111772) atpuz) (фр) фл = (ФР) 


И ыы 


Ax Ay Az |: At 


En el límite, cuando Ах, ду, Az y At tienden a cero, puede escribirse, 


сэл ME E ЛЭ A oue (7) 
— (ры + —(pu + (pu = - — (bp 
a. * ay Y az 1 at 


La expresión anterior es la ecuación de continuidad en coordenadas car 
tesianas para flujo monofásico a través de un medio poroso. 
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El siguiente paso en la deducción de la ecuación de difusi- 
vidad consiste en introducir la Ley de Darcy a la Ec.(7). Para el flu 


jo еп las direcciones х, у, 7, esta lev puede escribirse, 


МЕ ky ap 
* u эх 
y ў зу 
kz 3p 
ug = - + Pg 
м 32 
.Sustituyend» en la Ес.(7), 
a Куо эр 3 кур ФР a Кур эр 
— + — ве + eg) а 
9х u 9х эу у зу 92 u 92 
3 
— (Фр) 2 (8) 
3t 


Finalmente, la adición de la ecuación de estado en la Ec. 
(8), trae como resultado la ecuación de difusividad. Para obtener es- 
ta ecuación para líquidos, se parte de la definición del factor de com 
presibilidad isotérmica. Este factor, se define como el cambio en vo- 
lumen de fluido por cambio unitario de presión por volumen unitario. 


En forma analítica puede escribirse, 


1 av 
с = - |-— — — 
1) 3p JT 


o en función de densidad puede escribirse, 


1 80 


— س |د с‏ 


Po 
р 
с(р-ро) = zl | (9) 
Po 
P = Po "хаг А (10) 


La Ec.(9) es la ecuación de estado para un fluido de-compresibilidad 
constante. Extrictamente, el factor de compresibilidad es función 


de presión, pero para líquidos es razonable suponer 81 factor de com 


presibilidad, с, constante dentro del intervalo de presión de interés, 


De la Ec.(9), puede escribirse, 


Эр 
3р. 1 
эх cp эх 
Эр 1 др 
—- = —— (11) 
dy ср зу 
2+ Эр 1 Эр 
92 р Эр 


Efectuando las derivadas indicadas en la Ec.(8), sustituyendo los va 
lores de la Ec.(11), considerando la viscosidad congtante y despre - 


ciando las fuerzas de la gravedad, puede escribirse, 


32p 3?p э2р ЭР 2 ЭР , эр, 
k + ky + kz + c |kx( ) + куС--( к kz( у |+ 
эх? ay? az? эх dy az 


эр dky эр dky 3p 3k. Эр эф 
+ — + + سے‎ ——| = фис + سس ل‎ (12 
эх эх зу зу 92 dz 3t at 


Si además, c es un valor pequeño (en el orden de 1077 1рс-1), el medio 
es isotrópico,es decir, la permeabilidad y porosidad constantes en cual 
quier parte del sistema y los gradientes de presión son lo suficiente - 
mente pequeños para que los términos al cuadrado puedan despreciarse, la 


Ec.(12) puede escribirse, 


ә?р э2р 32p фис Әр 
92 = $ id + = سسس‎ — (1 
эх? эу? az? к э 


Esta es una de las formas más conocidas de la ecuación de difusividad 

que gobierna el flujo de una sola fase a través de un medio poroso, ex 
presada en ccordenadas cartesianas y Flujo en tres direcciones. La so 
lución de esta ecuación permite obtener la distribución de presión co- 
mo función de las tres direcciones х, у, 2, y a un tiempo t, cuales- 


quiera. 


Si el flujo es lineal, es decir, ocurre sólo en una sola di- 
rección, por ejemplo, en un plano horizontal,los términos en z e y son 


iguales a cero, y la ecuación de difusividad en este caso será, 


эр фис ap 
—_-. з —— = (14 
ax? k эт 


Esta es la ecuación básica aplicada a sistemas lineales de intrusión 


de agua, como se vará más adelante. 
Coordenadas Radiales 


Quizás de más uso en ingeniería de petróleo sea la ecua- 


ción de difusividad escríta para un sistema radial ortogonal, ya que 
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el flujo (o inyección) en un pozo que penetra completamente la forma - 
ción sigue condiciones de flujo radial. Para deducir tal ecuación en 
un sistema radial, puede partirse de la Ec.(13) en coordenadas carte- 
sianas y usando métodos puramente matemáticos (operador de Laplace, - 
etc), se hace la transformación al sistema radial. Sin embargo, puede 
hacerse una deducción semejante a la anterior para coordenadas carte - 
sianas. La Fig.6 ilustraun elemento del medio poroso en un sistema 

radial. El flujo se considera horizontal y radial en la misma direc- 


cción en que aumenta r, 


, 


Si up es la velocidad volumétrica de flujo radial, el flujo 
de masa por área unitaria será (pup). En forma similar a la deducción 
anterior, la masa que entra al clemento а travós del área lateral del 


cilindro de radio, r, durante un instante de observación, At: será, 


(pug) COrh) (At) 


y la cantidad de masa que sale del elemento a través del afea lateral’ 


del cilindro de radio, (r+Ar), será, 


[pup+ra(ou,)] [9(к+Аг)в] (4t) 


Por tanto, asumiendo que no sale ni entra masa al elemento a través de 
fuentes o sumideros, la cantidad neta de masa que se acumula en el ele 
mento durante el instante de observación At, será la diferencia de lo 


que entra, menos lo que sale , 


oh (ou, - (ou, * a(ou,)] (rran) }at (15) 
Por otro lado, la acumulación de masa en el elemento duran- 


te el instante de observación At, puede obtenerse en forma diferente. 
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A un tiempo t cualesquiera, la masa de fluido en el elemento es, 


(8rar)(h)(9)(0) 
t 


y a un tiempo (t+at), la masa de fluido en elemento será, 


(өрле) (в) ($) (р) 
trat 

por lo tanto, la acumulación de masa durante el instante АК, será la 
cantidad de masa al tiempo (t+At) menos la cantidad de masa al tiempo 
t, 

(ёгАг)(һ)(%ф)(р) - (6rar)(h)(9)(p) (16 

t+At t 
Las Ecs.(15) y (16) individualmente representan la acumulación de:masa 
en el elemento durante el instante At. Igualándolas, eliminando oh de 
? 


ambos lados, efectuando y transponiendo términos, puede escribirse, 


1 фрјежћт - $e), 
pupbr + rA(pu,) + NN | = س‎ {17 


At 


rár 


ya que el producto de los términos (Ar)(A(pup)) es muy. pequeño, pude 


considerarse igual a cero y la ecuación anterior puede escribirse, 


1 ACou,) А(фр) 
بے‎ рн зл o omm (1€ 
r àr At 


tomando el límite, cuando Ar y At tienden a cero, puede escribirse, 


1 2(pup) 3 
Т тен | e — (40) 21 
at 


r 3r 
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1 
— (трчи) = - (60) (20) 
r ar at 


Esta es la ecuación de continuidad para flujo radial. Si en esta ecua- 
ción se reemplaza la ecuación para la velocidad volumétrica de flujo, 


END dada por Darcy, 


Ka Эр 
Ч". u ar (21) 
resulta, 
1 3 rok, эр а 
س = )— — — س‎ (фр) (22) 
r эг y эх at : 


Finalmente, introduciendo la ecuación de estado para líquidos, Ec.(9), 


en su forma diferencial respecto а г, se obtiene, 


1 др др 1 Qk, эр Эр 12 фис Эр 
— — (Y ) + سسس | ن + — سس‎ a = — + 
r dr dr Ka dr ar эг k 3t 
u 39 
+ — {2 
Xp at 


Si se asume un medio isotrópico, es decir, que la permeabilidad y poro 
sidad sean constantes en cualquier parte y dirección del sistema, se 


puede escribir, 


1 ә 3p 3p ]? ap 1 3p ap 12 
— ا )س‎ )+ e| —— — + — + с = 


r dr ar | ar эг? г Эг дг 
фис Эр 
= — (24 
k ЈЕ 


Si además, зе considera el valor de с pequeño (en el orden de 1076 


1рс-1) y por otra parte puede asumirse un gradiente de presión peque 


= 97 ж 


ño ап cualquier parte del yacimionto, el término c(3p/àr), puede consi 
derarse igual a cero. En este caso, la ecuación resultante $3, 


| 3?p 1 3p фис Эр 
LOL وا‎ = == (25) 
| эг? r Br к эх 


Esta expresión es comunmente conocida como ecuación de difusividad para 


líquidos de baja compresibilidad. 


Debe tenerse en mente las suposiciones que incluyen las e - 
cuaciones deducidas, Ecs.(13) y (25): sistema isotrópico con porosidad 
y permeabilidad constantes» viscosidad del fluído constante, no se con 
sideran las fuerzas de la gravedad, flujo viscoso, compresibilidad pe- 
queña y constante, gradiente pequeño a través del sistema y una sola fa 
se. Ecuaciones similares pueden deducirse? para otros fluídos (gases) 


y cuando ocurre flujo de dos o tres fases (petróleo, agua, gas). 
SOLUCIONES DE LA ECUACION DE DIFUSIVIDAD APLICADA A ACUIFEROS 


Las soluciones de las ecuaciones de difusividad aplicadas a 
sistemas yacimiento-acuífero más comunmente conocidas, se refieren ala 
solución de las Ecs.(14) y (25), o sea modelos lineales y radiales ho- 
mogéneos en una sola dimensión y donde fluye una sola fase. Con el u- 
so de computadores digitales, día a día se incorporan a la literatura 
petrolera soluciones particulares (modelos) de las ecuaciones de difu- 
sividad en diferentes formas: una, dos y tres dimensiones, una dos y 
tres fases y cambio en ха ааааай: tales como porosidad, permeabili - 
dad, espesor, etc. En estas notas, se cubrirán las principales solu - 
ciones de las Ecs.(14) y (25) existentes еп la literatura, para siste- 
mas homogéneos y una sola fase y otras soluciones variando una o algu- 


nas de las propiedades tales como porosidad, permeabilidad, espesor, y 


compresibilidad. 
Acuíferos Lineales Homogéneos. 


Como se mencionó anteriormente en la sección de Generalida- 
des, algunas veces en la práctica, la geometría de un sistema yacimien 
to-acuífero se aproxima mejor a un sistema lineal que a un sistema ra- 
dial, principalmente cuando el sistema está limitado por fallas. La e- 


cuación aplicable en tal sistema es la Es.(14), 


3p фис Әр (14) 


ax? k at 


Con el fin de determinar las diferentes soluciones de esta ecua- 
ción, es necesario fijar primero el tamaño del sistema (acuífero), co- 
mo también las condiciones iniciales y las condiciones de límite (con- 
torno). Los tipos generales de soluciones, se refieren a acuíferos in 
finitos y a acuíferos finitos (limitados). Luego, a su vez, para cada 
tipo existen varias soluciones, de acuerdo a las condiciones iniciales 
y de límite. La Fig.9 presenta esquemáticamente las condiciones de 1í 
mite de las soluciones principales del sistema yacimiento-acuífero 11- 
neal, según los modelos representados en las Figs.7 y 8. La Fig.9 re- 
sume seis soluciones: dos para acuíferos infinitos y cuatro para acuí- 
feros finitos. Еп ambos casos, deben especificarse, tanto las condi - 


ciones iniciales, como las condiciones de límite. 


Condiciones Iniciales.- Tanto en acuíferos infinitos como en acuí- 
feros finitos, la condición inicial, es decir cuando el tiempo, (550, 


considera una presión inicial, constante para todo valor de х. Апа 


Pis 


líticamente,p(x,0) = pi. 
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ACUIFERO 


CONTACTO AGUA- PETROLEO, CAP, ` ~ 
O YACIMIENTO- ACUIFERO. 


LIMITE INTERIOR 


Fig. 7- Sistema lineol yncimiento- ocuífero infinito 


ACUIFERO 


х1 


А Ls 
CONTACTO AGUA -PETROLEO, CAP, LIMITE EXTERIOR 
O YACIMIENTO - ACUIFERO, 
LIMITE INTERIOR 


Fig. 8- Sistemo lineal yacimiento- ocuífero finito. 
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Condiciones de Límite (Contorno).- 


I - Acuíferos Infinitos 


Caso 1 - Rata de intrusión constante en el contacto agua-petroleo 
(acuífero-yacimiento). Es el denominado Caso de Rata Terminal Constan- 


te. Analíticamente,las condiciones de límite son:* 


kbh др 
a) (е) сар = constante з ГС сар рага t > 0 
( 3P ер 4 
) = = constante,para t > 

эх CaP  уъһ ? 

Эр ер 
—(0,t) = 

эх kbh 


b) р > р; cuando L + œ todo t. 


p(=,t)= Ра 


Caso 2 - Presión constante en el contacto agua-petróleo (acuífero- 


yacimiento). Es el denominado Caso de Presión Terminal Constante, Ana 


líticamente, las condiciones de límite son, 


a) (р) сар = constante,para t » 0 


* En las ecuaciones cubiertas en este estudio, e expresa la 
rata de intrusión de agua al yacimiento, dWe/dt. De acuerdo a la no- 
menclatura de la АГМЕ17, debe ser ey, pero para presentar las ecuacio 
nes más simplificadas se usa simplemente e. Cuando se trata deun sis 
ma de un yacimiento que fluye a un pozo, en lugar de e se usa q. En 
forma similar se usa u y К en lugar de y, y ky. El Significado de cap 
es; contacto agua-petróleo. 
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о para t > o. 
p(0,t) = constante 
apt0,tp) = constante 


b) Igual que la Condición b) del Caso 1 


II - Acuíferos Finitos (Limitados), хэй 


Caso 3 - Rata de intrusión constante en el límite interior (contac 
to agua-petróleo) y acuífero cerrado, es decir, rata de flujo en el 11- 
mite exterior igual a cero. Analíticamente, las condiciones de límite 


son: 


a) Igual a la Condición a) del Caso 1 


b) (e), = 0 para todo t 


3p 
) = 0 para todo t 
эх 
др 
— (6,5) = 0 
эх 
3p 
( ML,tp)= 0 
ar 
Caso в - Rata de intrusión constante en el límite interior (con- 


tacto agua-petróleo) y presión constante en el límite exterior. Analí 


ticamente, las condiciones de límite son: 


a) Igual a la Condición a) del Caso 1 
b) p = p, para x=L y para todo t 
p(L,t) = pi 


p(b,tp) = р; 
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Caso 5 - Presión constante en el límite interior (contacto agua- 
petróleo) y acuífero cerrado, es decir, rata de flujo en el límite ex 


terior igual a cero. Analíticamente, las condiciones de límite son; 


a) Igual a la Condición a) del Caso 2 


b) Igual a la Condición b) del Caso 3 


Caso 6 - Presión constante en el límite interior (contacto agua- 
petróleo), lo mismo que en el límite exteríor. Analíticamente, las 


condiciones de límite son: 


a) Igual a la Condición a) del Caso 2 


b) Igual а la Condición b) del caso 8, 


Varios autores 7, 8, 9, 19, 20, 21 han estudiado uno o varias de los 
casos descritos. Entre ellos, Nabor y Barham? presentan en forma my 
práctica soluciones de la Ec.(15) para los seis casos considerados en 


acuíferos homogéneos. Tales soluciones se presentan a continuación. 


Caso 1: : 
ev kt 
Ар = (ру-р) = 2 (26) 
kbh "$uc 
Caso 2: 
kt 
We = фбће(Ар)|2 — (27) 
чфис 
Caso 3: 


ej L kt 2L 1 2 
Ар = (рї-р) = (— + you e $ucL? 


donde n = 1, 2, 3, 8, ... 9 


Caso 5: 


2r2kt 
ey з * 1 DU. ZA : 
àp = (pi-p) = ь - i e duel (29) 
kbh x? n n? 
donde n es impar, n = 1, 3, 5, »». œ 
Caso 5: 
n?s?kt 
81 ^ 1 Ру 
We = ¿bhe(áp)¡L ~- 2 e puc (30) 
12 n p? 
donde n es impar, n = 1, 3, 5, ... 9 
Caso 6: 
n?s?kt 
з L kt 2L o 1 т 
We = фъһс(лр)|(—— + ——) - г. е euet? | (31 
3 фис т? п n? 


donde n = 1, 2, 3, 8, +... œ 


Las ecuaciones anteriores pueden simplificarse bastante, si 


se introduce el término tiempo adimensional, definido por, 


0,006323kt 


t. = 
D 
$ucL? 


donde los términos se expresan en unidades prácticas: k=md; t=días; 


p=fracción; u=cp; с = vol/vol/lpc; L=pies. El término Ь, para acuí 
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ros infinitos, puede tomarse un valor cualesquiera. Podría tomarse un 
valor imitario para simplificar los cálculos. En el caso de acuíferos 
finitos, (limitados), el valor de L es la longitud del acuífero. Ade- 


más si se definen los siguientes factores Fx (tp) donde k = 0, 1/2, 1» 


п?л?%һ 
BO. 1 jn 
Бо (Ер) = 1 - І O ч 2 
p 12 n n? Мо (32) 
donde n = 1, 3, 5 +... = 
tp 
Е1/2'%р) = 2 | —— (33) 
т 
242 
(ty) ( : ) : P : шшр: (34) 
F,(t s (t, + -—)- — 
У p 3 2 n n? 
donde n = 1, 2, 3, 8, +... = у 
Las Ecs.(26) а (31), en unidades prácticas, pueden escribirse, 
, 
Caso 1: 
887,58 eL 
Ap = (pi-p) = سس٣ )ےر‎ tp) (26a) 
kbh / У 
Саво 2: 
$bhLc(Ap) 
Ме = ——————-F (tp) (27а) 
2 
5,615 1/ p ' 
Caso 3: 
887,58 euL 
Ар = (pi-p) = ———Р. (tp) (28а) 


5 kbh 


~ 


Caso ц: 
887,58 ер 
Ap = (pi-p) = ————————Fe(ty) (29a) 
kbh ; 
Caso 5: 
ébhLc(Ap) 
Ме = 191437) (30а) 
5,615 
Сазо 6: 
$bhlc(Ap) А 
Че = F, (tp) (31a) 
5,615 3820 


Es interesante hacer un análisis de las funciones ЕС) Без. 
(32) а (34). Cuando tp tiende a cero (to » 0), la Sumatopía de la Ec. 
(32) tiende a la suma de 1/n? (n = 1, 3, 5, ... e), cuyo valor es 12/8 
Por tanto, el valor de F¿(tp) = 0, cuando tp + 0. Cosa similar ocurre 


con la Ec.(34), ya que la sumatoria de 1/n?(n = 1, 2, 3, ... ees, n2/6 


es decir, Fi(tp) o, cuando tp * 0. La Ec.(33), de hecho es cero, cua 


40 tp > O. 


De otro lado, cuando ty + =, F (tp) > 1,0 y Fy (tp) + (tg + 4 
es decir, llegan a ser lineales con tp, lo que es de esperarse si se 
analizan las condiciones de límite de los acuíferos finitos. Por ejen 
plo, para el Caso 9, la caída de presión llega a ser lineal con tiempo 
tan pronto se alcance flujo semi-continuo, o sea para valores de tiemp 
(Fig.3), después que el disturbio debido a la intrusión, llegue al lín 


te exterior. 


Para valores pequeños de tp, Ро (Тр) у 21 (tp? se acercan y lle 
gan a ser iguales al valor de F4/5(ty), tal como es de esperarse, ya 


que sistemas limitados se comportan como infinitos para valores peque- 


279955 
ños de їр. Tal cosa puede observarse en la Tabla I, donde se tabulan 
los valores de Py (tp) para tp entre 0,01 y 9,0. Además la Fig.10, рге 
senta gráficamente los valores de Fk(tp) para tp entre 0,02 y 10,0. 
De esta figura, puede observarse que el valor de las funciones Fo(tp) 
у F4(tp) comienzan a desviarse del valor de la función F1/2 to)» para 
tiempos adimensionales mayores de tp = 0,25, En general, pueden usar- 


se las siguientes reglas: 


1) si tp 5 0,25: 


Ро = F4/2 = Fi = 


2) Si 0,25 < tp < 2,5: 


Usar Tabla I o Fig.10 para F, y Fi. Para P,,5 usar la Ес.(33), 


3) Si tp > 2,5: 


= 
n 


1,0 


Ру = tp 


pi 
3 
Las ecuaciones presentandas para acuíferos lineales, se apli 
can bien sea para e o Ap constantes según el caso. El valor de e es 
la rata de intrusión de agua en el contacto agua-petróleo y en algunos 
casos puede considerarse igual a la rata de producción de petróleo, q) 
en BY/D, especialmente en acuíferos muy activos y cuando la presión del 
campo aún está por encima del punto de burbujeo. Si el valor de e о 
(Ap) varían con tiempo, es necesario aplicar el principio de superposi 
ción de efectos, el que discutirá más adelante y donde se presentarán 


las ecuaciones correspondientes. 
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Acuíferos Lineales Мо-Нопорблеоз. 


Mueller?! presenta soluciones de la ecuación de difusividad pa 
ra acuíferos lineales con variación de las siguientes propiedades: es- 
pesor, h; permeabilidad/viscosidad, k/u y (porosidad)x(compresibilidad), 
éc. Considera que estos factores varían linealmente con distancia a 
partir del límite interior del acuífero (exterior del yacimiento) has- 
ta el límite exterior del mismo. En base a esto, define el factor be- 
ta, 8, en la siguiente forma, 

у Valor del mismo factor variante en el límite interior 

Los resultados para varios casos de acuíferos lineales finitos 
se presentan gráficamente en las Figs.11 а 16. La Fig.11 presenta la 
solución para las condiciones de límite del Caso 3,(Fig.9), variando el 
espesor de acuerdo a los valores de 8 entre 0,25 y 4.0. Esta figura se 
usa para obtener el factor Е (р) de la Ec.(28a) (definido en este caso 
como Ар(Ер) en función de tp)» La Fig.12 presenta la solución para las 
condiciones de límite del Caso 4,(Fig.9) variando el espesor de acuerdo 
а los valores de 8 entre 0,25 y 8,0. Esta figura se usa para obtener € 
factor Foltp? de la Ec.(29a) (definido también como Ap(tp)) en función 
de ty. Las Figs. 13, 14 y 15 presentan la solución para las condicio- 
nes de límite del Caso 5,(Fig.9), variando respectivamente el espesor, 
el factor k/u у el factor фс, de acuerdo a los valores de B, entre 0,25 
y 4,0 para el espesor y entre 0,1 y 10,0 para los factores к/р y $c. 
Estas figuras 20 usan para obtener el factor Foltp) de la Ec.(30a)(de- 
finido también como QC) en función de tp. Finalmente, la Fig.16 


presenta la solución para las condiciones de límite del Caso 6, (Fig. 9), 
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variando el espesor de acuerdo al valor de 8 entre 0,25 y 8,0. Esta 
figura se usa para obtener el factor г: (р) de la Ес.(31а) (definido- 


también como Q(tp) en función de tp 


Los resultados descritos, fueron obtenidos por medio de compu 
tadores digitales de alta velocidad usando una solución numérica (dife 
rencias finitas) de la ecuación de difusividad para acuíferos lineales 
Debe tenerse en cuenta que al variar k/u, puede simplemente variarse in 
dividualmente la permeabilidad, la viscosidad o ambos a la vez. Igual 
sucede con el producto фс, es decir, puede variarse la porosidad o la 
compresibilidad individualemente о ambos a la vez. Así por ejemplo, ob 
sérvese el comportamiento de la Fig.1%, de la cual puede decirse que 
para todas las posibles combinaciones de permeabilidad en que varía B, 
a determinado tiempo adimensional, se obtiene el mismo valor de la in~ 
trusión adimensional. Esto puede justificarse si se tiene en cuenta 
que su cambio en permeabilidad no varía el volumen total de agua dispo 
nible para intrusión; únicamente desplaza el eje del tiempo cuandc el 
agua fluye del acuífero. También puede observarse que un sistema con 
B = 0,1, es decir, un sistema que disminuye 10 veces su permeabilidad 
hacia el límite exterior, la diferencia máxima de la intrusión adimen- 
sional respecto a la intrusión para un sistema homogéneo (8 = 1,0), es 
de un 25 por ciento, valor relativamente pequeño si se compara con el 
cambio es permeabilidad. Para un 8 = 10, es decir, el inverso del ca- 
so anterior, o sea que la permeabilidad aumenta a partir del límite in 
terior, el camhio de la intrusión adimensional máxima respecto а un a- 


cuífero homogéneo es del 33 por ciento. 
Acuíferos Radiales Homogéneos. 


Este tipo se refiere a acuíferos marginales, donde el agua en 
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tra al yacimiento lateralmente y su movimiento puede aporximarse а flu 

jo radial. Las Figs.1 a * representan esquemáticamente el modelo de 

flujo más comün en este tipo de acuíferos. La solución más conocida pa 
ra este tipo de acuíferos es el modelo resuelto analíticamente por van 
Everdingen y Hurst?. El modelo se cumple tanto para petróleo que flu- 
ya radialmente a un pozo, como para agua que fluye también radialmentea 
un yacimiento según lo indica la Fig.17, partes a y b. El espesor del 
acuífero, h, es pequeño comparado con el radio del yacimiento R,. En am 
bos casos, se considera flujo horizontal radial y demás suposiciones de 


la Ec.(25), 


э2р 1 Эр фис Эр 
FAS ii s (25) 
3r? r дг k at 


Generalmente, con el fin de resolver esta ecuación más fa - 
cilmente, se transforma en otra ecuación en función de radio y tiempo 
adimensionales, rp y tp» definidos en la Fig.17 y además un término de 


presión adimensional, 


(pi-p) . (pi-p) 
: Р = Ap(rp,tp) = 2 E 2 
—— (35) 
2пкћ 2akh 


donde q y e son las ratas de producción al pozo y rata de intrusión al 


yacimiento. 


Introduciendo los términos грубу y P en la Ec.(35) resulta 
la siguiente ecuación adimensional para un sistema acuífero (o yacimien 
to) homogéneo, 

э2Р 1 ЭР ЭР 


+ = (36) 
атр тр PFD dtp 


En forma similar el caso de acuíferos lineales discutidos ante 
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Fig. 17 – Modelo radial para los cuales se ha resuelt 
cimiento que fluye aun pozo, (b) acuífero 
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0,006323 Ко! 


tp: 
Ø со Mori 


PARAMETROS 
ADIMENSIONALES : 


to с 0006323 Ки! 
Ø Cw Ши Rf 


o la ecuación Че difusividad: (а) yo- 
que fluye a un yacimiento. 
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riormente, existen soluciones de esta ecuación de acuerdo al tamaño del 
acuífero (o yacimiento) y de acuerdo a las condiciones iniciales y de 
límite. La Fig.18 ilustra las condiones de límite para los diferentes 
casos resueltos: dos para acuíferos (o yacimiento) infinitos y cuatro 
para finitos. Las condiones iniciales y de límite para estos casos 


son: 


Condiciones Iniciales. Tal como en acuíferos lineales, la condición 
inicial, tanto para acuíferos (о yacimientos) infinitos como finitos 
(limitados), considera una presión inicial constante, pi, al tiempo 


tz0 para todo valor de г, Analíticamente, 
Unidades Dimensionales Unidades Adímensionales 
р(г,0) = pi АР(гр,0) = Р(гр,0) = 0 


Condiciones de Límite (Contorno)- Las condiciones de límite o con- 
torno se refiere a las condiciones impuestas al modelo en el límite in 
terior (pozos, o contacto agua-petróleo) como en el límite exterior ra 


dio r, si el modelo es finito o г. > * para el caso de modelos infini- 


Los diferentes casos incluidos en la Fig.18 han sido amplia- 
mente estudiados 7,8,19,2,22,23,25,25,26,27, Entre estas referencías, 
una de las más conocidas para la solución de todos los casos, excepto 
el Caso 6, es el estudio de van Everdingen y Hurst’, vEH, resumida y 
ampliada luego Chatas?, Además, una sola solución general del Caso 


1 ha sido obtenido principalmente por Mortada?8, Theis?? 25, 


21 


y otros 
26 y el Caso б ha sido resuelto por Mueller А continuación se pre 


sentan las condiciones de límite y las respectivas soluciones para los 
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casos incluidos en la Fig. 18. 


1 Acufferos Infinitos - En este tipo de sistema, existen dos ca- 
sos: Rata Terminal Constante Presión Terminal Constante. Se han de- 


nominado Caso 1 y Caso ? respectivamente en este estudio. 


Caso 1 - Rata de intrusión constante en el contacto agua-petró 
саво 2 E = 
189 (acuffero-yacimiento) o límite interior. En el límite exterior, 

la presión tiende a pi cuando el radio tiende a infinito. Analítica - 


mente, las condiciones de límite son: 
Unidades Dimensionales Unidades Adimensionales 


a) e = constante para t > 0 


Aa P | para t 0 


.. 
у эг JRi 
8р E ey ЭР : 
(81,0) = — (1,5) = -1 
эг 27khRj arp 
b) р > pi cuando г + = para t > 0 
p(r,t) + pi, cuando r ¥ ® pírp»tp) + о cuando гр + = 
Existen varias soluciones para este caso. Las más comunes 


son: (A) solución exacta de van Everdingen y Hurst?, que permite cal 
cular la presión en el contacto agua-petróleo (o er el pozo), o sea 
а un radio Rj o сү, (o en unidades adimensionales a un radio ry = 1) 
y a cualquier tiempo t; (B) solución del Integral Exponencial, Ei, 
о de Theis??, que permite calcular la presión a cualquier radio r o 
rpya cualquier tiempo t, y (C) solución de Mortada??, que también 


permite calcular la presión a cualquier radio r o rpya cualquier 


"заимрол (OJBJNIO —OJUIWMNIDA) SOWIISIS озод POPIAISMJIP ap.uomonoa ој эр sauonnjos-gi $14 


02)-- 9 214045405 4015214 
7 кшз э14045405 4015214 
(2) (а) 2) — s. ороцәз озаџпоу 
IN СМЕ 
(22) (802) — У 214045402 9915324 
314015405 3 
оп}; эр оцон 
2) (8)(2) — € ороллаз 013jIn2y 
531У!дүн 
0:а))пзү - ойзатшиэод 
$9 3151$ 
(22) (8) (2) — г 319045402 uptsaJg 
90015002 
(8)(2) | | оо 
| 30138153 2011910] 
ЕРЕ] 0505 эиши әр seuonpuo) озуу старой 


тәр оџошој 


- 53 - 


tiempo t. А continuación se presentan estas soluciones en unidades 


prácticas: 


A) Solución Exacta de Van Everdingon y Hurst, vEH. 


141,2 eu 
p(R4,t) = pi - ар (tp) (37) 
xh 
donde y 
-u? d 
x. [* о" 
Ар(5у) = Ap(1,tp) = — سس‎ (38) 


12 5 и? (2200) + pa} 


Para la obtención de Apltp? pueden usarse las siguientes e- 


cuaciones simplificadas de acuerdo al valor de ty! 


Si ty « 0,01: 


t 
ap(ty) = 2 — а 1,1285 t (39) 
т | р 
51 ty > 100: 
1 
Ap(ty) = — (1n(ty) + 0,80907) 
2 
= 0,5 1n(tp) + 0,40454 (40) 


at 


1,15129 log (tp) + 0,40454 


Si tp > 0,01: usar la Ec.(38). Afortunadamente esta ecua 
ción ha sido calculada como función de tp- La Tabla 11 presenta valo 
res de Ap(tp) para ty entre оу 5000 y la Fig.19 presenta los mismos 


valores para ty entre 0,01 y 10000. 


Edwardson y otros??, presentan las siguientes ecuaciones 
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TABLA II - Valores del Cambio (caída) de 
Presión Adimensional, Ар(їр), para el ca 
80 1A: Acuífero (Yacimiento) Infinito, Ва. 
ta Terminal Constante, Ес,(38). (Referen- 
cias 7 y 8). Fig. 19. 


гтэ. 


tp Ap(ty) tp Ap(ty) 

0 0 20 1,9595 
0,005 0,025 25 2,0621 
0,001 0,035191 30 2,1469 
0,002 0,089587 35 2,2192 
0,003 0,060359 40 2,2822 
0,004 0,069433 45 2,8381 
0,005 0,077384 50 2,3883 
0,006 0,084528 55 2,4339 
0,007 0,091063, 60 2,4757 
0,008 0,097115 65 2,5142 
0,009 0,10277--- 70 2,5499 
0,01 0,10810 75 2,5833 
0,015 0,1312 80 2,6145 
0,020 0,15031 ‚85 2,6439 
0,025 0,1669 90 2,6717 
0,03 0,18176 .95 2,6979 
0,04 0,20767 100 2,7229 
0,05 0,23007 150 2,9212 
0,06 0,24997 200 3,0626 
0,07 0,26800 250 3,1726 
0,08 0,285855 300 3,2629 
0,09 0,29989 250 3,3389 
0,1 0,31423 400 3,4051 
0,15 0,3750 550 3,1635 
0,2 0.82425 500 3,5158 
0,3 0,50282 550 3,5632 
0,4 0,56462. 600 3,6065 
0,5 0,61686 650 3,6463 
0,6 0,66220 700 3,6833 
0,7 0,70242 750 3,7176 
0,8 0,73867 800 3,7497 
0,9 0,77172 850 3,7799 
1,0 0,80215 900 3,8083 
2,0 1,0222 950 3,8352 
3,0 1,1663 1000 3,8608 
4,0 1,2749 1500 4,0627 
5,0 1,3625 . 2000 4,2062 
6,0 1,5360 2500 4,3175 
7,0 1,4996 | 3000 4,4086 
8,0 1,5556 -3500 4,4855 
9,0 1,6056 4000 4,5522 
10,0 1,6509 .. 4500  . 41,6410 
15,0 .1,8293 ` 5000 4,6637 
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para el cálculo de Ap(t,) como función de tp: 
Si 0,01 < tp < 500: 


370,529 tp + 137,582tp + 5,69549tp \, 
Ap(ty) = ا لے‎  ——— 
328,834 + 265,488 Vtp + 45,2157tp + to Y р 


Máximo error respecto al valor exacto: 0,08 por ciento 
si tp > 500: 


1 1 
(г 


1 
Ар(тр) = E In(ty) + 0,808558 |1 + 
2 2t ut 
р р 


Máximo error respecto al valor exacto: 0,01 por ciento. 


Nótese аце el primer término entre paréntesis de esta ecuación 
es igual а la Ec.(40). Las Ecs.(51) y (42) pueden ser de gran utilidad 


cuando se trata de trabajos en computadores. 


B) Solución del Integral Exponencial, Ei. También se denomina solución 
de la linea fuente o solución de Theis?9, Aunque esta solución en 1а 
tecnología petrolera es más comunmente conocida como solución del inte 
gral exponencial o solución de Ei, parece que Theis fué el primero en 
demostrar como usar esta solución a analísis de flujo no-continuo. Par 
ticularmente en hidrología se conoce comunmente como solución de Theis. 
Esta solución permite calcular la presión a cualquier radio r o rp pa- 


ra cualquier tiempo, Al mismo tiempo puede emplearse para el caso de 


rp = 1 (r = Аро гу). 


La solución del Integral Exponencial introduce la siguiente 


definición de tiempo adimensional tor» 
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0,006323 kt tp (53) 


t = —— 5 


Dr 
$ucr? 


y 1a solución está dada por, 
181,2(e)g, и 


p(r,t) = pi - Ap (ny ty) (n9) 
kh 


donde, 


1 1 1 Pp? 
Aplrpsťtp) = с -Еі(- )| = | -Ei(- ) (45) 


A su vez, el valor del integral exponencial está ne? por, 
5 


eu 


-Б1(-х) = du (46) 


(3 
=} 
o también, en forma de una serie, 
Яо (am 
(47) 


Е1(-х) = ln(x) + 0,5772 + £ 


n=1 (n)(n!) 


Los valores de Ei se encuentran tabulados y graficados en el 


Apéndigce A, para diferentes valores del argumento X. 


si, 2 


usar la Ec.(40) para el calculo de Аре» ty). 


С) Solución de Mortada?9. Permite calcular la presión a cualquier ra 
dio r o ry para cualquier tiempo, usando una ecuación de la forma de 


la Ec.(44), donde el ap(rp ty) se obtiene en forma diferente. van Ё» 


‚ verdingen y Hurst presentan una ecuación de &p(rp,tg) del siguiente 


típo, 
-u2 
2 {= (1-679 6р) - (34 CU YoCurp) - Y, GDJ$S(urp))àu 
àp(rp,tp) = --| ыы (58 
5 u2[J2(u) + Y12(u)) 
1 1 
Sin embargo, no presentan una solución numérica de esta ecuación у 56- 


lo la usan como paso intermedio para la solución dada en el Caso Á, es 


decir, para гр = 1, la Ec.(48) se convierte еп la Ec.(38). 
р р 


Mortada?0, usando en parte consideraciones matemáticas y en 

parte computadores digitales y analógicos, presenta la solución de es- 
te caso. Las Figs.20, 21, y 22, incluyen valores de Ap(rp,tp) como fun 
ción de tp y para valores de гр = 1, 2, 4, 8, 16, 32, y 64. Edwardson 
y otros”, dan también valores de Ap(rp,tp) como función de tp para va 
lores de rp * 2, 3, 5, 8, 10, 15, 20, y 30. Tales resultados se inclu 
yen en la Tabla III. Finalmente, la Fig.23 presenta valores de Ap(rp, 
tp) como función de tp, = 0/25 » Para гр = 131,15 1,23 1,3 ; 1,» 


1,5 5 1,6 ; 1,7; 2,0 ; 3,0 ; 5,0 ; 10,0 ; 20,0, . 


En base a los resultados anteriores se presenta la Fig.2%, q 
relaciona гр “omo función de Ap(rp,tp) para valores de tp entre 1 y 
4000, Incluye también en esta figura, para comparación, los valores d 


cambio de presión adimensional, usando la solución Ei, Ec.(50). 


Mortada sugiere las siguientes ecuaciones para el cálculo de 


Ap(rp,tp), 


Si ty < 0,01, 


Ap(rp,tp) = Зее (49) 


Ver parte inferior de la página siguiente. 


59 


"000050 
z£000*0 
STI00'0 
T6900*0 
90810*0 
11150290 
19060" 0 
EHA 
99nz*0 
081870 
665490 
159550 
0552°0 
1658*0 
вепо* Т 
Бета Y 


og = ба 


21000"0 
11000290 
98500*0 
96810°0 
20920290 
26990"0 
петт?о 
tzoz“0 
928210 
89Th*0 
zhzs'0 
261990 
888120 
921670 
1660“1 
186211 
0628 T 
189891 


02 = M 


90000*0 
T5000*0 
oze00*0 
99600*0 
£8T€0'0 
02190"0 
26060390 
881170 
264210 
114650 
908890 
8089 “0 
42810 
9088 90 
6t50* T 
*"6LlT 
921691 
LTTS*T 
560131 
015891 


ст = Пл 


10000“ 0 
52000" 0 
toroo*o 
£0900'0 
6,610*0 
STSh0'0 
18080“0 
245170 
25422170 
888270 
ШААМ) 
162590 
626999 
#918°0 
2866"0 
IETT 
"96241 
£0£h*I 
169551 
649151 
5606" 
660152 
9252“ 


от = la 


є9200*0 

20000250 Ls600*0 

ST000*0 £L6T0'0 

14190050 senno“o 

Or500*0  ZTZLO'O 

10000'0 s9800%0 8£660*0 

"1000" 0 18452020 HETO 

LLOOD*O 241ч0"0 nozz*o 

17000050  19500*0 500т*0 691670 
10000950 Е2210"0 А) 096Е*0 
88000*0 895800 erecto 112550 
186000 061900 14.06 “0 00с9*0 
L8L00*0 £1960'0 161650 TIOL “O 
09920*0 тело 0805 *0 92580 
90250*0 тъне“ 0 ЄЕТ9 0 151650 
оотт'о engeto 26140 288191 

289190 nc95*0 0n68'0 SSLC* T 

9842°0 #519°0 9010*t 0095*1 

s£9£*0 826120 00451 666591 

sono ч 12850 560691 1669*T1 

9165 *0 8900*t 0н68“1 1268“1 

АТМ 7298: £TE9'T 116052 

L888*'0 «826191 21281 є62242 

181051 8S9n*T *"L96' Tl ТАА 

ZLOZ*T 1189"1 ч99т1 2 2045 2 

¿Ene T тёбе*т єв06“2 SzIL'& 

60481 511691 LBI} Z teze“z 

148951 арыг ESE Suz0*£ 

S884°T LESZ'C ГАТА А 149176 

188691 055822 sh96*z 169656 

ч0єТ72 08652 13411385 161596 

ЅТЕЕ Е соов“ г 201656 SSTL'E 

CETA DEA "65672 BESRE 165826 

в = бї $ = $ = a z da 


“(14 PID2US42J94)'91uP3SUu0) триүшдә >з 


*o3rurjur 


(озцәтштоед o) odeirnoy :91 ose) ‘la эр задотед зэзиэлэзта Laed [euorsueurpy одшатј әр 
ugroaunj ошоо * (04 a)dy *reucrsueurpy чотѕәла ep (терүгә) отҷшеә тэр зэлотел - III VIEVL 


-60- 


— 


оро /ojua 
101рол ошај 5! 5 


En 


4 


(92 отоџалајон) «ОТО 


[DUOISUAWIPO одшэ! эр ч 


opun owo? (6; 


оу dy 


¿01=Y * 211 -0) "орла 204 әр 
20Á үэр 04492 {әр AIDA D OIPDI= 0: *SjUDFSUOI (UOISNIJUL 9p)jouruJaj оцон оноо! 


0 * dvNOISN3YlIQY OSWAIL 
су 4. 


2 


100 2 су $ 2 1990 


¡DUOISUDUIPD upise1d op огдшоо 


i 


-02 by 


$ 5 
бу NOISZUA 34 OIBAVO 


5 
Ф 


зїй 


(буда) су “IVNOISN 


- 492 0190913404) (60101149) 69 0 t =0 ‘очаниюок Jap 
о1роз /ojuonut2oK јер 03439 јор 1104 0 0!рол=6 "91013405 (80:58:40) Эр )јошшла DOY ‘O! 1 
101р0л 24194515 Y (оцол5иэширо ойшэц эр upiouny ошоэ (901) dy ° jpuorsueupo uprsa1d ap їдшо2-12 ©з 


9,* зүмог5мЗунау Оаа 
4 8 v € 2 4 $ v € 1 


одот 2 soe 2 oot 


-6р- 


"Qv м0(53н4 30 OIGWVO 


~ 


(690) dy “IVXOISN 


-62-- 


үәр/о:ро1 орәшооК үер 024495 
10!ро2 024156} 


{ap и} 209 o о!ро 
jouoisuawpo odi 


4 


X92 0198319498) 
Y 9¡UOSUO) (UPISNAJU! p) Ouwa} 01033 
ipuoisuaupo 4915229 эр 


prouny ошоо (d 


qVNOISN3WIOY OdW3i 
t 


[2 


-22 b 


(°) dy “IVNOISNZMIOW NOIS3Ud за О1841УЭ 


^ {с г Di9u019303) юнючөцодхэ jo1bojur үер орош А ороџоџу 
Ч 291) ¿samp А цэбирлелз цоггор opojau (ә 10d (0:рол ofni} ооа popi^isnjp эр шогоопоэ ор эр UPIMIOS —'ez 613 


[6923] 294 "IVNCISNZMIOW Онц 


ооо! 098 o'i 100 
E 7 = A ое ИИ ] LP : : 747. 1:20 


.. WION3NOdX3 1489 1N! 
130 могп705 


(®, 0)) dy “IVNOISNINIOV NOIS3Hd 30 ойу 


59) pouorsueunpo odws А 0/100:24201р9 01003 5р чо!оџп; осоо 119) 505) 41) "igocisuewipo 215829 ор CUED yT ЗА 


Md 39 O97 


esa 


(69)23'!3 NOMOS === 
012vX3 КОЗМА y 


GROSA HI 


10,140 


URNA СТЕ 


- 55 - 


Si tp2 500: 
1 v? 
АРЬЕ.) = ——|- Е1(- DT (50) 
D 2 yty 


Comparació Soluciones del Caso 1,- Mueller y Witherspoon?? pra 


sentan una comparación del valor del cambio de presión adimensional, 
. 

АрСоу» t5), usando la solución de van Everdingen y Hurst (ry=t), la sa- 

lución 40 Mortada y la solución del integral exponencial. Los resulta 


dos aparecen on la Fig.23. Como puede observarse, esta figura presen» 


ta Ap(sp,tp) como función de tpp, definido por ia Ec.(53), es docie, 


круга. Le solución de Mortada рага r,y=1,0 es igual a la solución de 


van Eviediugen у Hurst, Las curvas entro грз1,0 y r,):?0,0 fueren obte 


nidas de la solución 43 Hovtada, haciendo la transformación dcl tienpo 


adiüenalonal tp al ticapo adimensional врь,Ес. (43). Раза 205 radios 
no incluidos en el trabajo do Mortada, se obtuvo la solución por compu 
tadores, en forma similar a la técnica usada por Mueller an otro astu- 
31021, La curva infevríoe, en parte indicada соп linoas inteorrumpídas, 
ropvescata la solución del integral exponencial, obtenida asigaando vi 
loves а tpp y calculando Ap(ry,ty) de acuerdo a la Ec.(15). Es intnva- 
sante observar la diferencia en los valores de Ap(rp, tp) de acuerdo a 


si че usa la solución бе van Everdingen y Hurst, Movtada o o1 140 


E 
e 
R irvíalx) зо — ~ (x)enfc(x) 


т 


ro 


2 
evfala) = ~ | e-u du = 1 ~ ent(x) 
о 


Exponancial. A medida que disuinuye el valor de tpp, aumenta la dife 


soluciones de УЕН y Могсолд en сопрага- 


ción con la solugión БХ, La Fig.25 ilustra cl error са porcentaje, 99 


по función de tp, de la aif moia de las soluciones de УЕН (тр= 1 y 


top = tp? y del Ei, tonto хоп Т pacto a la solución do УЕН cono res- 


pucto a la solución БІ. Claramente puede observarse que a medida que 


a 
o 


disminuye бр sunenta ol error, de donde pus que la solu - 


ción del ini 


spueto a la so 


pal exposscacial pezsanta un yor согог 


lución de УЕН a medida que Tal ёгсог debe tenerse en 


ogpo2ialnmon^a va valores de tp рот азћајо da 10. 


-2032ní4 observar de la Fig.23, como а modida 
qua aumenta el valor de vy, da diferencia de Ар Orp> ty? entro las dis- 
tintas soluciones, oou wespecto a la solución Zi, es menor y puede con 


cálculo 


мега pA 


favo-yacimieuto), о sua за ol linite interior. Соло 2n 


el radio tiende а in 


el lisite exterior, 


camente, las 


s Ding 


= лодљЕтлаће, 


JACO 


constante 


pí1,ty) 


Un otra 


(^p) tnnte, para t » 0 


Ap(Rsg,t) вс. 


АРС ty) 2 constante 


5200127305 $5} auy N 40: 0:5494Иро обшэ;; ap ие 
| | ! 1p 


“отэчәчоахә 10:56: јар А }SANH А usDuipao^2 UDA эр роріл:ѕпир эр иотэолоа Oj ғр 
0:240) ошоо 10419 ар 9(10:492:404 ~ са ба 


Э Py co О 


Ши 


CO id MON GO 


gemuesvocxo збаци : 


say А зөбоитрзалл UDA ! 
ёр 0012105 pj 0 0422652; 26222 сү: 
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b) Igual а la Condición D) del Caso 1 


En esta caso, vau Everdiagen y Hurst” ргезезфай la solución pa 
va са1си1ах la cantidad total de agua que pasa а través del contacto a- 
gua-petróleo (o que fluye al pozo), manteniendo una presión, р, cons- 
tante en el contacto agua-potróleo (pozo) durante un tp. La solución 


en cantidades prácticas es, 


0 
Ma(R,t) = e s 1,119 фе h 82 -- OP QCty) (51 
360 
donde, E 
" 9 (1-е буди 


ыы (52 


008.) s 0(1,5,) a 5 
| ý 5 Јо u3 (1200) + 1200] 


Алмаз, 0 presenti 01 angulo en grados, corsespondicatoa al avso de 


cilvcunfapcaeia en que sl acuífero es activo al yaciniento. 
En forma similar al caso anterior, 
$i t, < 0,01: 


0055) = 


NE T 
= 2,1288 Ves 88. 


Si tp > 0,01: usar la Ec.(52). Afortunadamente cota ecuación 
ha sido calculada como función de tp. Los valoras 48 обе) ве presen- 


tan en la Tabla IV po 


а tp entra оу 2х1032 у en las Figs.26 у 27 рака 


valores de Хү entre 0,01 y 1010, 


Liwacdson y orros?? аваз ecmelones para 


а) calaqulo de Q(tp) сл func 


Si 0,01 < to < 200: 
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TABLA IV - Valores de Intrusión Adimensional, Q(tp), como Función 
de Tiempo Adimensional, Үр, Caso 2: Acuífero (Yacimiento) Infinito, 
Presión Terminal Constante, Ec. (52) » (Referencia 7 y 8). Figs. 26 
y 27. 


to об) їр (tp) — fp Уб fa 21) 
0,00 0,000 33 18,011 79 35,697 223 83,497 
0,01 0,112 34 19,429 80 35,058 230 85.023 
0,05 0,278 35 18,245 81 36,418 235 85,555 
0,10 0,404 36 19,259 82 20,777 240 88,052 
0,15 0,520 37 19,671 93 37,136 245 89,575 
0,20 0,606 38 20,080 8^ 37,895 250 91,08" 
0,25 0,689 39 20,388 85 27,851 255 92,589 
0,30 0,758 чо 20,89% 85 38,207 260 94,090 
0,40 0,898 ui 21,298 37 38,563 265 95,588 
0,50 1,920 52 21,701 88 38,919 270 47,981 
0,60 1,150 43 22,101 89 39,272 275 98,571 
0,70 1,251 цц 22,500 20 39,626 280 100,057 
0,90 1,859 45 22,897 91 29,979 285 101,540 
0,99 1,569 46 23,291 92 52,331 292 103,019 
1 1,569 47 23,08% 93 40,684 295 104,995 
2 2,147 ^a 24,076 94 41,035 $60 105,968 
3 3,202 49 24,165 95 41,385 305 107,537 
u 3,893 50 24,855 96 41,735 310 108,90 
5 1,5392 51 25,298 97 12,084 315 110,367 
6 5,143 52 25,633 58 42,33 320 111,827 
7 5,743 53 26,020 99 u2,781 325 n 113,284 
8 6,314 54 25,405 100 33,129 330 114,738 
9 6,869 55 26,791 105 44,358 335 116,189 
10 7,411 56 27,174 110 46,574 340 117,638 
11 7,940 57 27,555 115 48,277 355 119,083 
12 8,457 58 27,935 120 59,968 350 120,526 
i3 8,964 59 28,315 125 51,648 355 171,966 
14 9,461 60 28,691 120 53,317 360 123,403 
15 9,949 61 29,068 135 54,976 355 121,938 
16 10,:38 62 29,443 150 56,525 370 125,270 
17 10,913 63 29,818 155 58,265 375 127,699 
18 11,386 64 30,192 159 59,895 380 129,125 
19 11,855 65 39,565 155 61,517 335 130,550 
20 12,319 66 30,937 160 63,131 399 131,972 
21 12,778 67 31,308 163 64,737 395 133,391 
22 13,233 68 31,679 i72 65,336 190 135,808 
23 13,68u 69 32,048 175 67,928 1405 136,223 
24 14,131 70 32,417 180 59,512 410 137,635 
25 14,572 71 32,785 185 71,090 415 139,055 
29 15,013 72 33,151 190 72,661 420 140,453 
27 15,550 73 32,517 195 74,226 425 141,859 
28 15,383 jM 33,882 200 75,785 "30 103,262 
29 16,313 75 34,247 205 77,338 135 1:44 , 664 
30 16,742 76 34,611 219 78,895 450 146,064 
31 17,167 77 34,974 215 80,128 HHS 147,461 


32 17,590 78 35,336 220 81,965 450 118,856 


TABLA IV - (Continuación). 


АГ 90) to 
455 150,249 25 
460 151,640 830 
155 153,029 850 
470 154,416 850 
475 155,801 860 
430 157,184 870 
485 158,565 875 
„90 159,9%5 880 
495 161,322 890 
500 162,698 900 
510 165,444 910 
520 168,183 920 
525 169,549 925 
530 170,914 930 
540 173,639 940 
550 176,357 950 
560 179,069 960 
570 181,778 970 
575 183,124 975 
59 134,473 980 
590 187,166 990 
600 189,852 1000 
610 192,533 1010 
520 195,208 19020 
625 196,544 — 1025 
630 197,878 1030 
сњо 200,542 1040 
650 203,201 1050 
660 205,854 1060 
570 208,502 1070 
675 209,825 1075 
680 211,185 1080 
690 213,78% 1090 
700 215,817 1100 
710 219,046 1110 
720 221,670 1120 
725 222,980 1125 
730 224,289 1130 
780 226,008 1190 
750 229,515 1150 
760 $232,120 1160 
770 235,721 1170 
775 236,020 1175 
780 237,318 1180 
790 239,912 1190 
800 242,501 1200 
810 245,088 1210 
820 247,668 1220 


QCty) 


248,957 
250,255 
252,019 
255,320 
257,953 
260,515 
261,795 


- 263,073 


265,629 
268,181 
270,729 
273,275 
274,545 
275,815 
278,353 
230,988 
283,420 
285,948 
287,211 
288,073 
290,995 
293,514 
296,030 
298,543 
299,799 
301,053 
303,550 
306,055 
308,567 
311,066 
312,311 
313,562 
316,055 
318,545 
321,032 
323,517 
325,760 
326,000 
329,480 
330,958 
333,133 
335,906 
337,152 
338,376 
310,843 
343,308 
345,770 
318,230 


1225 
1230 
1240 
1250 
1260 
1270 
1275 
1280 
1290 
1300 
1310 
1320 
1325 
1330 
1350 
1350 
1360 
1370 
1375 
1380 
1390 
1800 
1810 
1420 
1925 
1830 
1980 
1450 
1450 
1570 
1975 
18480 
31490 
1500 
1525 
1550 
1575 
1600 


1800 
1825 
1850 


349,460 
350,688 
353,144 
355,597 
358,048 
360,496 
361,720 
352,942 
355,386 
367,823 
370,257 
372,709 
373,922 
378,139 
377,872 
380,003 
352,132 
384,859 
386,070 
387,283 
389,705 
392,125 
394,543 
355,959 
398,167 
399,373 
401,786 
00,197 
106,606 
409,013 
410,214 
411,418 
113,820 
116,220 
422,214 
528,195 
434,163 
140,128 
415,077 
452,015 
457,945 
463,853 
169,771 
575,669 
181,558 
487,427 
493,307 
459,157 


3250 
3399 
3350 
3400 
3450 
3500 
3550 
3500 


522,520 
528,337 
534,145 
539,945 
555,737 
551,522 
557,299 
563,068 
568,930 
573,585 
580,332 
586,072 
591,806 
597,532 
603,252 
608,965 
614,672 
620,372 
526,068 


.631,755 


637,437 
643,113 
648,781 
260,093 
671,379 
682,640 
693,877 
705,090 
716,280 
727,449 
738,598 
799,725 
760,833 
711,922 
782,992 
794,042 - 
205,075 
816,020 
827,088 
838,067 
819,028 
859,974 
970,903 
801,816 
092,712 
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TABLA IV - (Continuación). 


Q(tp) Q(tp) tp Q(tp) tp Q(tp 
1 
3650 903,594 7100 1525,729 2,5410 »,066x10" 6,0x10  5,368x 
3700 914,459 7200 1546,011 3,0x у,817х " 1,0% © 6,220x 
3750 925,309 7300 1566,265 Y,0x " 6,257х " 8,0x " 7,956» 
3800 936,144 7400 1686,90 5,0x " 7,699x " 8,0x "40 7,909» 
3650 946,966 7500 1706,688 6,0x " 9,113x "5 1,0x10 8,747» 
3900 957,773 7600 1726,859  7,0x " 1,051х10 1,5x " 1,288» 
3950 968,566 7700 177,002 8,0x " 1,189х " 2,0x " 1,637» 
1000 979,344 7800 7767,120 9,0x " 1,326х " 2,5х " 2,103» 
4050 590,108 7900 1787,212  1,0x109 1,462x " 3,0x " 2,505» 
4100 1000,858 8000 41807,278 1:,5:-" 2,126x " Y,0x " 3,299» 
4150 1011,595 8100 1327,319  2,0x " 2,781x " 5,0x " — 5,0873 
4200 1022,318 98200 1347,336 2,5%: 3,927х " 5,0х " 1,868» 
1250 1033,028 8300 1857,329 3,0x " 4,05x " 7,0x " 5,643; 
4300 1053,72 3400 1387,298 Y,0x " 5,313x " 8,0x " 6,818) 
14350 1054,409 8500 1907,243 5,0х " 6,58х " 9,0х " 7,183: 
4400 1065,082 8500 1927,166 6,0x " 7,761х " 1,0x10% 7,658) 
4450 1075,743 8700 1987,085  7,0x " 8,965х " 4385 1,17 à 
4500 1086,390 8809 1956,942 8,0x " 1,016x109 2,0x " 199.3 
4550 1097,02% 8900 1036,7393 8,0х ", 1,10nx " 2,5x " 1,92 > 
8600: 1107,645 9000 2006,528 1,0x107 1,252x " 3,0x " 2,29 1 
8650 1118,257 5100 2026,%38 1,5х-" 1,828x " 4,0x " 3,02 ; 
8700 1128,855 9200 2050,227 2,0x " 2,393х " 5,0х " 3,75 
1750 1139,439 9300 2055 ,996 2552. t 2,951x " 6,0x " 1,573 
1300 1150,012 95900 2085,75 — 3,0x " 3,517x " 30x № 5,19 
+850 1150,574 9500 2105,473 Кох " 4,610x " 8,0х " 54893 
h400 1171,125 9500 2125,164 5,0% " 5,689x " 9,0x ",, 6,58 
1950 1161,666 0700 218%,878 5,0x "  6,758x " 1,0510 7,28 : 
5000 1192,198 9800 2165,555 7,0x " 7,816х " 1,5x " 1,08 : 
5100 1213,222 9909 2185,215 8,0x " 3,266x " 2,0« " 1,12 | 
5200 123%,203 19000 2203,801 9,0x " 9,911x " 
$200 1255,151 12500 2608,967 1,0x108 1,095x107 


5400 1276,037 15000 3,564,780 1,5х." 1,60%х " 
5500 1296,043 17500 3633,368 2,0х " 2,108х " 
5600 1317,709 20000 095,800 2,5x " 2,607x " 
5700 1338,886 25000 5005,726 3,0x " 3,100x " 
5800 1369,225 30000 5899,508 4,0x "  4,071x " 
5900  1379,927 35000 6730,247 5,0x " 5,032x " 
6000 100,593 30000 7050,096 6,0х " 5,984x " 
6100 1821,22" 50000 13363,099 7,0x " 8,928x " 
6200 11441,820 60000 31047,299 8,0x " 7,865x " 
6300 162,383 70000 :2708,358 9,0х " 8,797x " 
6:00 1182,912 75000 13533,8457 1,0x10 9,725x " 
6500 1503,808 30000 19350,121 1,5x " 1,829х " 
$600 1523,872 57700 15975,2389 , 2,0x " 1,880x " 
6700 1544,305 100000 17586,2845 2,5 У 2,328х " 
8800 1554,706 125000 21560,732 3,0x " 2,771х " 
6900 1585,077 150000 2,538х109 4,0x " 3,645x " 
7000 1605,18 200000 23,308x " 5,0х " 4,510x " 


о 
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Fig. 26-Intrusih adimensionol, Q(t5), como función de tiempo odimensional, fy, Sisterao radial infinito. 


Presion(o coída de presicn) constonte,es decir, presion terminal constante. Intervalo de to: 
1072 o 10% (Referencios 7y 0). 


m 


ma rodial infinito. 
stante.Intervalo de fp: 
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Fig.27-Intrusión edimensione:, Qlty),como funció 
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> presión) con 


НА 


Presiónto cci 
10% o 106, 


орел 


| 


1,128238 X + 1,19328t, + 0,269872t ү 


Q(tg) 2 ——————— == 


1 + 0,616539 Мо МЕТА (: 


+ 0,0035529%(tp)2 


Máximo error, respecto al valov axacto, 0,02 рог ciento 


Si tp > 200: 
-в,29881 + 2,02566t 


QUEJA аана ата (5: 
1n(tp) 


Máximo error, respecto al valor exacto, 0,07 por ciento. 


II - Acuíferos Finitos (Limitados) - En este caso, г ticne un va- 


los máximo constante denominado го» correspondiente al lim'te exte- 
rior, El valor adimensional de este radio гс será гора pQ/Rg ( en 


el caso de yacimiento тор = 7/0). 


Existen cuatro soluciones para esto tipo de acuífovos: dos 
de Rata Terminal Constante y dos de Presión Terminal Constante. EET 
ban denominado Caso 3 y Caso ^ à las dos primeras y Caso 5 y Caso ба 


las dos últimas. А continuación se describen las condiciones de limi 


te y las soluciones para cada una de ostos casos: 


Саво 3 - Rata de intrusión constante en el límite interior 

, 

(contacto agua-petrólen) y acuífero cerrado, es decir, rata de flujo 
en el límite exterior igual a cero. Апа1ї їсаасатс las condiciones 


límite son: 


Unidad 


a) Igual a la Condición 3) del 
Caso 1 


b) (e) = 0, pava t > 0 
те 


En otra forma: 


ta o о 

os га 3 O para t > 

3p ОР 

зөө (бозї) AR (гар? tp? 0 
де дур 


Este caso tesbión fug resuelto por van Everdingen y Hurst. 
Presentan la Ec.(55) para ol cálculo de la presión a cualquier radio 
v y a cualquier tiempo t, 


151,2 (Daz u 
p(r,t) * pi - ——— ¿ppt р) (56) 


kh 


donde, 


(rap)? бер) 


(rap)? - 1 


dp (70:54) з 


р)" "(тор)" атар) - (ap)? - 1 


п те Bonne m db 


«| гойд” x 1)? 


[14 Can Yrs Cagnp? 74033270 Саа ep)) 


donde ap son las raices de la ecuación, 


J; Ca сыр! Таи) - 91 Cap) Үү (n r.p) = 9 (5i 


De la Ec.(37) no s^ hon dado soluciones ^n fora tabulada. 
Es prácticamente la misca saulución dida por ма кас2 3, 


El caso 3 también Ба sido vasualto рог УЕБИ7 para obtener la 


ба en el gostuko agni-potróleo (o on el pozo) а cvalquier tiempo 


141,2 кади и 
PERA menm pte) (59) 
kh 


donde, 


14 "tp ВН M татар) -2 (rap sí 
Ap(ty) = Apli, pp) + - uM —————— өнд 
al. 2-1] "T ез? 


eo -а2 


(59) 


ia 


donde лду Cono autos, son Май saicas da la £c.(53). ба с,(00) по us 
: 


ads qua la Ec.($7) pora гра! 


Si ty s 0,01: шил» la Ёс.(39) рога calendar Ap(ty). МЕН prascatan la 


soianión de la Ze.(90) рага rap entre 


ргошсайан an Ла Tabla V y la Uiz.2. esloves de сар се ha 
caleulaco Ар Ку) usando la Гавс100: УС рд) de ИВЛ, dszarvollada por 
Названа, Brous y па ађтози 3, de сепегдо а la siguioate ecuación у зас 


gûn se duduca en cl Арјадјев В, 


ME o) r tpa ~ РОСА) (51) 
"tp 


on basa al 


doade tpa es el 


Para un сїоач1о, la violación сике tp y tpa 9$; 


ty = try EDA (62) 


2 e 


Los resultados de los cálculos de la Ec.(61) para гор entre 20 y 25000 
па presentan también en Ja Tabla y y se han graficado, раса гор entre 


100 y 25000, en las Tigs.29 y 30. 


о - Rata de intrusión constante en el línite interior 


(contacto agua-petróleo) y presión constante en el límite exterior. 


Analiticamento, las condicionss da limite son: 


Ünidadas Dinensionates añes Adimensionales 


a) Imal а la Condición a) del 
Саво 4 


b) p = pi para r = р, y para 
todo t. 


píra,t) = pi PO, tp) 3 9 


Баба caso ha sido vosualto analiticamente por Hurst? considerando un 


„2+ 


onte у Musiat considerando una línea como fuente. 


cilindro 


En uubos casos, la solución permito calcular la presión ea el contacto 


“26100 (о en el pozo), 


1891,2(6)7 4 
рр) а рі = ——— 
kh 


Apíty) (63) 


donde: (1) Pava la solución de un ciliadro como Fuente: 


2 B 
© ”م‎ MACH rap) 
in( Pop у -2 £ مم ر پر پوق و‎ л ی‎ (64) 


n=1 2132 2 у 
afl (ав) A, ven?) 


áp(tg) 3 ^p(i,t, 


li ecuación, 


donde ад con 1 


Ју (ав) Y Са Хер» - Y (аа) Jalan 749) = 0 (65) 


(2) Pava la solución considerando una linea cono fuente: 


a 


Adi- 
Radial 

1 
Refe 


e igua 


~ 
о 
2 
a 
E 


or 


Ü 


o 
£u 
Ê 
Ф 
эн 
E 


de 
с 


a 
te Inter 
re 1,5 y 10. Ес.(60), 


n 
P 


u 
e 
Ec.(61), Figs 29 y 30. 


F 

1 f 

nte en el Li 
t 
1 Autor, 


> 
В 


9 
Ccnsta 
о). Valores de rap en 


s Valores Calculados por e 


Cerra 


{его 


ог (Acuí 
os demá 


i 
Пее 


а 
) 
te Exter 
Fig. 2 


017 
» 


4 
a3 
^n 
o 


AO tv 
ч юмо 
му чэ wn со 


a. om a... 


ел 


nonoa 
Пээ 


мэ 30 UT 


~ 
mun c 

WOOO O MOM MM O 
^ ^ ^ nn ^ ^ 


оороо ч 


жап цо суо о 
динамиках 


үз wi vl wt rl ка 


MOH © с A Фо о Оо 
e Qo mico со е 
Qt Cop ~ ~ са ouo 


ci» ж-а а и ди. ж A AA жол 


200220000000 veia 


(4 C4 CI CY ш (0 СУ 6 
тз Oo eb O I D M u 
FS De са (0. CI t0 o 62 


Фо Фо т Ww 


со 
ооло 
(o iD Ф tO 


oo 


o 
ДРОМ оу cina 


ољзосоосво 
оо тоо ос 
OSOAN 


стою 
rara ow 
Бум 
> 


ooooco 


o 
o 


тэгээ 
юном ~ «а с0 т ҮҮ OO foci 2 
лошолаасо ~ о готово 
nn. a. ^ тэгээ 
cOoOoOCOOo0o000000020202002020wn00 
ос оооооооо 
eo OI OM O O 7 O 
~ + an +. st n n 
oo O 0 Ф O v UN ом 


о сәг o 
оз CH iO лу QU 
сч OM m л 


эсооооососо сэ Ф 53 du 
d) O (Nr «о 00 (Чч t 10 о суп сз! ом O 
DARAS са Сч © с: 0 zt > 92-03 
оосоосооососооосоосооооО 


ж e сэ со со AN о Q0 ш сул чэ ei ~ DD 
24 | bO rt O O O CO nmi Үд 0) N 
с 5 хо о Ч се во СУ OY со оо OM O 
[57 ^ en monent «ns. 
" мээн 
a 
Ф оаа оо осал ошо O 1 O 10 O O © 
я я Mon n s m hm ^ я 
e о O «i e («ci e 0 Ф TOW d) 60 
плана edes ҮЭ ҮЗҮҮ ҮҮ et nt 
^A о сак я езщ (S е Fe x0 00 0 HO O O 9 (00 O 0 m0 C0 
Р му cO O сиса 60 10 с AA O £O D 4 
о | ~ ыз мо D t о о «О ге i5 D- CS O 00 O Ф n c (0 zr 10 00 
а. АК E A A 
и 4 Гаа HHH үү У чч Єч N MN N оч бэ 
a 
Ф ох O 10 O i0 ош ст 000000000000 
ER a an a НСВ ИО ^^ a ^m ^ mnn o 
+ со оо о О Ф 54 nuovo жож мз гч су үх сээр Qi ом 
HH HH e4 4 Та еч лч ча N N оста АГ 


ч e 
о одао о) 02:40 etu) (1: 00 Сэ tS са 1D OY 03 алал 00 et 
E 


Ї 
MIT гох 
5. E 
ЦАР: eed a ele а eb HH хэхэх 
a 
Ф сососошоюФоооооооосэт ооо 
= а ВИНУ ^ ^ nn ^ a eh аа ~ ~ 
Y БҮТЭЭЛ ооо 
"MEE сыч o 
mas Sou 
= x сїш Оо х о 
villanos Nano 
а кач о ^ e e е 
“|4 нана NANG 
A 
Ф onom ороо 
f a anne anan 
m ما صا چ‎ vono 
ANNO 
очи зове о EEE ш со an man 
0 с с O ci oO sine Мороко сова 
E A EN N сч N I OLI MM MO uoi г< со p Ch O тач 
SS LÀ хо» аара же ГК Жө гар ж, Ж 
n loa rt et aa үл лэн ne СЕУ 
a 
9 | neo ч coo 
A al м тала РЭГ 
m | со сэ mmo O] 
позе 


| оос сафо ж к лаге Сеш ОУ БЕШ т (oO (o in Oo m C 

Pd xao c co O rim тот 0 Ci if O 0 Оо OLN су озор Тт 
411 оооооячаачанаялмаянят Ж m oO Ф eb o ut 
~ аа па паа лао ЗЫР ЧЫ ЫР А so na na. o. 
< HHH HH O4 wi E TO r tf cf el vie мес MN еч сч 

" 

a о « C со cross OP OOM OO x ID сз O а ©› ИЗИ) 

Ф ata a ИУ ОВОГ ^^^ a m o n 

я Planas ч өм оч 06) 0 9 (з 69 со 00 лг тю MO scho d i0 


0S n m nb nd 


200 


Pep = 


100 


Tap 


eD 


аз мз U 0 ау илан шу мут ао ко 


КЕСЕГЕ мэ эз нь ии га Ч 
[EE суа я CO OF oF r 


Гэгээн 
ХА A үн А 04 e OT OL СЧ 


X X XXXXxXX 
соосчоомло 
со 


О Hen Na 


“a о ооо о о 
42 CACA bor wy X 
z өэ Gn 
D. na ona 
4 чат охо 
HHHH N N 

> ю 
ОБЕ: БЕБЕ ОБЕ кп ЕЕ 

= a 
Ol хххххихяхххххх ххх ухх хХ 
elónonooóooóooonuodnoooo 
ОО а а A ^ ^ h^ ^ n^ а 
з сч £4 ю шог о соо CD ord ed СЧ NOM то 


O бу M c4 DM ‘tb © zou цу P in uà 0 D 10 ма VO O 
iO мо ол uro ul YD uso 

о о э 0) 2 сч m д о оз О ct 29 со со 00 оз O со 03 со о 
< OAM FF додо о UD з? чо но у 00 1D Te OU Yl 
eon noon 


N m 
УЕ АУЕ E E 5E = Е? oct 
ті = 
Ө|химиххиххх XXMXXKXXXM 
2120030040000 Фф о © © © © OO 
жак mo EE E жол е ма а ае а 
оконная FIO O т о ти O e ооо 


a O оо. 


Гэгээ 
БО) лоо OY су эр e D m сч cao c СУ ана 
MN WM Gi со O OO © O HAN + чу O 0 D 


eben Coe ETES ЕЗБЕ 


00000000 


Ol mor or о сус мэ r есл 
= dp с та C1 со FF O 


БЕ 


о €» peo сз €) (S MD O wd FIO ау Ф O ver 2F 
И i ее р «0 eC ЧО O O E> e O з O 
e DE арт, са mana © eb еа oet Са оо m су т m О 
р 
. RA ET E E RPA Umm DAI а есет Жар 08 
Дд been domi ANN єч єч сч N OC © сом 9 
cu O can Фа © 000000 
(nn a. . . « nn. ^v^ n on ж а чу ө ж, 
„| ereren ion из n to DS 19 ло o (I OO m OG d O 6 0 0 O O 
Cl elei rf үл el cd vi vict eb ON са с ср CEOS (1 0 0) 3 ao 10 0 


соп‘ 95 и х0*6 
09H tz и X0'8 | 

146554 u х0*5 09H" 22 u хот: L9L'LZ 

Liv TZ u XO 085"02 u #08 LƏL ET aig А 
LL6' 8T a хоје 088581 и Х0'5 ¿2L ot Bu : 
1445 u *0°ё OS 9T а KOH LYE iT А B 5 
LLX* 8T „ XS*T O95 hT и X*0*£ 184551 De и atete E 
145521 gOTXCT 088561 u XS*Z 95:95 IAS ET и ЕЕ ^ 
154561 u *0°6 Сен‘ Ст u хет: LST HT iet n REED : 
165513 u X0'8 094911 а XST Loc ez 183001 u НЕ й 
11951 u XO*L OSa OT gOIXO't L3S ZT hp Herm Eu Жы 
165531 u хо'9 082501 „ X0'6 184551 758, : toto t M кулон 
145501 u хо 099701 x08 {560% EA а 3059 ETO! а KO'Z 
158501 a X0'h 093*6 L LST OT P ios “ ЖОР залж peo 
LLE* OT u XOE 08956 9 LeL' 6 ев n келе S IE Eno 
1110 и XS* 0986 s Leets 188 m kots 8665. а ®0°8 
410701 0926 u : 18618 un pp. ER v EUR 
86898 56056 и хосе 18598 P ud E ховс "iud оо 
22356 er6*8 u AS Z 18458 те: 8 Мо BET A 
0LS*6 02858 a хоса 10853 TOT ааыа Е Корс 
215'5 té9 °8 u XSST 81618 Lor pA 2 ios UN 
85856 94578 „0т®0*+ 92.8 094 шэн : 
21618 и хо'9 4258 u X0'6 18118 eesti geo эь 5 
18416 4 хој5 «9876 в хоје ero*s ЭЭС u х0°6 нэ : 
88156 а Xo +. 862308 u 50° 21814 отс. " ха леч c 
12055 u Х9“6 61218 ССРИ 2495. xo A por = 
81898  ,otmo%z 82118 j TLI L 8:659 | “015055 лаг: 

Gaig da 
0002 = де 


Y IG Y, 


48^4:522}24) (0p01:92 J0121X9 2. 01393X9 anut 19 чә 9125 0 1026: от} Эр 0402 А ALUDISUOO (ugisnajur) чоор 
«cad ар вЫ (оро) осин (9122) 924515 18:0:5с84Ирс обшец гр шоочу SO) [сос:зчәчлро 40:59:4 гр ссаш22-82 Pu 


2, 'HvNOISN3A:GY семь 
со Yo 


зма 30 01993 


мау NO 


лумс њи 


nág' 


y 


XNOISNZAIGV мо!538с 30 CAY? 


SIONAL, to 


м 


ПЕМРО ADIME 


toms sadial 
ово) y оне 


DES 


icr {contorto одем p? 


el Поло extorior iyual a cero. 


з 
2 
2 
5 
3 
2 
= 
ES 
2 


2 
32 
ЕЕ 
5 Е 
28 
== 
32 
"E 
Se 
58 
e$ 
26 
25 
БАС 
26 
Чи 
zS 
5 

H 

E 

E 

Е 

3 


limitaduífinito). Rata ds 


Cambio da pra: 


Fig. 29 


ratu de finjo en 


foro currado, es Чеси, 


Curva teg, 

$ 2000 
T 5000 
u 10000 


y 25000 


E 
TIEMPO ADIMENSIONAL, to 


jo tompo sets arion 
rior (cont. 


ensionel, Ар (0), como función 
onslente en of I 
tuju ea ol Minitu exterior igual a aro. 


a 45 nr 
nat 2 
cerrado, 05 decin cita 


2 
2 e етп J bp) 


Ap(tp) z Ap(1,tp) = 11(гор) Е —————————— (66) 
2 2 
гар n=1 bi Ji, Pop) 


donde bq son las raices de la ecuación, 
Jon rap) = 0 (67 


van Everdingen y Hurst" peescutan cúlculos de las Eos.(6h) у (66) como 


función de tp , para diferentes valores de бор. Para valores de т.р 
antne 1,5 y 50, se calculó següa la Ec.(64) y para кар entre 60 a 3000 
según la Ec.(66). Los resultados se presentan en la Tabla Vi y en las 


Figs.31 а 3%. 
gi t, < O,0i:usap la Ес.(39) para el cálculo de áp(tp). 


D 


si ty > (ор), ocurra flujo continuo, es decir, el valer de 


4215) es máximo y toma en valor constante e igual a, 
Ap (tp) = 1n (rap) (6f 


ión constante en el límite interior (contacto 


Caso 5- 


agua-petrólco) y acuífero cerrado, es decir, rata de flujo un el lf- 


mite exterior igual a сего, Analíticamente, las condiciones de limi 


a) Igual а la Condición а) del Caso 1 


b) Igual a la Condición b) del Caso 3 


La solución para ceste сабо, dada por УЕН7, permite cl cálcu 


lo del fluido total que pasa por el límite iufevice (contacto agua - 


en un tiempo dado t, 


"а 
Ф 
К] 
E 
10 
A 
о 
E 
a 
fa 
o 
E 
о 
о 
~ 

A 
m 


Ap 


Саво + 


1 


се 


Valores 


Presión Constanteen 


imiento) 


nterior y 


te 1 


Acuifero (Үгс 
mi 


Li 


el 


Constante en 


) 


Үер entre 1,5 y 3000. 


га 


D 


Pep 


л 


86 - 


CI EO IO f O сэ Єч О 10 O (CO Q (0 O 
©з їл C) тр оз NID O CO 0 i0 O NOD 
AMNION DOO O el e A N сч Сч 
ж“ na ^m a na от ол nt n n 


O O оо Ф ct et ct et noe on ed 


о € zc d 0 QO CF (0 9 O ct coin o 


nn. na л © ^ m m ко na 


мн NN NNN O со 0 AN 


0,981 
1,013 
1,041 
1,065 
1,087 
1,108 


оо сч о о 


1 
тэгээс 


HN MN CY CN CY 


м с о о смооюо 
ona сч ап со сузсо 
МО ~ ~ ~ DAORODO 


2000 оос хүү 


ослооосло 


ҮЭ Үү s ҮЗ үү Оч Оч см 


мол со 
оно 
зал ор 


nm O л © бугу 
Q9 Qu 00. 
10 9 оо со 0 


000000 


из ою о чу сш оч 
сэр 2F 0 UY ОО ~ ~ 


сз сэ >ос2ооооооооооо 
оо д) оо сч сч C 0 TODA AHO 0 
C) cf cb Ш ау оз O сч cf (D р t со со O 
сч 04 єч N сч CI со сб со о 00 €0 00 со 0 
оосоооооооооооо 
000000 

чэ ао £O f» оз Co. O WN zh O со O CF O 
0000060 v ct cd c хэ а ма 


осооососоосооосоосооосооо 


ососоосоооосоо 


ю «о C 0 O о мч доо 
mH HH ond ont 


€ со 0 10 (OL г О ANO 
NID O та O ANNANN 
vC с CIN OL N N N N сч 
nanna h^ n^n ^ л n 
md eb dod ond edo ocn! nod 4 
ососсоооо 


эээ 
ою то ог 


za) ~ ~ 
ооо 
2000 


eo 
EI ло 
oo 


Вечна 


мэ © ш со о о 


шо له‎ ~ 6 0 
= 


сз о OQ 10 O O со 10 O O оо 
02 OO O O c d e cien e 
OUND O O O N O OY O N O 
^en nnn бол т .. 
оооооооооосо 


сл оо сососооо 


^om m eon ^ n «ол єл 


сч єч єт сч со о оао ш О 


TONNO 
о (о DMN оо а ©: Ф Ф О 
«© tO о OOO O i9 (D (0 O 


© чэ сл соосу 


осоосооооососо 
onon 
бу € OQ Qo O t ctr WO 0 Оо, O 


ОФ Ф c c e et ei СИ 0 


со о Ооо 
о а со оос 
со с) оо 


2000000 


с сш Фм 0 
ХОГ 


2000200 


олоо тоо d «v DOM «Фал суту t e 
(p iO Coe D 60 v ЛБ D Qd C0 Q O oe on 
со о ома n (0 O фо хер eb NNO са 
оо ч с сс сч сч ст єч сз OL UL єз ҳо 00 оз e со 0) 0) 1 00 00 
OwowonwonounooooooooooooooooQg 
DD Q Qc t GOO00Oq:d000«0410q9€ounonooo 
ҮЭ eb хэ HNN NN (LO) со 3. 2F 0 tO 0. 
_ | са e io (S en p o cS р er (0 о о (0 n o O O o 10 © © 
S Фф m Od Oi Oz OO OO iuc oi moo 
Э x BON m eran ао S PS P Po C 0 цо ою © c) 0 OY ло ©, 
5. кызра QU da SA SAU IA нга AU н" Оки 
<] сс е еч сч єч ч ч сч са єч єч єч N N N су су N єч єч сч N N єч сч сч 
»" 
e 
9 
я оао сч со ос оз озо осш о осоо ооо о 
€» со 2 170 1 A QUO w озо шо со О тоо O 
Vnd ХА mH A са doe Сч NN OD HF 
“alos о HOM O N та св єч m 10 C4 er vi + O г> D 
Миг. СР КО ОУ СЭ С? оссо 
тээг. он 
ye 17 Е тер аера е ат ач ы Ж час E 
' 4 HNNNNNNNNN NN NC NN N OI (N N Cé N ээ 
" 
A 
» | 
ч срослооос 0000500000 
ENANA со бэс осо CO чао о о 
l e NN AN 
^ чооон чоно тоо опо ог O rt єч ч Oi O 
e постасосмолнпоочлосо~-веосфтоосово 
о ху iQ te (- о оул О «д «а чч 0 СЧ сч NN С сч NN O (0 00 02 #0 
ха ак nm UR алиа A кеи а е 
9 1 ера oe! ка (LG CL сч са сч сч NNN с 6ч Ч СУ ме 
1 
a 
Ф 
я a ONG о а ош ош соз соо Ou DID O O O O O OQ 9 O 
в | AS A A HAN ИЕ К-ДА P D (0 0 O er O1 99 3 10 
eo үү үү 9 
"да | өх 
jou 
m ~ о 
7 ~ ~ 
a «t Ка) 
сэ 
ue 5 сочоссоооососооооосооосооосооооо 
H ы Co( «о оз ду OY OO єч 24 оосчмоновомнодоо 
“аза лама CL C (0 00 о Фа о © 
| помзомсово“ ва доо 
брге сак O 0 ч 0 Hom o о ~ Фаз озб ал OD 
«о sz MMMM iF abc d an г. pop C CS ~ 
" AA e nl vd nion ci ond ont EE хэ 
с 
9 ош Ow сиз © мэ © мз 000000000000 
e под чэ uy (QD (M о о CD ON HF O мо С ОЧОР O 0) COO O о 
qoe 4 с CN CN СЕ OL 00 0 


509 


inuación). 


VI (Co 


TABLA 


800 
900 


1 
ч 

‚580 
ч 
e 


ooo o 
2) aloooooooooooo 
кросу ته‎ 


e idc 


o Pot co ex ro e m ОТО) 
Qo: l2 | No ор Эр о ~ 0 со 


я га 


са 


о 
я о] ос доо у суо сш) ооо 
из 


H јача су сч тога inco 


^) азо ef бо QN en D 054 2 


~ |оопгеосотиноо сео г с 
о S jou roma салга p сч 
о | > |јоостогосоомнса се 
o. M^ m ж hie m коб nnn n 
n 4 јоугтвоталла 4 
a 
Ф oo oooooooooo 
я а јоо ооогооооо» 
e jon о сч олоо ооо ош 
пене ние NN 
| 
1 


сз 
Ф оо ооо 
5 nica тоне бл те en 
mam em д 
i 
o тэ, | со со 
e lio mme 
DEED 
” A ane 
N NN са 
а | Ч 
Ф 
f oj соо о 


4 ге (0 D O 


1 


о 
о 
e 
+ 


о 


с» 


er сч сч 


з,2%8 


7500 


гч б< су 
со ога ч) 
шоо 


сэ хо сэ 0 


oo 
УЖА 
+0 
1 


09 


- 89 - 


46 оеш сі O) од с mina + n N HD So 


М | суб) @ oui cc o oto сз O (O 0 сэ c сб ni nd 
oj oci уа N NO OO 
ЊУ ^ e^ A an he c c n n n e... m on 
< јело ко ت ما م‎ оо о ف م‎ 0 (0 
оооооооосоооосооосооооосоо 
өсоссдоосооосоооФосоосососооо 
e Їорбоообососооооосооосоосоосооо 
Owunouwownooooooooooooood2 
сч єч M c) d it I0 (0 СУ © ОО QU or (0 о O NONO 
хя Үү үү HH сс осот 
~ | чоо 0 Hin 0 OQ r C (0 f О 0 0 O en rt 
ÀA | inu ccm з со сч тал «о г< гу о O O O ооо 
= 1 чом шог о о со оо део Фо оо Ф 
p" аи ACA UR NU, a a Та E DE E Эн А та 
< мю ш AD үз їл (D صا‎ 10 10 иэ мз 10 иэ мэ NI 
соооосооосоооооооооооо 
оооооосооосообоосооосоооо 
alooóooooóóooooooooooo 
u l^ooooóoooóoooonooooooo 
ха Єз: о zb чэ O tS O OY © el CELO о F O ооо 
"EE 
~ |=0о оюшосолшт 
У | оте тооооосоо 
5 |гоо оосу PS с OS 
а, ^ ^n ^a^ ве ^^ 
«T хо ч из иу: Шу из з 
ооо ooooooo 
о|ооо 200000009 
» looo 2920005090 
ооо 009292000 
e њо оо бз 
nad 
~ | мога Nox 
5 | вочо® con 
3 |өсоөоллй аза 
& ^ ^ ne^ ... 
<] ES min 
осососоосооооосоооосооооооо 
n|ooooooooooooooooooooo 
» loononoooooóooooooooooo 
eA сч сч сэ сэ дап ос томно O 10 O «0 O 
vA Hoe e LOL 0 0 dr 
~ | о св сч соо со ол со O t Сї io mono і iu» 
а |хо г. стал о © є- Ne чч хэ лю д > со ою O O ооо 
© فی و ا‎ o م‎ O ا‎ 00000010 D 10 10 (D (О (D 
й Ша а ёсе; той A ^ 6 A A & ааа аа 
4 |oossr srta چ ج جج چم‎ чох 
oooooo 
a 2500000 
e "»ooouo 
ono 
"не 
па «зо о азо о 
шо олоо оо Ф 
ES 


2000000000 
i» O CO O 2 a e 


820 


" 
£a 
я 


е 


t5 


са 
m 


o 

M 
o 
Ф 


n 


“= 


un 


о 
= 


" 
п 


е 


оо 


“ 


m 


o 
2 


а 


х 


32? 
3 


н 
" 


e 
e 


" 


"m 
on 


~ 


ю 


e 


e 


" 


i 


" 
" 
" 
m 
м 


ю о 


ю 


eo 
10 C+ o 
mur 
naonn 
о 


Er 
ч 
" 


o 


шо 


Nom 


мг OO 
о м ост 
e oc 0 co 


о шю 


“ 


u 


ооооо 


CESE EOE 
edo d vd 


и 


6,766 


› © 


10 


" 


n 


aron 
00000 
€ ооз 0 0 Ф 


о «о iD iD D 


" 
" 
v 
“ 
n" 
" 
" 
" 


coooooooo 
А ад math n ^ 
охооо м 0 
HH 
гап OO 
с-оососооо 
p. ta D» (~ 00 со 
nan 
00000 
а 
Боје 
= 
х 
оооооо 
»соосооо 
му tO t 00 Ф ті 
ор HO 0 QU of 
iO су p C. CÓ 0 со 
10 00 0000 
эээ 
(D 4D (0 о ‹о مه‎ 0 
ocooooo 
non GOO Ф 
d ww от оо 
m 
сз On соо du 
божо ао: 1 iO и 
D 0 an i0 i0 م‎ 10 
mr 
ч) о о О о Фо 
0000000 
ne nn na on 
у ому о о © © 
со xf on O ш 
3 6/00 6 09 
0 су сл су Оо» Ф 
1027000 
ай аф ад оо 0 
осоооо 
муо о о о 
ч «зо апо 


хюсюносочч 
ODO Gm сэз O сэ чо O J O оо O 
эз MAN SN IO из ша итэ из O O O 


DD O P DP СМ 6 P ~ ~ ~ 


9 
7,355 


и 
" 
$" 
т 
ч 
n 
м 
т 
" 
n 
" 
2 
" 


n 
n 


оосоооооо 


o 
о о чуо шо Ооо 
e Л 


ч мз zb со со О гу OY? ( iO э 10 ш маат 
I0 ID аю Єч оз с? бо O уо со оз O отто Ф 
DADO eb CONCI op таи ју та 
РРА РИЧИ a ^ n a o 
ао ао fo сз ку POP CS PS DS D P гү 
осососоососоосюсоосооосооосоо 
таз жо т сз озо о сз с) о чече жю 


como D cde O y dudo 

сом р coe (oH т оо [~ 

f- C 09 03 QD үз c4 сч о со 4) сс о сз (0 СЭ e 
о ~ 


озогу по ~ C DS ~ ~ 


“ 
и 
n 
п 
" 
n 
" 
" 
и 
и 
ч 


ч 
и 


"n 


ooooooooooo 
E AI ја 

Q9 DO N QU O10 O QT 
AANAND 

тош зз о оо суо кєз суя о о Ф 
HOM FID CG O NID съ сч tt cour 2 
хо г о т о OHH HNN MON Сч N сч 
^o^ ө n^ ^ nim а ал а а о о а ке 
мэд чу вз б f» сз (ч 0 D P СУ P p D М [9 
S EXE = БЕБЕ E BE 


1500000000000 


o 
Moe eo^ n^ ^ ^9 а o h^ nn se 
MO cb D) «ро го) OVW о i0 O INO 


apo €) чэ сэ со UO 1- O 10 O ПУ 

Q m «4+ C ж ID 10 O 00 по (O O о ©: оз 

с о) оз сл ср с сз Ср OO © ©з OOOO 
~ 


о о о x) CS (9 со (~ 0 ES СУ 


м 
n" 

n 
n 
" 
ч 
" 
м 
r 


1 
т 


000020000000 
Гэгээ 
DON F O (0 m (o eb оч сэлт 
SS 


му O CY оз саг о WD 
ус о о Ф о о 
о досоооо 
"^a^ nn. 


с су су озо 00 со 


t 


ОЕБЕЕЗЕЕЯЕ СЕЕЕРЕ 

a 

x 

Дд.солосоолооуосооососоосоооо 

ка non mom mom him пл A m ЭЭ 

4 Үү ed doen NC сз о) со до шоо (S 00 Ф О NA 
"че 


хз O сч со сч о 
oon Od oce 
OQ со aan о ст 


c. Berg C P D 


= о O10 соо ct ODM 
о, > —очово со 00 
&oooooouoo 
.. 
~ ~ 


7,775 
7 


(о по ге Om ES ~ P P 


n о 
ОБЕ БЕБЕ эээ 
< ч 
х * 
фоооосоосооосооосооосоосоо 
Эа а ае а а а юв т m^ enhn n = 
d£ ш с сч б сг ср зз «з мэ бо © оз с) о ОЧ? 
voe ex ба CI O4 (0 КО CE мә (D г 00 0 elc rt 


Ton Ow d 

€) суму et OO рх c m m o Т 

аз QUON OO uox mI Qo tc 
~ 


ON 


үз O CI с) 


E 

ПОББЕЕВБЕБЕБЕВБЕВЕБЕЋЕЕЕ 

wi 

х 

Дсооосоооооосоооооооо 

i$ с 00 € (CO мг foco (2 мз «э O O O ооа 
днс N N со со токо 

(€. pm peo € t Оз n ed CS O сэл ОО 

«з ал л л нэ сэ 

єл 00 эчу ко '0 2 


n «ож а nn «on sn nn. 


NN 


v 
" 
" 
- 
т 
" 


0000000020000 
а о Domo © оо 
И ко сој 


1854: зерә) 08: D G'| sp ¡DUOISUSIMIPD огроц зомзцха энш je ое 2)40)5002 42152.83 А ациосооэ { чо! пи) мотор > 
-^id ао ором оршу) орсиын (0:24 ошај 5 Фр Чошовшдшнро ошєц ар силу ©шоэ [oucisueunpo voiso;d эр окисо0-1205:4 


го 550 edo 500 


aio). 


1011 fe чим им) Ру 


опџциоо oppis3 


xou (91) dy 


{84 "si3398)"06 © Ol эр |200150әштро паром JO:J2/x29331 ve DLUD¿SUO) 0015920 А aJUD¿SUO) (опции) upto 
-элро ap ороз (0:804)) оргуш (0 00: Cues} чохогөсэцарэ Ode әр UDI} OO) (2:0:508шир2 0016910 ёр 0:4ш02-22 5:3 


б, "лумозаато“ байа 


„а è 9 T 2 „а 0008 соса СОС» 0002. ссо 053 069 207 002 ос. ce оз 


у 


AE 


pr зарын! 


яасын 
E 


ME 


i (onunso? op 


2153) 


n 
« 


M 


f 
| 


(91) 40 “зүхогэм208014-110153н4-30 0183 


кете 


тт 


t 
[s 
НЫ 
F 


4854 'а9јен) 006 2 OO! эр (04015 


миро строен 10:08:38 арми, 19 68 243045002 4916229 Á aLUDISUCI (90:80:41) 


:sejd ap сашо2-СсС Dhu 


o 
о 
E 


Ф 


e. 
Go 
Фе 


0 


n" 


Gay Q QO 
во 
y 


900 


NUO? Ср2155 


A a 


39 отму 


dy "NNOISN3WIO Мона 


(8) 


16914 А ajun; 
2 124015096. 


про 1) 


1882 


suo) (UoIen:tut) vorooap 


од 


мэр Ч 

¡O ~ о 

ооо 
гаг 


~ 


10 
^» 0 0 t) 
ММ ММ МА 


00 
20 M 


ONO may 
оч 


[24 


ро 


3 


24 


A b ta 


«mmo 


AVI 


id Old 


5384 


ау NO! 


VNOISNJ* 


ta 
1 


dg 


) 


(8) 


~ 97 - 


де 
W(R t) = We 5 1,119 $ e hR ——— QU) (69) 
3609 
donde, 
( )2 1 2: 2 
гер ај e e^?n^D Jilan Ten) 
Qt) = QU, tp). — zie» рше ен 10) 


=1 42112 2 А 
asi а (32 (ap) - Јубар vep?)] 


donde ар зол las raíces de la evcación, 


11 (ап Tep} Yalan) - Yilan тер) Jolan) = 0 (71) 
Si ty < 0,01: Usar la Ec.(52) para calcular Alty). 


Si tp 2 0,504527, ocurre flujo continuo, es decir, el valor 


de Q(tp) es máximo y toma un valor constante calculado рог, 


Q(tp) i SKS SEE (72) 


Los valores de Q(tg) son dados por vEH7, para nap entre 1,5 y 


10 y en la Referencia 31 para entre 20 y 10%, Тоз resuitados aa 


Pep 


representan en la Tabla VII y Figs,35 a #6. 


Caso 6 - Presión constante en el límite interior (contacto a~ 


gua-pstróleo) lo mismo que en el límite exterior, Analíticamente, las 


condiciones de límite son: 


а) lgu2l a la Condición a) del Caso 2 


b) Igual а la Condición b) del Caso №. 


M 
Este caso fué vosualto por duellor?! y permite nl cálenlo del Fluido 
total que раза por el limite interior (rg = 1) en un ticmpo dado t, 


usando la Ec.(51), donde sl valor de Q(tg) fue calculado por computado 
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TABLA VII - Valores de Intrusión Adimensienal, Q(tp), como Fun- 
ción de Tiempo Adimensiounal, tp, para Valoves de гор entre 1,5 
у 10%, Caso 5; Acuífero (yacimiento) Radial Limitado (Cinito). 
Presión en el Limite Interior Constante y Rata de Producción en 
el Límite Exterior Igual a Cero (Acuifero Cerrado), Valores de 
гер entre 1,5 y 10, Ес. (70), Referencia 7 y 8, los demás Valo- 
res, Referencia 31, Figs. 35 а 56. 


0,638 


0,250 0,578 0,50 1,02" 0,125 1,791 

0,275 0,715 0,55 1,083 0,150 1,997 

0,300 0,751 0,60 1,140 0,175 — 2,184 

0,325 0,785 0,65 1,185 0,200 2,353 

0,350 0,817 0,70 1,298 0,225 2,507 

0,275 0,848 0,75 1,229 0,250 2,648 

0,500 0,877 0,80 1,338 0,278 2,772 

0,425 0,205 0,05 1,395 0,300 2,888 

0,450 0,932 0,90 1,480 0,325 2,990 

0,475 0,255 0,05 1,984 0,350 2,08% 

0,500 0,983 0,10 1,526 0,375 3,170 

0,550 1,028 0,11 1,605 0,800 — 3,257 

0,600 1,070 0,12 1,679 0,425 3,317 

0,650 1,108 0,13 1,757 0,050 — 3,331 

0,700 1,153 0,1% 1,811 0,475 3,439 

0,759 1,174 0,15 1,870 0,500 3,991 

0,800 1,203 0,16 1,924 0,550 3,581 

0,900 1,253 0,17 1,975 0,500 3,656 

0,100 1,295 0,18 2,022 0,650 3,717 

9,38 0,603 0,11 1,330 0,20 2,106 0,700 3,767 
0,10 0,606 0,12 1,358 0,22 2,178 0,750 3,809 
0,45 0,613 0,13 1,382 0,20 2,281 0,800 3,843 
0,50 0,617 0,14 — 1,902 0,26 2,29" 0,900 3,89" 
0,60 0,621 0,16 2,532 0,26 — 2,950 0,100 3,928 
0,70 0,623 90,17 1,059 ,30 2,589 0,130 3,951 
0,80 0,62% 0,18 41,953 0,38 2,299 0,120 3,967 
0,20 — 1,968 0,38 2,891 0,480 3,985 

0,25 1,987 0,42 2,525 0,160 3,993 

0,90 1,135 0,46 — 2,551 0,180 3,997 

0,80 — 1,499 0,50 2,570 0,200 3,999 

0,50 1,500 90,60 2,599 0,220 3,999 

0,70 | 2,613 0,280 14,000 
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TABLA VII (Continuación). 


+ Е о о фо бо о (о له‎ о № 


ی ی ی ی 


à ruo 


nu 


eD = 4,0 eD 
tp tp) Ер 
2,00 2,442 2,5 
2,20 2,538 3,0 
2,40 2,788 3,5 
2,60 2,893 5,0 4,5 4.193 
2,80 3,038 5,5 5,0 à 839 
2,00 3,170 5,0 5,5 — 8,792 
3,25 3,33% 5,5 6,0 5,074 
3,50 3,493 6,0 8,5 5,355 
3,75 3,645 6,5 7,0 5,605 
4,00 3,792 7,0 T5 5,854 
4,25 3,932 7,5 8,0 6,09% 
1,50 5,088 8,2 8,5 5,325 
4,75 4,198 8,5 9,0 5,517 
5,00 1,323 9,0 9,5 5,756 
5,50 4,560 3,5 i0 6,955 
5,00 5,779 10 11 7,3506 
6,50 4,982 14 12 7,796 
7,00 5,159 Не) 13 8,035 
7,50 5,343 13 1» 2,339 
00 5,50% 14 1 8,620 
8,50 5,653 15 15 8,870 
5,00 5,790 16 18 2,338 
‚50 5,917 19 20 9,37 
5,035 20 22 10,07 
6,246 22 24 10,35 
6,425 2% 26 10,59 
6,580 28 28 10,80 
6,712 28 30 10,98 
8,825 30 34 11,76 
6,522 33 38 11,98 
7,00% i8 42 11,61 
7,076 “2 45 11,71 
7,159 Hh 50 11,79 
7,272 50 60 11,91 
7,332 60 70 131,96 
7,877 10 20 11,98 
7,434 80 90 11,59 
7,4604 20 100 12,00 
7,4081 109 120 12,0 


окр) 


6,851 
1,127 
7,389 
7,902 
8,397 
8,876 
9,351 
9,791 
10,23 
10,65 
11,06 
11,46 
11,95 
12,58 
13,27 
13,92 
18,53 
15,11 
16,32 
18,41 
13,83 
19,24 
20,51 
21,45 
22,13 
22,53 
23,00 
23,47 
23,71 
23,85 
23,92 
23,96 
24,00 


uo 
45 
50 
55 
60 
70 
80 
80 
100 
120 
150 
160 
180 
200 
240 
280 
320 
362 
400 
500 


101 


24,94 
25,53 
28,11 
26,67 
28,02 
29,29 
30,59 
31,61 
32,67 
33,65 
34,50 
35,48 
38,51 
40,83 
42,75 
44521 
45,36 
46,95 
47,9% 
48,5% 
48,91 
49,14 
49,28 
49,35 


19200 
13000 
16000 
20000 
23000 
30000 


50000 
50000 
100000 
110000 


2003 
2460 
2520 
3230 
3560 
3940 
370 
4629 
4770 
4920 
4970 
4980 


13000 
16000 
20000 
24000 
30000 
50000 
50000 
60000 
80000 
100000 
130005 
160000 
200004 
240000 
200000 
500000 
500000 


9680 

11700 
13300 
15200 
16500 
17800 
18500 
19200 
19700 
19000 
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Pep = 500 Гей 1000 
to QCtp) tp QCtp) 
100000 17500 з х10" 5890 
130000 22100 вх" 7649 
160000 26500 5 x" 9350 
200000 32100 1 x10 17699 
240000 37300 2 x" 33200 
300000 94700 $ x 48000 
400000 55600 ох" 62500 
500000 65000 5 x" 76400 
600000 73100 6 x" 83300 
800000 86200 7 x" 103000 
1000000 96000 8 x" 116000 
1300000 106000 9 x", 128000 
1600000 113900 1 x109 140000 
2000000 128600 1,4х " — 183000 
2400090 121000 2 x" 239000 
3000000 123000 2,"x " 271000 
4000000 125000 3 x" 311000 
5000000 125000 4 x" — 263000 
5 x" 601000 
7 x" 048000 
8,4x ", 5870900 
1 x10 480000 
1,8х " 495000 
2 x"  ^99000 
2237 500000 


Ба 


tp 


D 


X оо‏ ت ص حم CIN NIB кэ‏ چ I) эхээ мл‏ دم ш‏ چ ہی ہہ نو ر 


x1 


[ 
о? 
" 
" 


2000 


обр) 
52700 
77000 
91100 
105000 
119000 
133000 
145000 
106000 
225000 
276000 
2256000 
3570029 
510002 
610090 
804000 
10530909 
1230002 
1520090 
1535000 
1650000 
178006000 
1530000 
1950000 
1590000 
2009000 
2000000 
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TABLA VIL (Continuación), 


ср = 4000 


tp Q(tp) 


1 
| 


х106 н,06х10° 
" 


1,6x109  2,26x10? 3 $,77x108 
2: 7X4 2,78x " вх" 5,31x 7,83x " 
2,4х " 3,30% " 5 x" 6,76x " 8,86x ", 
з x" 4,06%" 6 x" 7,76х " 9,93x10 
uox" 5,31x " 8 x", 1,02x " 1,10x " 
5 x"  6,5üx " 1 x107 1,25x " 1,51x " 
6 x" — 7,Tux "| 1,2x "| 1,49х " 2,11x " 
8 x", 1,01x109 их " 1,72x109 2,70x " 
1 x10 1,28х " 1,6х" 1,95% " 3,00% " 
1,2x " 1,48% " 1,8x " — 2,17x " 3,37x " 
1,4x ®" 1,57% " 2 x" 2,%0х " 3,56х " 
1,6x " 1,87x " @ lux Qd ac "03x " 
1,3х " 2,07 & " з x" 3,515" 4.u0x " 
2 x" o 2,27% " u x" — 4,60% " Ob 
2,9x "  2,58х " 5 x" 5,89x " 5,30% " 
3 x" " тох" 5,824 " 
ox" " 8 х" 617% " 
5 x" 5 б с 6167x " 
7 к" " 1 x10? 7,01x " 
8 x" " 1,2x " р 7,50х " 
9 x" " сүйх " Bux "o 7,824 " 
1 x108 " 1,74 "  1,66x " 9 x " - 8,204 " 
1,2x " " 2 x" 1, 3х " 1 xio? 9,06% " 
qux " Ч A "o пи 1,3x " — 1,12x109 
1,7x " " 2 x" — 2,59x " 1,64 "o 1,32x " 
2 x" " ох"  3,00x " 2 x" 1,55х " 
Ziar t n SN 3,47x " 2,%х " 1,75х " 
2,6% " " 6 x" 3,79% " з x" — 2,00x " 
3 x" " тох" чтуу " goa" 2,33х " 
Тр У !! 8 х и 5 wx" у 63 2 
5 x" u 2» АТ " 8 ox" 22x " 
1,%x ” 8 x" 2,92х " 

2 х " 1 x10? 3,02x " 

3- ow У ),ux " 3,10x " 

Оо 3.1254 " 

3 x" 3,12х:" 


TABLA 


Pop = 1000000 


tp 
s x10)  u,s1xio? 1 1011 

3,5x " 3,59x " вх" 5,37х * “Зи " 
цих" 4,07x " Two 6,22x " 1556" 
5 x" 5,03% " 8 x" 7,07к " 2х" 
6 x" 5,98x " 1 x1010 8,75x " 2,4x " 
7 x" 6,93х " 1,9х " 1,21x10? з x" 
LEE NEU 7,87х " 27) x: 1,70x " uox" 
B.ux " 8,25х " з x" 2,49% " 5 x" 
а x" 8,00x " вх" 3,27x " 6 x" 
1 x109  9,72x " ux 3,58х " 7 x" 5,11x " 
Lx " 1,34x100 7x " 3,80x " 8 х" 5,80х " 
Lop” 1,88x " 5 x" и,03х " C 4 x1012 7,13x " 
з x" 2,77% " 5,ux " 4,32% " 1,2x " вучја " 
ох" 3,68х " 8 x" 4,76» " 1,Ux " 9,67x " 
5х" 5 49x " 1 x" 5,90х " 1,5х " 1,03x10! 
$ x" 5,33х " ПАРИ 5,70x " 
7 x" 8,15x " в x" 6,17x " 
8 x" $,95x " ix " 6,42x " 
9 x" 7,75% " 9 x" 8.85х " 
1 x1010 3,52x 1 10l1 7,81x " 
1,3x " 1,08х 1,1х " 8,15х " 
1,6х " 1,29x " 1,3х " 0,38х " 
2 x" 1,56x " 1,6x " 1,11x1010 
2, х" 1,80x " 2 x" 1,32x " 
з x" 2,1йх " 2,9x " 1,51x " 
uox 2,62x " з x" 1,79x " 
$ x" 3,03х " ох" 2,08х " 
вх" 3,36x " 5х" 2,33x " 
7 x" 3,69х " 7 x", 2,8x" 
8 x" 3,87% " 1 х1012 2,92x " 
9 x" н,06х " 1,3x " 3,084 " 
i x10 4,22х " l.óx C 3,09x " 
1,2% " "85x " 2 x" 3,11x " 
1,6x " 74x " 
2 x" n, BBX " 
2,4x " Ny Ox " 
3 х" 4,98x " 
Зак! 4,29x " 


5,00x " 
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тез y utilizando métodos nunmíricos para resolver la esuación de difusi- 


vidad. Los valores de Q(tp) para тер entre 1,2 y 100 se presentan en 


las Figs.47, 88 y 19. 


2eniferos Radiales Мо-Ното 


En forma similar que para Acuíferos Lineales, Muellen?A pra- 


soluciones para los Casos 3 y 5 de acuíferos radiales, donde í- 


gualmente varian linealmente las propiedades espesor, h, permeabilidad! 
: t 
viscosidad, k/u y (porosidad)x(compresibilidad), фи, entre el limito in 


terior y el límite exterior. Define igualmente el Factor B como la pro 


piedad variante en el limite exterior dívidada por la misma propiedad 


variante en el lími.e interior. Сото puede ob3arvarse de las solucio- 


nos da Jos Casos 3 y S, habrie пп número ilimitada de soluciones, ya que 
pura cada valor d^ Top existe una solución especifica. Mueller consi- 
deró los Casos 3 y 5 sálo para гор=10,0. Pate valor ez tipico dei ta~ 
maño de muchos acuífaros en el mundo. ba Fig.50 presenta los resuita- 
dos, ар (р) como función de tp, correspondiente al Caso 3 


del =spesor, según el valor de В entre 1/3 y 30. Та 


los resvltados,Q(tp) como función de tp, correspondisuz 


ra variación del espesor, de acuerdo al valor de р entre 


el a 


Las F3ys.52 y 53 presentan soluciones similares, 
5, donde varía el factor premeabilidad/viscosidad y eli tas 
dad)x(compresibiiidad) respectivamente pora valores de B асте об y 


1000 y 0,1 y 10. 


Lo que so discutió respecto а los factores k/u y Фе para ei са 
so de acuíferos lineales, se cumple también para acuíferos radiales. чо 
{а similitud entre las Fign.1" y 52, би comportamiento es алаіоцо 


tes 


por las razones discutidas al tratar la Fig.1i^. 


t 
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ОЗЮМ ADIMENSIONAL, Q (tp) 


INTE 


Fig. 47-Intrusión acimensiono!, 006), como función de tiempo odimenzionol,to. 
; Sistemo rodio! limitado (finito), Presión, tonto en el límite interior (coniocio 
Ї сбив petróleo), coma en el limite exterior permanecen constontes. Rodin odi- 

racacionol de 1,2; 2,0 y 3,0. Tiempo odimensionol ce 0,0 a 1,6,/Бо!огел- 


o6 08 O9 м 1,1 
TIEMPO ADIMENSIONAL, lo 


iis :: "edi : ШЕЕ į: 

48-Intrusión аф по!, О(10), como función de tiempo odimensionol, to 
“Sistema гайг! limitado (finito). Presión, tonto en cl límite interior (contoc- 
о agua—petrolzo), como en el límite exterior permanecen constantes, Ro- 
io adimensional de 1,2 а 6,0. Tiempo odimensionol, 0,0 a 18.(Rete- 


rencio 211 


7 10 11 12 13 14 16 18 
! TIEMPO ADHAENGIONAL, to 


lnirusicn duna sj omno mensionol; to, E 


(finito) Presión, tonto en el límite interior оша! 


ig. 49 
Sistema redha! limitado Ч 
о), como en el linite слепог permanecen constantes . Rodio | 

| 


agua „реше 
adimensione! de 4,0 a 100. Tiempo adinensianal de 0,0 о 180.1 {Rofe ~ 
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SUDERPOSICION 


En las secciones anteriores se discutieron las diferentes so- 
luciones más comunes para sistemas lincales y radiales, 4о046 la rata 


ite inteeior, perma- 


de intrusión (o producción) о la presión en el lix 
necen constantes соп tiempo, es decir, los términos e (o q) ур (o Ар= 
рї-р) son constantes con tiempo, о sea los casos de rata terminal o pre 
sión terminal constante. Sin embargo, ésto en muy ravas ocasiones oc 


tre en la práctica, pues es bien conocido que tanto e (о q) сомо p (о - 


Ap = pi-p), cambian a medida que la intrusión o producción də un уасі. 


miento aveaza. Con el fin de aplicar las ecuaciones a estas situacio 


nes se ha hecho uso del principio matemático denominado Teorema de Su 


perposición o simplemente Peincipio de Superposición. Matemáticamente 


el principio o teorema de superposición estublece que la conbinaci 


ticulares do una ecuaeión diferencial homo 


neal de soluciones ра 


бп de la ecuación diferencial. 


lineal ss también una solue 


imponen condiciones de 


pio es muy util en sistemas donde 


olución general es la suma de las scluciones particulare 


torno). La 


obtenidas tratando una de las condiciones de limite a determinado tiem 
po. 


La obtención де la hi 


Ei principio de 


tocia de presión, dada una historia determinada de rata de intrusión ( 
prodicción), o viceversa, puede obtenerse de la rata de intrusión (o pr 
ducción), dada la historia de presión. Además sj se coucce la histori 
de presión, puede calcularse la intrusión total de agua al yacimiento | 


la producción total de petróleo del yacimiento) durante un periodo dad 


sta filtima aplicación del principio de superposición es más comün cua 


do se trata de problemas de intrusión de agua a un yacimiento, El pri 


mer caso, о sea cálculo de presiones, dada 14 historia de rata de pro- 


ducción, es más común en problemas relacionados con pruebas de pozosi8, 
о en problemas de acuíferos cuando se conoce la rata de intrusión. Por 


ejemplo, en yacimientos con un empuje hidráulico activo que producen por 


encima del punto de burbujeo, rata de producción de petráleo puede 
asumirse aproximadamente igual а la vata de intrusión de agua. La ra- 
ta de intrusión tanbién puede obtenerse їе datos de presión-producción 


de petróseo y haciendo un balance volumétrico de los fluídos en el ya- 


cimiento. 


Cuando se conoce tanto la historia de pre storia 


гоп como la 


de Ja peta do inteusión puedo determinarse la denominada función de In 


3,3% 


fud E PR : :a - 3 
fluencia definida como el cambio en presión del yacimieaco рог ра 


za do iutruzión unitaria. 


La Pig.5" ayuda a doducirv fórmulas del principio de superposi 


с 


Se conoce la histroia de la rata de intru - 


ción para el primer caso. 


én (o producción), e, уе e .. C, que comienzan 3 actuar respecti- 
І > 0,02» ба» n 1 


vamente а los tiempos, 0, t4, t5 e Бао En la práctica, la hosturia 


de la vata de intrusión (o producción) se presenta en forma continua i 


nea contínua), sin embargo,con el Kin de poder aplicar las 


a anteriormente (rata terminal 


попка ве), ze aproxima en Forma es 


la presión ini - 


calonada, tal como indica la Fig.5". Adomás 


cial (or.zinal)del yacimiento. 


Supóngase que se desea determinar саїла сока) de pre 


por сот 150) à un tiempo, Tp, comprendido entre x; y 9. 


isujente Ја pr 


Ya que en este intervalo, la rata de anteusión (o ргодисслби) е5 cons- 


tante, сү, 8e aplicarán directamente las Ecs.(28a), (284) y (290), para 


diccid 


Ms] 


intrusi 


о 
= 
5 


Ga 


iz 
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га 
un рогоүо уз 


Historia и Presi 
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АР 
x. «у 

E 
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el caso de acuiferos lineales, y las Hes.(37), (59) у (63) para el са 


cuíferos radiales: 


с 


зо de 


Рага 0 < t < tys 


Acuíferos Lineales, 


287,58 ey ub 


Ар = (рї-р) = — > ву (091 (73) 
kbh 
Acuíferos Radiales, 
141,2 е и 
Ар = (рі-р) = — A ap [ (1) p) (37) (59) (63: 


kh 


Al tienpo, Хү, la rata de intrusión (2 producción) aumenta por цаа саг 
tidad (e2784). Este aumento en la rata de intrusión lógicamente causa 
24 una caída adicional de prosión y la caída total da prosión Су рог 
consiguiente la presión) se caleulará, para un tiempo t, comprendido en 
tve t9 y t1, por la superposición de dog efectos: ey que actua por el 
tienpo (t-0), más (е2-с1) que actua рот el tiempo (44:1). En bass a Заз 


ecueciones anteriores, puede escribirse la caída total de presión: 
Para tq < t < ta! 


Acuíferos Lineales, 


887,58 ey ub 887,58(85 - 0478 
Ар = (pi-p) = === SI + ——7— A AE 
kbh xbh 


Acuíferos Radiales, 


141,2 ey u 
Ар = (рі-р) = Lii Аль СЕЗ] 
kh 


-e4)u 
- X---ap[CGeotiy] ( 
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capo te, la rata de intrusión (о producción) cambia en una can 


Si al 
tidad (eg-e73), éste nuevo efecto causará otra caída extra de presión, 
y la caída total de presión (y por consiguiente la prasióv) para un tien 
po, ts comprendido entra ty y t2» será la suma de tres efectos; еј que 
actua por un бістро (t-0), más (ez-e4) que actua por el tienpo (t-t4), 
más (ед-2,) que actua por el tiempo (t-t2). En base a las Ecs. (73) y 


(37) puede escribirse; 
Para tg < t £ tg: 


Acuifovos Lineaios, 


087,58 су pb 887,58(02-01)н9 

1 TOES 2323 { 1 

Ар * (pi-p? = A Ру р) Sy AA О + 
kbh kbh 


887,5f(e3-e9)ub 


Е — Fe (Go 1235] (75) 
kbh 
Acuíferos Ra les, 2 
141,2 ey y 151,2(e 
Ag а (pi hp a ТҮЛСЭГЭ) oM + 
kh | kh 
181,2(e35-65) . 
‚дм 2 -ар((1-5293| (77) 


kh 
Si al tiempo tg, la rata de intrusión (o pecducción) cambia en una can 
tidad (94-04), рата calcular la caida total de presión a un tiempo, t, 
mayor de t3, basta agregar esto efecto a los auteriares, y escribir e- 
cuaciones similares a las Eos.(76) y (77), agregando otro término de- 


bido а (54-723). Asi podría continuarse, para una secuencia de en ra- 


tas de intrusión tal como iudica la Fig.55, y ја ватда total de prec- 


~ 130 - 


а un tiempo ty puede escribirse, 


c 
o 
с 
5 


Para tn > ta-1' 


Acuíferos Lineales, 


887,59 и L 
2 (ta-t4)p] + 


„381,58 и С * (e2721) xd 


Ар = (рї-р) = 
kbh 


+ (03-87) тиг + ... + Centen) тиг (7 


» 887,58 Ч 
Ар 3 (pi-p) = cm | [Cedo] + (76 
kbh 1° 
237,58 pk а 
Ws (pip = е 1) Fpi tantj-10] 79 
kbh j 
donde eg = 0 y to = о. 


La Ec.(78) puede arreglarse en forma diferente, 


Ap = (рїзр) — 


987,58ub | 
СА Ру 05 Уу) + ту 0К-5436)| * € 
kbh - 


à 
| 5 . Чи > 2020 
a Ся) 0 و‎ =: Ру Стат 
+ An Fg (Си EOM ( 
887,52 uL P | 
Ap = (р-р) 5 0707 АЕ; ry [C n7ta-j?0] D гү [Cn 5241-15) 
kbh је 
(81) 


donde ty + 0. 
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Acuíferos Radiales, 


101,2 y 
Ap = (pi-p) = irr е1 ap [Cra p] * (05764) ap[(tn-t425) + (63-01) 
ар Султан) + +5 lenten? ap[Ctp-to-1)0) (02) 


151,2 eq u 


др = (pi-p) = ——— ——— apí(tadp) + (83) 
kh У 
141,29 P : 
Ap = (pi-p) = ме 77 I (ел-е4.1) ap[tta- 12) (аа) 
А j 1 
kh 131 
donde во = 0 y to = 0. 
La Ec.(82) puede arreglarse en forma diferente, 
141,2 n | 
ap = (рі-р) = ene M Ap (Ен! ~ PEG tul] У sa во Сато - 
k 


= КОА Fore enol ^Ap[tg-ta-2)p] 7 по e-2 | + 


+ од ApÍCta-tn- 125] (95) 
| 141,2 џи P 
Ар = (рї-р) د‎ = LO ср+1-] пр (катка) > Арсы бъл) Ч 
kh j=1 
donde ty = 0 


Debe tenerse presente en las ecuaciones anteriores (:2.(73) en 


te) que el término p es la presión al tiempo t o tp. Ya que pi 


noce, puede calcularse p, una vez calculado Ap. 


- 13% - 

El principio de superposición con el fin de obtener la intru 
3161 {o producción) total, dada la historia de presión, es muy similar 
al caso anterior. | La Fig.55 ilustra la historia de presión de un yà- 
cimiento (о contacto agua-petróleo). De nuevo, con el fin de poder a 
plicar las ecuaciones vistas anteriormente para el caso de presión 
terminal cosntante, la historia de presión se aproxima en forma escalo 
nada, y se determinan las caidas de presión (presiones diferenciales о. 
decrementos de presión) Ару, ^p2» Аруу see Ард» (Sobre la forma de ob- 


tener astas presiones diferenciales, véase p. 145). 


Sa desea calcular la intrusión total a un tiempo, t, compron- 


dido entre 0 y tz. En este caso, el Арі #5 constante durante este tiom 


po. Por lo tanto pueden aplicarse las Ecs.l Ta), (304) y (31а) para 


et caso de acuíferos lineales y las Ecs.(91) y (69) para acuíferos3 ra- 


diales. 
Para 0 < t £ tg! 
Acuiferos Lineales, : 
¿bhLe 
We = m АР Ру [(Е)р] i (87) 
5,615 
Acuíferos Radiales, 
0 
We = 1,1199chf --арү (00001 (51), (63) 
3609 
Hasta ol tiempo t4, la peosión se ha mantenido constante сл 
р торхї-Арү» Si a óste tiempo t4, la presión disminuye en un Apo fp = 


ра“бр1- 492)» habrá un aumento extra 6n la intrusión para tiempos ров 


so desea obtener la intrusión 


torioras. Si el tiempo, t, а la 


total está comprendido entre t, y t9, tal intrusión se calculará por 


- 133 - 


la suoeeposición de los efectos: Ару que actua por el tiempo (7-0), паз 


Ард que actua con el tiempo (t-t,). En esta forma puedo escribirse, 
Para 0 < t < 54: 


Acuíferos 11008165, 


фъһ1с $bhlc 
We = ——— ap, fx[(005)] + ----т ôP? rkf(t-t,2p) (88) 
5,615 5,615 


Acuífnros Radiales, 


0 0 
Ho = 1,119 échRT ————- ^p, q[tt55] + 1,1198cRÍ ————Ар2 0106-5195) (89) 
3509 3609 


Similapmonte, 51 se desea calcular la intrusión total a uu tiempo, t, 
comprendido enire 09 y 23» Y teniendo en cuenta que el tiempo tg, ocu- 
veo una caída extra de presión Ард (la presión ahora será р = pi-^pi- 
-Apo-hpa), el :focto total sorá la suma de tres efectos: Ap, que actua 
por ua tiempo (1-01), Ар; que actua por un tiempo (t-t,) y Ард qreactua 
рох un tiempo (6-02). Рог tunto, puede everibirse, 
Роза t4 < ts to: 
- 


Acuíferos Lineales, 


$bhLc форе 
Ма = Ару Ру [00р + Аро reÍCt-t42p] + 
5 5,615 


Acuiferos Radiales, 


- 1364 - 


шоо 9002) + 1,11096884 Pc aftct-t42,) * 
| | 


We = 1,1196chRÍ 
(91 


+ 1,4198ehR] -— 9 арз ЗООГ 
350? 


Si al tiempo їз, ја presión disminuye айп más, por una canti- 
dad Арц, para calcular la intensión total à un tiempo, t, mayor que 
t3, basta agregar 21 efecto del Ард а las ecuaciones anteviores en for 
ат а Арц caídas de pre- 


ma similar. Así podria continuarse hasta 11 
sión, y la intrusión total a un tiempo, tns será, 


Pava ty < tn-1! 


Acuíferos Lineales, 


hp, РК! Сеј АРа шиг 


тэнчооал i (3 
ува n , 
Ма = о Е: Ар; ка 4120] (3 
„015 21 


" 
o 


donde to 


Acuiferos Radiales, 
| у n'p 2 м/р 
о ДАр 5144) ) + Ард Q (ta t.) } + 


, 
ма = 1,1199chn2 
: 3609 | 


+ o... АРр QfCta- t4. 35] (94 


0 а 
че = уууллвола олон z anjal ван 130) 
36 jsi | 
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ANALISIS DE YACIMIENTOS 


ON EMPUJE HIDRAULICO 


El mátodo de análisis de un acimiento (о campo) que presen- 
y q 


ta intvusión de agua (empuje hidráulico) depende de los datos disponi- 


bles, de lo que se desce determinar y de las suposiciones sobre солрог- 


tamiento futuro dol yacimiento (predicciones). En la presente sección, 


se discutiván diferentes mótodos de utilizar la historia del yacimiento, 


con el fia de det :rminar el tipo de intrusión de agua; ademas, se discu 


te el análisis de la historia con cl propósito de llevarla a axpres iones 


analíticas con el fin principal dc hacer predicciones. Posteriormente, 


-unos de los métodos existentes de predicción de avuer= 


3 


зе añalizaván 


do а lo que se desce determinar: presión, rata de intrusión, eta. 


ón de Agua 


Le Expre ас la Int 


Una de las formas más sencillas para determinar el tórmino de 


inteusión We, es por balance de materiales. Lógicamente, ¢ 


mê todo 


asume conocidos todos los términos de la ecuación de balances de máte- 


A 
rialss, excepto Wo, es decir, asuma conocido el petróleo originalmente 
en el yacimiento, М, el factor m {сп caso que exista сара de gas  ini- 
cial), datos de producción y datos РУТ, 


Cuando se conoce con bastante exactitud los límites del yaci- 
miento y se dispone de un buen mapa isopaco, el valor de N puede deter 


actitud por ol método volumétvico, 


sante e: 


ми обе con 


7758 АМ9(1-8,) 


р нм (86) 


Boi 


= 130 - 


douge N se expresa en barriles normales, DN; А en acres; h en pics y 


Bg; еп BY/BN. La saturación de agua connata, S, y la porosidad, 9, 


м 
pueden determinarse por registros o de análisis de núcleos. Boi pue- 


de conocerse de análisis РУТ o por correlaciones. 


Otro mátodo9$ para determinar N, consiste en calcular собе 


valor de la ecuación de balance до materiales a diferentes tiempos, con 


siderando Wez=0 (m, en caso de existir, дође ser conocido), El valor 
de N calculado aumentará con tiempo y la extrapolación gráfica a t=0, 


será el valor correcto de М. 


De lo anterior puede aZirmanse qua el método do aná 
la historia ^el yacimiento con el fin de exprosar la intvusión de agua 


por una fórmula determinada, depanderá de si se conoce o no con sufi- 


ciente exactitud el valov del petróleo originaimente en el vaci 


Por “auto, la зе dividirá en dos partes principales: М cono- 


cido y N desconocido. 


Original 2 


53-09 conoce este valor, pueden calcularse una seris © valores 


de Uo, por balance de nateriales (m debe ser conocido en caso 48 yaci - 


de Ме сото 


mientos de petróleo con capa de y del análi 


-petróloo, pueda ob5teas- 


función de tienpo y presión en el contacto арч 


se una expr:;sión рага voprezontar la intrusión de agua. Una vez obteni ° 


da tal expresión, puede роодосїгсо una intrusión futuca, de acuerdo al 


método usado, como so verá posteriormente. 


La siguientes ecuaciones de balance de materiales! за usan ра 


ra caleular la intrusión de адил, Ма: 


= 1397 = 


mai 
mi 
a) Yacimiento Че petróleo no-satuvados (por encina del punto 
de Burbujeo), 
We = MBoi - (N-Np) + Wp - Wi (97) 
ya que por encima del punto de burnujeo, el 32 dupende de la coupresi- 


ribirse, 


bilidad efectiva Ces puede ase 


Во = ој + Воз celpi-p) (98) 
donde, 
lg = (39) 
indo esta expresión en la Es.(98), puedo esepibivso, 
Ма = HyDo ~ 1301 eo(pi-p) + Up - Yi {100 
b) Yaciniontos de petvólco satuvados, 
Ba 
ue > HO - bii cglpi-p) - m ВЕ - 1) - 
Bgi 
2 Mplat + (Rp-Rgj Bg} t Yp v МЕ > Giñig (101) 
Un este caso, 
са = (102: 
с) Yacinientos d^ дола 
Bg 
ме = бр Bg ~ Om c 1) + Wp- Wi (103 
Byi 


Bas orviginalicote en el ya- 


Come so moneion$ antes, 


cimiento, G, con el fin de poder calcular We. 


Las cucaciones anteriores son generales y consideran inyección 
de agua, ga5,0 ambos, según el caso y lo mismo que compeesibilidad del 
agua y de la formación. Si alguno o algunos de estos términos no exis- 
ten o se concideran iguales a cero, se eliminana de las ecuacionnas 20- 


pr 


spondientes. Recuérdose que la saturación de agua connata, Sy, 38 


considera constante е inmóvil. 


Aparte del uso de balance de materiales para la devenninación 
de Чо, se han desarrollado otras fórmulas de acuerdo al comportamiento 


de presión con tiempo en el contacto agua-petróleo. ба iu 


agua а un yaciniento, generalmente puede dar lugar а tras tipos 
puestas de presión con tiempo: flujo continuo, flujo continu: -n difica 
do y flujo no-continuo. El método de analizar cada caso, pava Y cono- 


cido, se describe чп las siguientes secciones. 


{intrusion 


o Continuo, - Este tipo de intrusión зз expresa 
mediante la fórmuia de Schilthuisd6, donde la rata de intrusión es fu 


ón lineal de la caída de presión en el contacto agua-petróleo, 


дне 
ын = Cg(pi-p) (104) 
dt 
Integrando, 
t 
We = Cg (pi-p)dt = Cg £f(p,t) (105) 
о 


conde. Cg es la constante де intrusión de agua да Sceuilthuis, goneral=- 


mente expresada en bis/tiempo/lpc (donde días, ma 


ses tpiaestees, etc); f(p,t) es el valor del integral; pi 95 la pra 
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sión inicial (original) en el contacto agua-petróleo y p la presión al 


Como puede observarse, para calcular We, debe conocerse la his 
toria de presión en el contacto agua-petróleo(p como función de t) y 
la constante Cg. А su vez para concer esta constante es necesario ћа- 
cer uso de la historia del yaziaionto y calcular We a varios  tienpos 
por la ecuación de balance de materiales, Una vez calculado We se ha- 
ce uso de la Ec.(104) para calcular Cs, escribiendo dic/dt en forma 


de diferencias finitas, 


àt 
Салон (106) 


(рі-р) 


Goneralicute ol st es constante y en tal caso, в. Ade se tema сото 


с cl We al tyempo posterior menos al Че 


promedio de la difevencia ent 


al ti 


mpo antarior, 


1 
Аме = (Wej - Меј.) (107) 


donde j es el período a la cual se desea determinar el valor de Cy. 


22 


En este Forma se colcula una sevie de valores de Cs (un valg 


da tiempo) y de acuerdo а su variación puede concluirse respecto al 


Cs es más о menos constante con tiempo, puede 


їрє de intrusión: 


considerarse que La inirus ón sigue condiciones. de. flujo continuo; por 


otro lado, si varía mucho con tiempo 50 concluye que la intrasión sigue 


otro tipode flujo (continuo-modificado, о no-contiauo). También дође 


tenerse en cuenta que puede presentarse el caso de obtener una varia - 


cic 


n 
о 
з 


siderable al comienzo de la intrusión, donde ocurre un flujo ni 


tamente no-continuo, poro Juego que el radio llegue а un límite (si el 
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acuífero es finito), el flujo cambia а semi-continuo y el valor da Cs, 
ticunde а ser constante. En este caso, puede tomarse un valor promedio 
de Cs en la parte donde presenta valores más o menos constantes y usar 
tal valor para calcular We a tiempos futuros, según la Ка. (105). Por 
supuesto, debe asumirse la presión al tiempo futuro, cosa que puedo ha 
cerse aproximadamente, si la predicción no es para un tiempo muy largo, 
extrapolando la curva conocida, de presión como función de tiempo, Со, 
mo puede verse, este mátodo de predicción no es muy exacto y sólo рле- 


de tomarse como una aproximación de lo que puede ezperarso, Más ada - 


ate se discutirán 


s mejores de predicción. 


Ll eilculo del iategral f(p,t) de la Ec.(105), puede culenlar 


se fácilmenie de la uistovia del yacimiento, de la cual зе сооосоп una 


serias du p iones ру, Di» P2, Pa Pu a les tiempos 0, t4, 52» t3 e 


... tp coPsesponálontas. па esta forma a cada tiempo habrá ема рисо1би 
diferoncial Ар, que puede definirse de varias maneras, tal como 56 
explica en 2} Apéndice C, donde se deducen una serie de Fórmulas para 


el cálculo del integeal f(p,t), teniendo en cuenta, edenás Ја la defi- 


nicióu de Ap,si los intervalos de tiempo АЧ) = ( 


-t 5-1) en que se ha 


dividido la historia del yacimicuto, varian o tienen UN veloce ^ens- 


tante. 


En general, el integral £lp,t) de 2а Ec.(105) puede calouler- 


зе por, 


n 
f(p,t) = == E (2р:)-017- pj-1)5tj (308 


Cuando 


~ 181 = 
tiempo son iguales, At, = constante = At, 
At f n-1 
i(p,t) = == (2n-1)pi - ра - 2 Е Pj (116 
2 441 


Otras ecuaciones para el cálculo de £(p,t) se presentan en ol Apándica 


C. 


Intrusión con Flujo Continuo-Modifí También se conoce como 


el método de Hurst Modificado3?, En este caso, la тата de intrusión, 
además de ser proporcional а la caida de presión, (pi-p) es tambián 


función de tiempo. En forma simplificada, se expresa por, 


Me Cnípi-p) 
а s а ió (111 
dt log(at) 
integrando, 
t (pi-p)dt t (рї-рдк 
рї-р pi-P E 
ме = бр а я бујне = Ср £(p,*) (112 
log(at) jo at tlog(t) 


donde Су es іа constante de intrusión de agua de Hurst, gener: nte ex 
presada ев bls/tiempo/lpe (donde tiempo puede ser dias meses, teimestes! 
etc); f(p,t) es el valor del intogral;a (o а' + logí(3)) es una constan 


te que depende de las unidades de tiempo; pj Y P fueron derinídos en L 


Ec.(105). 


Como puede observarse, 21 este caso se requiere el conocimien 
to de dos constantes, Cp ya (о a'). Para su determinación, se haca us 
de los resultados obtenidos en el caso antucior. Si Cs varía ^n ta 
forma que la intrusión de agua no puede representarse por la iórmula д 


Schilthuis, tal iuformacióa (Cu como función t) se usa рага la determ 
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nación de Cp y a. 


De la Ec.(111) puede escribirse, 


ам 
ае. 
at n Ch 
AO (113 
(pi-p) logat) а' + Лод (х) 


pero el lado izquierdo de esta ecuación es igual al valor de Cg da la 


Ec.(104), Рог tanto, 
Cg log(a) + Са log(t) = Cp (114 
Cy a! + Су log(t) = Ch (11 


Ya que se Съ se сопосо сото función de t, los mejoras valoves de las 
constantes Cp y à pueden determinaese por ol método de 108 mínimos cua. ~ 
drados98, Aplicando la Ec.(115) a cada uno de los tiempos, tj, рага 


los que se han calenlado, Сеј, y sumando talos ecuaciones, se obtiane la 


Лапоміпада primera ecuación normal, 


log(ti)) = ис; (11 


Multiplicando cada una de јаз ecuaciones incluidas en la Ес.(115) por 
el tj correspondiente y sumando como antes, se obtiene la denominada se 


gunda ecuación normal, 


n 
a! Е (Caj tj} + E (C54 tj шигээ 1 ts иг 


Resolviendo simultaneamente las Ecs(115) у (117) se puede 4 


y at ( у por tanto а), 
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n n n n 
E Cs; tj EET log(ti) - : Са) : ©з} log(t3) 
Cn = — I — — (11 
n n n 
n X Cg. t ~ E Cgj E tj 
à CS و‎ E 


Una vez obtenidos los valores de Су y a! (o a) se regresa a la Ec.C112) 
y 5e calcula Me, conociendo p como funcióa de t de да histopria dul улс 
cimiento. Si los valores de Ма caleulzdos сп esta forma comparan favo 
rablemente con los 72010800 de Ме calculados por balance do materiales, 
la Ec.(112) озрговеста adecuadanoate la intrusión de agua al yacimica- 
to y podria culenlaraa узадо We en el fuiuro, siempre y cuando зо 257 


sea más 


похса р al tiempo futuro. En lugar de calenlar W2, y quiz 


práctico, se puede vorificar si las constantes Cp y a! 
cuadamente la lutrusión, calculando la prata de iatrusión, ез dWe'dt, 


por la Ec.(111). Si esie valor compara Favorablenante cou la vata de 


intvusión calculada рог balance de materiales, al calcular Cy, En. 108) 


los valores Cp y а! (у а) зе consideran satisfactorios y Ve pueda cal- 
cularse a partir de la Ec.(112), Si lo» valores de со W^ no comparan 


favorablementa, la intrusión dabe expresarse en forma diferentes» 
А р 


Para el cálculo de We, 22 necesario determinar el integral de 
la Ec.(112). El Apéndica C incluyo ia deducción de iórmulas para el 


el cálculo de este in! ural, de donde se han tomado las siguientes, 
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1 (рї-ра) ЭГ n-1 (pi-p3) 
Elp) = —|—— —— (tg-tg.1) + 2 e tp rigor) (126 
2 log(atg) j=1 1og(atj) 
donde tg = 0. Si los períodos de tiempo son iguales, Абу = constanta * 
At, 
(pi-pg) n-1 (pi-pj) 
£p, t)se РАК | emo г E (1? 
log(atg) 151 log(at;) 


usión con Flujo Wo-Continuo.- Generalmente se conoce como el 


mútodo de van Everdingen y Hurst. Ка este caso, la intrusión total de 


agua а un yacimiento se expresa por las ecuaciones discutidas anterior 


mante, с5.(23) y (95), para aculferos lineales y radiales respectivas 


monte, 
19397 n 
ме = — ээ E Ар] (8: 
8,615 141 
9 n А 
Wo = 1,1 8фелАД === E Ар) 41СЕ6:54-170) (9 
3609  js1 
donde ty > 0. En general, estas eovaciones pueden escribirse, 
Че = ^pj 01 65::54-12) = Cy f(p,t) (i2 
donde la constante Cy (denominada comunmente в), está definida por: 
Acuíferos Lineales, 
$bhLe 
Cy = oM — barrilecs/lpc. (1: 
5,615 


Acuíferos Radiales, 
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Cy = 1,1194chR? 2440. baeriles/lpe (129) 


y 3609 
donde с es Ха compresibilidad efectiva dal acuífero, es decir, incluye 
además de la сотргев15111444 del agua, la compvesibilidad de la formas 


ción (с=сү+с„). Si el acuífero es lincal Q(-tj 0) en (Cra 714-15] 


En la deducción de las Ecs.(93) y (95) se hizo referencia а 


como función 


la Fig.55, donde aparece la historia da presión (pres 
de tiempo) como una linea continua y са forma escalonada, пой el fin 


de poder aplicar al principio de superposición disczutido. Ка la prác 


tica, se conocen una sopie de pensiones ру» Ру» P2 +з, Ра а los tiom- 


pos 0, t4, t9 s.. tne El mátodo más conúa y posiblenenta el міз usa- 


do para obtener los tecrementos de presión, Ар» es el propuesto рог 
van Everdingon y оџеоз30, anunciado en la siguiente forma: la caída de 


presión a cualquier tienpo, pueda tomarse como la mitad бе la caída de 
presión еп el intervalo antevior, más Ja mitad de la caída de presión 


amente, para la prinevca caida de pre 


en el intervalo posterior. Mi 


sión, а mitad dol intcevalo ancerion será cero. En forma analítica: 


1 A 
Ару = 0 + e (psp)? PD 
2 
1 
apa = (pispa) + --- (рА-Ра) = pido) 
2 2 
1 i 1 
^ пољ - q л -pa по-1- ET 
Pa > (24-р) 5 (ropa) y Pa 22) 
1 1 1 
Any = === (озера) + =e (pa pi) = lpas Py) 
2 2 2 
. 


1 1 1 
POE BR Р Aldo. -р.) = (ps.97P4 
Ap; 5 (р 4-2 P 4-1) ив (0j 17? ) : ОР 1-2 ру) 


'ор0џо102 


DEE 
зоо vn ep upisaıc эр 


410616: 
NI 


ч a лт 


1 1 1 
Арр = тү Pn-27Pn-1) ж Ту PaPa)? = тү (Ра-а“Ра) (125) 


La Fig.55 ilustra las ecuaciones anteriores. 


Los valores de Apis también pueden obtenerse "a ojo", para lo 
cual se gráfica p como función de t y se traza una curva continua. Lue 
go, "a ojo" se toman los Api de acuerdo a la forma de la curva tratan 
do de compensar las pequeñas áreas triangulares que quedan por debajo 
y por encima de la curva al tomar cada ёру. Finalmente, se leen nue- 


vos valores de py y t; (j = 1, 2, 3, ... n). 
Pj i 


La obtención de Fy(tp) у Q(t) fué ampliamente discutida, йе 
acuerdo al tamafio y condiciones iniciales y de límite impuestas al a - 
cuífero. Recuérdese que para acuíferos lineales y acuíferos radia - 
les infinitos, los valores de Fx(tp) y Q(tp) se han determinado como 
función de tp, Tablas I, IV y Figs. 10, 26 y 27; en canbio para acuífe 
ros radiales finitos (limitados), Q(tp) es función de tp y rep» Tablas 
VII y Figs. 35 a 39. Por tanto, la aplicación de la Ec.(122) exige 
conocimiento de: (1) tipo de acuífero: lineal o radial. (2) tamafio: in 


finito o finito, (3) si finito: acuífero cerrado o acuífero comunicado 


a una fuente de abastecimiento de agua. Por otro lado, 1а 
del tiempo adimensional, tp, requiere el conocimiento de una constan- 


te, Атр, para converitr tiempos reales a tiempos adimensionales: 


Acuiferos Lineales, 


0,006323kt 
tp о 5 Айр t (126) 
$uci? 


Acuíferos Radiales, 


м es جک‎ = + 


a 148 ж 


0,095323kt . 
pos ee emos Аср t (121 
фус RÈ 
i 
Como puede observarse, tanto la constante Cy como Атр requie 
реп datos del acuífero: $, ©, К, (o L), valores comunes а Cy y А. Pa 


ra ecuiforos radiales, a medida que el producto écR2' aumenta, Atp dis- 
i Н 


minuye у Су aumente. Por tanto, lopos que definen а Ар» parcial- 


mente definen а Cy. Aden , los valores k yu se encuentran en Ар,» 


ућу 6 (о b) за encuentran en Cy. һа mayoría de esta información se 
conoce con poca exactitud o es necosaeio agumirso, a no ser que se per 


fore un pozo dentro del acuífero, cou lo cual puede dotorminarsa parte 


de tal información con т jov exactitud. Generalmente, фу k y h 


man igual a los valores correspondientes promedios para la form: 


la zona de petrtleo; Y y Ca se toman como 108 valeres d:l agua да for- 


1» 111444 de 


mación y mucnas vaces no зе toma on cuenta la compras 


moción del acuífero; R 0, b y Lb (este último valor рака ej caso de 


couniferos 14164422), se obtienen de acuerdo al tipo de estructura уат 


о 


inic 


ato-acuífero y del conociemiento del tamaño y forma del mismo» 


En esta forma, Су Y Ар» pueden calcularse de infovcación conocida, pg 


vo tales valores deben tomarse como "valores guias", ya que ics "valo- 


ros verdaderos" es посол: detevainarlos de 1а historia del yacimien 


to, es decir, Су у бер зоп incógnitas en la ре. (122), Además, si el a 


3 = 


cuifevo es radial y finito, debe conocersn el valor de гор» Con el fin 


de evaluar Q(tp). 


Para probar si Че de la historia 441 yacimiceto „Ме calculado 


por balance de materiales, en función de persión en el contacto agua - 


petróleo y tiempo), puede eopreseniarse por la Ec.(122), es necesario 


hacerse por ensayo у ёггог (tanteo): en caso de que ој acuífero pueda 
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considerarse infinito, 89 prueban diferentes valore с Aty y si ста + 
cuifero es finito (limitado) deben probarse diferentes valores de Сор» 
y para cada uno de estos valores, probarse valores de Абр. La exacti- 


tud de estas constantes se vorifica por la consistencia dal valor de 


Си calculado por, а; 


(Balance de Materiales) 


ЭР” 128 
Cy = AT (129) 
I Ару 0154-51-10) 
j=1 


Si esta razón permanece más o menos constante para detovainados valores 


de Atp y Гор» 28 porque Че puede expresarse pov la Ec.(122) y el mejor 
valor de C, puede tomarse como un promedio de los valores calculados 


con la Ec.(126). 


El procedimivunto de análisis más recomendado cs el siguiente: 


Si Сор puede sonsiderarso infinito o determinarse cou suficiuute 


citud, зе calcula un valor aprozjuado de Atp con los datos conocidos 


(o conocidos y aucmidos). Luego ва toman valores mayores y Menores de 


se calcula n valores du C, de la 


^tp у pava сада uno do зыр s valoras 


ya.(128). Se “tetesminan un Cy promedio para сада valor de ауу 


1 n 
Cyg com X vj (79) 
n jal 


y са calcula Ja desviación normal vcorrespondicnte, 


1 n 
СЕН (130) 
ү ур] vj 


¿ión nownal como función de los valo- 


Finalnente, 5e a la desvi 


res de Аїр considerados y donde ocurra un mínimo on la g, se toma co- 
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Mejor Valor de Aty 


Constante de Conversión de про, 4р 


Pig.56 ~ Geüfico esperado de desviación normal, Y со- 
mo Función de la constante de conve, 
Аку, para determinar el mejor valor de Atp. 


Гага сада curva, Кор = constante. 


mo el jov valor 43 Аїр» tal como ilustea la Fig.36. А veces pue- 


de ser más práctico gráficar о como función de Log (At y 


‚ ya que puede 


nimo de 9» 


visualizarse mejor el m 


Con el mejor valor de Ар obtenido, 59 determina el Сур со- 


rrespondivate quo sepi шоог valor de Cy. Si no se conoces rap: Se 


у рата сада ray asumido se hacen los cálcu- 


asumen diforontes valore 
los indicados antes. Finaluente се tendrá un gráfico de q como fun - 


curva pura cada valor de 


dorados, un 


ción de los valoras da Atp con 


Үер asumido. Та curva que presenta el mínino valor de 0, sovi la dul 
mejor Үүр У en el punto de mínimo se toma el ^tp correspondiente como 


ol mejor valor de tp: Con este valor se regresa y 898 calcula al Сур 


correspondiente, que será el mejor valor de Су. 


so Desconoce el Patróleo o C 


En este caso, habed una incógnita más, о 564 ol valor de N. 


Como eu el caso anterior, s? hace uso de la ecuación da balanco da па 


tesialas, pero en forna ditoronta. Esta ecuación so ordona en taJ for 
ma de agrupar los tórminos tal quo se obtenga la ecuación de una linca 


vecta. ha pendiente de tal linea recta es la constante 42 intrusión, 


C5, Cp о Су y ul intercepto, el valor de N. Las variables son cauti- 


dades que puedan eulcularse de datos de presión y de producción. Una vez 


inada la fórmula por la que vá 4 exprosarse la intrusión de чума 


tal que se cumpla la lincavi id de la ecuación, los mejores valor>s de 


за coustante de intrusión y del peiráleo originalmente en el yacimien- 


: ЯА 3: 
to, se cbitionen por minimos cuadrados ج‎ 


и de la ecuación general de balance de averia- 


128, Гс.(101), puec: 


прве + Срба) nj] t Wp ~ Ме - Wi - бу Pig Teen 


Moso mecene neme: - ————— M————— Á— (131 
ps 
(0,1 4) жээ > - 1 + Bej celpi-p) a, 
Bas à ко : өм 
ге Swe 4 A Ce cA 
SONA р 
(1-75) 


donde e, está dado por la Ec.(102). Si se defines como D, el 


дор de la ecuación anterior, puede escribirse, 


Np [Be + (т„-к„{ В] АЛ CiBig Ме 
= MERE TENE 


D D 
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se define, 


шин (Rp=RsidBg) + "p - "i 7 Ci Pig 


Y = мы (133. 


р 
y el valor de Че se expresa por, 
We = C f(p,t) (134 
donde C puede ser Cs, Су о C, y f(p,t) son los valores correspondien 


tes al integral en las Ecs.(105) y (112) o a la sumatoria de la Ес. (122) 


Por último si se define, 


X = E(p,t)/D . (135 


La Кс.(132) рисле escribirse, 


үст OX FN (136 


Para yacimientos de petróleo no-saturados, үсѕ.(97) y (100), 


puedo escribirse, 


ORARE Че 


м а d 


(Bo-Boi) (Во-во: ) 


{13 


donde с. está dado рог la Ec.(99). Si We s? reemplaza por la Ге, 4139) 


y se define, 


н 


— M (13 
Bog Celpi-p) 


f(p,t) 
— (14 


Boi Seofpi-p) 
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Las Пс8.(137) y (126), pueden expresarse рог la Ec.(136). 


Finalmente, para yacimientos de gas, Ec.(103), 


G= (141) 
Si зо veemplaza Wẹ por la Ес.(138) y definiendo, 
Y = (142) 
f(p,t 
X = не (143) 
в, 
бат 
Bai 
La Рс. (191) puede escribirsa, 
| у = CX + 6 (1457 


Esta ecuación es similar a ја Es.(135), excepto que aparece 0 en lugar 
de-N, Por tanto, cualquiera que 524 de tipo de yacimiento, la ecuación: 
de balance de materiales pueda expresarse рог medio de una lin«sa recta 
Los valores de Су Но С von constantes; el valor de Y es función ded 
tos de producción : Мр» Gp: Rp, ЦЭН de datos de inyeccióniWi, бр» Big 
de presión (datos PVT) : Dti, Boi, Dt, By Вр Rgj; del tamaño de la c 


p 
ра de gas, m, de la saturación de agua connata, Sy y de las compresihi 


lidades de la formación y del agua de formación, Cp Y Си» ®" la попе 


de petróleo. A SU vez, producción y presión son funciones de 
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Рос tanto, Y es función de tiempo. El valor 4e X depende del tamaño 
de la capa de gas, saturacion de agua, compresibilidades, presión 
(datos PVT), f(p,t), el que a su voz depende tiempo y datos especifi- 
cos de acusrdo al tipo de intrusión, según se discutió anteriormente. 
En esta forma, para cada tiempo de la historia del yacimiento pueden 
calcularse valores de Y y Х. Aplicando estos valores а la Ec.(136) o 


(144) pueden calcularse los valores de Су Мо б, 


Al igual que en el caso en que N es conocido, a continuación 
se presentdá un analísis para la determinación de C y N, de acuerdo а} 


tipo de intrusión. 


Intrusióu con llujo Continuo.- Ба intrusión de agua so expresa 
рог la Ec.(105). El valor del integral f(p,t) de esta ecuación, se 


calcula a partie de la historia de presión por las Ec.(108) a (110) u 


otras ecuaciones, tal como 


explica en el Apéndice C. А cada tiempo 
ас la historia se calcula un valor de Y y X. El mejor valor de C (са) 
N o G, зз ohbtienoan aplicando mínimos cuadrados, según las ecuaciones 


deducidas en el Apíndico D, 


с = (145) 
+ n Y 

N (o G) = У - C X Xi 196 

n ki 3 1 3 : ) 


En lugar de la Be.(115) puede usarse la Ec.($D) del Apéndices D. Аде- 
más, se recomienda al calcular los mejores valores de C y Нос, tos 


mar el promedio de los valores calculados aplicando minimos cuadvados, 


primero considerando linea de regresión de Y sobre X y luego línea de 


regresión de X sobre Y, sogún ecuaciones incluidas en el Apéndice D. 


Una vez calculado el mejor valor de C y Но С se regresa y se 
calcula Ме por balance de materiales у рог la Ec.(105) de Schil:huis y 
se comparan los valores. Si comparan favorablemente dentro de los lí- 
mites de error pernisibles, la intrusión de agua sigue flujo continuo 
y los valores de Су Мо С calculados representan рага el yacimiento, 
la constante de intrusión de agua de Schilthuis y el valor de petrólco 
o gas originalmente en el yacimiento, Si los valores de We varian mu- 


cho uno del otro, la intrusión de agua no sigue flujo continuo. 


El análisis de este tipo de inteusión puede hacerse си Forma 
diferente y quizás sea la forma más práctica, especialuente si los cál 
culos se hacen manualmente (ealenladoras. de escritorio). Es el máto- 
do до Graficaciones Sucesivas, sugerido por Havlena y одвећ“. Una vez 
calculados los valores de Y y X se gráfican en papel normal. En «esta 
forma se imparte un significado dinámico а cada punto calculado, уа que 
dependiendo de la forma del gráfico rosultante, puede concluirse respe 
to al yacimiento (tipo, tamaño) y а los valores que doben vaniarse par 
continuar el análisis. Si el grafico rosultante es una linca pocta о 
muy cercana a ella, se calculan los mejores valores de C y No С рог m 
nimos cuadrados y se chequean con los valores obtenidos del gr fico, y 


que la pendiente es igual a C y ol intercepto es igual a M o б. 51 10 


puntos se apartan mucho de una linea recta, la intrusión no sigue #19] 


continuo y dube considerarse otro tipo de intrusión. 


A veces cl acuífero es muy pequeño y la intrusión de ayua pue 


1 
de expresarse por 04 


We = C'(pi-p) (147) 


donde 0! = We, y W = volumen de agua cn el acuífero. En este caso, 
para el cálculo de X se usa (pi-p) en lugar de f(p,t) y cl gráfico de 
Y contra X se analiza en forma similar. Si ocurre lincaridad, зе de- 
termina C' y N о С por mínimos cuadrados y se compara con los valores 
obtenidos del gráfico. Del valor de C' puede detesminarse el volumen 


de agua del acuífero, W. 


El sentido de variación de Y y X eu este caso ocuvre en for- 


>> 


па diferente. El valor de X disminuye а medida que la presión disuinu 


ус, 9 sea, a medida que el tiempo aumenta. ato es debido а que el de 


nominador de Y y X aumenta más rápido que (pi-p) a medida que р dismi 


ivo (pendiente de la linea recta), 


риуе y ya que C' es siempre posit 
luego Y debe disminuir a medida que la presión disminuye. En esta for 
¿fico varían en sentido inverso, o sea que Y y X 


ma los puntos del gra 


disminuyen a medida que p disminuye (tiempo aumsata). 


De lo anterior puede concluirs» que el acuífero es muy peque 
fio y se calcula N o ба росїїт do la acuación dc balance de mateviales 


(че se considera igual a сего), el valor aparento de N o б calculado 


(realmente so calcula el valor de Y) disminuye con tiempo. En la prás 
tica puede оситргіг que el comienzo de la intrusión el compertamisato 
dol acuífero siga flujo no-continuo, la variación de Y y X sea nocmal, 
os decir, aumentan con tienpo, y luego, si cumple la Ес. (147), los va- 


loves de Y y X disminuyen con tiempo, El punto de cambio de pendianto 


ocurre al cambiar de flujo no-continuo a continuo, 


odi ficado. La inteusión se expres 


jo Continuo 


por la Ёс.(112) y para el cálculo 401 integral t(p,t) de esta ecuación 


además de requerir la historia de рге 
lor de la constante de conversión de tiempo 


calculado es el siguiente: 


(1) 


- Se calculan valores 


historia del yacimiento, Есв.(133), 


(2 


en la forma indicada más adelanto) 


lores do X para los mismos tiempos еп que зе 


(3) - £a vopite 
(4) - Para 


de Су Y o G por minimos 


(5) - Раса coda valor de а, 


N о G obtenidos си (^) y los valores 
una serie de valores 


calculau de Үс 


(6) - Pura cada valor de а 
(1) y Үс calculados en (5), se determina За 


contaje de desviación normal, % 9, 


Para un valor considerado (determinado 


de la солзїспїё a , ye 


el Paso (2) para otros valo» 


cuadendos, Ясз.(115) y (156). 


con los 


por la Ес.(130), o 


y los vi 


es necesario conocer el va- 


a . El Procedimiento de 


de Y рага los tiempos en que existe 


(139) o (1^2). 


aproximadamente 


calculan va 


calculó Y. 


cada valor de а , 59 ootionan los w2iores valores 


урћадјав D, 


mejores valores 


46-0 y 


de X calculados on (2) y (3), se 


(унэ. 


pos de Y calculados en 


isviaeión nocmal, y 01 por 


о = (148) 
! n 2 
ha = m —— } (159) 
T. 
(7) - Se geafican los valores deg (о $a) соло Zuncióu de а. 


El valor de а 


correspondiente а Ма desviación normal ма, es el 
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mejor valor de а para el acuífero. Dependiendo del iutervalo de 


variación de a , puede graficarse el valor de а! (logía)) como función 
de o, obteniéndose un mínimo mejor definido. Sea cual fuere el gráfi- 


co, la forma de la curva resultante es del tipo ilustrado en la Fig.58$. 


Para determinar los mejores valores de C y N о G, зо gráfican 


Ф 


stos valores como función de о у donde ocurra el mínimo valor de о, зе 
leen los mejores valores de C (Cp) y N o G. Es también recomendado grá 
ficar C y H o 6 contra tiempo real у analizar su variación. Los mejo- 


res valores de C y N o G deben presentar la minima variación con tinm- 


po. 


(ву Finalmente, con el mejor valor de Но G se calcula Y, por 
balance de materiales y con los mejores valores de Ch y a (о а!) se 
calcula también И, por la Ec.(112) de llurst Modificada. 51 estos valo 
ros chequean satisfactoriamante, la intrusión de agua puedo represen - 
tarse por la Ec.(112); de Lo contrario, la intrusión no sigue condicio 


nes de ¿lujo continuo-modificado. 


La primera aproximación del valor de а , puede obtenrse a 


partir de los resultados al analizar la historia del yacimiento рага 
flujo continuo. De este análisis se toma el valor de Но Су se caleu 


la el valor de а рог el método discutido en caso de conocer N o G. 


Con valores de a mayores y menores al obtenido en esta forma, se cale 


la el Paso (3). 


El caso de flujo continuo modificado, también puede analizar- 
se por graficaciones sucesivas de Havlena y Odeh. Se calculan los Pa- 
sos (1) (2) del procedimiento anterior, se gráfica Y (ordunada) como 

р » 5 


Función de X y se analiza la curva resultante. Si indica una linea re, 


= 459. = 


ta,se calcula С (Cp) у М por mínimos cuadrados para el valor correspon 


diente de а. 
muy grande y si 
queño, Fig.57. 


el Paso (2), se 


Si presenta curvatura hacia abajo, el valor de а 
la curvatura es hacia arriba el valor de а cx muy pe 
En esta forma, modificando el valor de а , se repite 


gráfica de nuevo у asi sucesivamente hasta obtener la 


mejor linea recta. Finalmente se obtienen los mejores valores da 0 


(Cp) y No б y se prueban si estos valores representan adecuadamente el 


acuífero, según el Paso (8). Puede suceder que el gráfico resultante 


no presente uno 


prichosa. Tal s 


de los tipos de curva mencionados, sino otra forma са 


ituación puede deberse a uno o varias de Jas siguientes 


vazoanes: (1) datos огга405, (2) cálculos errados y (3) el acuífero no 


se comporta según este tipo de intrisión. En tal caso se recomienda a 


nalízar la historia considerando flujo no-continuo. 


En general, este mótodo de flujo continuo-modillcado es poco 


usado en la ргабїїса y directamente en el análisis de un yacimiento con 


sideran flujo no-continuo. Sin embargo, ya que la intrusión Ме, ca 


función de presión y de tiempo, en este aspecto se asemeja el flujo no- 


continuo, discutido а continuación. Adomás es lógico suponar, que si 


un yacimicnto cumple flujo continuo-nodificado, con mayor razón ве ёс» 


pera que cumpla 


flujo no-continuo. Esta puede sev una d? las (бомгаз 


рог las cuales no sea muy usado en la práctica, 


Intrusión con 


flujo ^Continuo.- Este caso es bastante similar 


al caso anterior ya que la sumatoria f(p,t) de la Бе, (122), ¿iguala 


te depende de una constante de conversión de tiempo Ар según se di 


ü- 


cutió anteriormente. Además el valor de f(p,t) depende del radio udi 


monsio p . 
onal Pap 


EL método más comín de análisis de este caso se debe 


a van Everdingen, Timmerman y МсМаћоп35, y generalmente se conoce со 
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Fig.57 - Curvas típicas de Y como función de X pura di 
ferentes valores de а en el análisis de intru 
sión de agua con flujo continuo-modificado., 


mo método de VTM. El mátodo sugerido рог Havlena у Odeh"? es tambien 


41 


muy usado. Melluen y Stanley”? presentan formas diferentes de аай 


tarse mejor un método que otro, los métodos mencionados 


dos comenzando por el método УТМ, 


ado УТМ.- Este mútodo sa adaptó al análisis del caso anterior de 


flujo continuo-modificado, donde f(p,t) solo depende de una constante 
de tiempo, a (о at). En cazo de flujo no-continuo, f(p,t) dependerá so 


ов, 


lo de Atp en acuíferos infinitos y además de бор en acutfevos fi 


E 5 


Si за trata do acuíferos infinitos donde f(p,t) depende solo 
de Atp, este valor se calcula aproximadamente con la mejor información 
conocida del campo de acuerdo a su definición, Ecs,(125) у (127).5е to 
man valores de АЁ monoras у mayores al valor calculado. La variación 
de Atp depende del criterio individual y de la exactitud de la informa 
ción del compo usada para el cálculo inicial de Atp. Una vez seleccio- 
nados los valores de Atp, зе cfectuan los Pasos (1) al (7) del procedi 


miento explicado en el caso autevior,tomando Аїр en lugar de а, y 


por supuesto, el v.iov de C obtenido, será la coustanto de intrusión de 


zan Eveordin:con Hurst, Cy, El Paso (8) se moditica en la siguiente 
y » v E 


форма: 


(8) - Los sojoras ves de Atp y C, obtenidos en el Paso 


autevrior, priden сомигорагаа nultiplicándolos entve si y comparando es 
to producto coa el valor obtenido de datos de compo de acusrdo а rus 
Actinieiones, Ec. (172) кЕс. (127), 

(1,119)(0,006323/xh00 


Cary (Су) = om — ———— (150) 
360 y 


(19 ,65х1079 кв 92 
CAM MV) = mm سسس‎ (1 
р 


Si kh09/p 5e соцоса con un buen grado de exactitud, esta хе 


lación contribuye x un chequeo de Айр y Су. 


ЗЕ по осцоо error en los cáleulos y la infoecrmación ns ге- 


presontativa dol compo, se considera пис los valores do las constantes 


Му y Cy Topersontun adecusdamente la intrusión de agu 


sie (ajustar), si es necesario, los datos de рте- 


sión de la historia del yaciminato, calculando 9, tauto por balance de 
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matoriales, como por la ecuación de van Everdingen y Hurst, vn" tal for- 
ma que los valores de Ма calculados por ambos métodos, varie dentro de 
una aproximación normal. Las presiones corregidas sevän las usadas en 
el resto de los cálculos del yacimieuto. £s entendido que si la intru- 

sión calculada por ambos métodos varía demasiado, puede ozuppir que la 


información del yaciniento о los cálculos están errados. 


Si el valor de Хер es desconocido, 08 necesario asumir valo- 


ves de csta constante y para cada valor de rep repetir los Pasos (1) a 


(7) descritos. Para cada valor de r se grafica 0 como función de At 
ер 5 D 


( о loglàtp) у, дуу То G. La curva que presente la mínima desviación 


normal, on este punto se leen los mejores valores de АКр» Cy Y Нос. 
El velor de Pop eorvespondiunto será rl mejor valor de Үер, Finalmente 


el Paco (3) anterior. 


Hoc ena y a forma cimilar 2 los Casos anteriores, la 


intrusión de agua “on flujo no-continuo pusde analizarse por el mítodo 


40 
4% Una vaz calculados 108 


de graticaciones sucesivas de Uavlena y одећ 
valures ne Y y X se grafican sn papel normal. De acuerdo 8 ја forma de 
la сцгуа resultante, puede contiauarse el análisis en forma 165109, La 
Fig.58 presenta dos casos principales que puedan pesultar. Además de 
estos савоз, pueda ocurrir que 195 puc tos fraficados 59 presenten an 
forma desordonida, lo que indica que los cálculos, 1а historia del уйл 
cimiento, O subas cosas, estén ervados. 


si los puntos graficados presenta. una línea vucta, Curva n, 


o 


es porqu las constantes Год» Cy y No С son cobractas. Los mejores УА 
loves de C, y 19 а, pueden obtenerse por mínimos cuadrados y tambien 


del mismo gráfico (pendiente = C. intorcepto = Но 0) para los valores 


de Atp У Үер vorrespondienta5. si los puntos presentan curvatura ћазја 


1 
Valor de N o Q 
X 
11 4.58 + función de Х para 


o5 sumatoria v(p,t) en 
21 análisis de inteusión de agua con lujo 
no-continuo. 


arriba, Curva E, o ha 


abajo, Curva F, indica que el valor de f(p,t) 
es muy pequeño o изу grande pospoctivamente. Esto а su voz, indica que 
cl valor de Yap, Atp, o ambos, sou bajos o altos. Por tanto, deben au- 


T 


mentaesc o disninuiesa según el caso y repetir los cálculos con estos 


nuevos viloros y asi sueesivacente hasta obtener la me 


De nuiva, Cy y Мо G ss calculan por minimos cuadrados y зе 
del gráfico para los valoras de АХ, y Pep correspondiontes. Pinalaeuto, 
Si los puntos presentan nua curvatara en forma do S, Curva G, ¿indica 


que el acuífero puede representarse 


jor pov un sistema lineal. 
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Con el fin de saber el sentido de variación de los valores de 
413 У Гар еп caso de obtener curvatura hacia abajo о hacia arviba, Fig. 
58, debe tenerse presente la forma de las curvas de Q(tp) contra tp, 
Figs.35 a 96 (acuíferos cerrados). Por ejemplo, si зе fija г.р у se va 
ría Аїр, debe observarse que para que intervalo de valores de tg ocurre 
el máximo de variación de Q(tg) con el fin de aumentar o disminuir (р, 
de acuerdo a la curvatura. Por ejemplo, para гер = 15, Pig.39, la máxi- 
ma rata de aumento de Q(ty) ocurre para 100 < to х 500, Q(tp) contra tp 
se va haciendo mas y mas horizontal hasta alcanzar un valor constante 
(flujo continuo) de 112,0 para ty = 1000,0. Si para el Atp considerado, 
los ty caleulados caen en su mayoría por encima de 500, ocurre un aumen 
to muy pequeño de f(p,t) y entonces Atp debe disminuivoe. Lo coutrario 
deba hacerse si la mayoría de los t5 fueron menores de 100, ex decir, 
At debe aumentarse. En resumen, una vez fijo Tep? debe examinarse muy 


cuidadosamente la curva de QU) contra t, correspondiente, con el fin 


de determinar si Ар debe aumentarse o disminuirse para ојљепек Linea- 


ridad en el gráfico de Y contra X. 


La curvatura tambien puede corregirse fijando el valor de áto 
y variando Tep’ Si se aumenta este valor, equivale a aumentar mas räpi- 
damente el valor de Q(tp) - у por consiguiente el valor de f(p,t) = a 


tiempos mayores que a tiempos munore 


» trayendo como cousecuencia la co 
rrección de la curvatura que se presenta hacia arriba. Si зо disminuye 


^ cun Д trario. 
Pop» 99urre lo contrario 


De la discusión anterior puede concluir que varias combina- 
ciones de Atp y rep pucden proporcionar una línea recta y presentar la 


solución correcta. En estos casos, para determinar los mejores valores 


de At, y rop debe hacerse uso de pruebas estadísticas: "prueba de consi: 
D eb р р E 


: А А А 222423 
tencia", explicada mas adelante, e "іпсотуа1о de confideneia" + 


En algunos yacimientos, el gráfico de с como función de ftp 
(o log(Atg)) no presenta un mínimo definido, sino un intervalo de бер 
donde el valor de 9 varía relativamente poco. Esto se debe a que el ya 
aimivato es incensitivo а cambios de Atp y en este caso, là obtención 
de este valor es solo de interés académico. Sin embargo, en algunos са 
sos se recomienda la "prueba de consistencia" 3; зе toman varios valo- 


res de Аїр alrededor del punto de g mínimo; para сада Aty ве grafica 


с N o G como función de tiempo real a cada gráfico, por mínimos 
у Y р , 


cuadrados se determina la mejor linea recta y SU pendiente. stas pen- 


dientes se д 


e" 


sAfican como función de Atp en papel normal. La interse- 
las curvas resultantes, corr .sponde а} valor más probable de 


ricomente, debe ocurrir la intersección pava una pendiente i- 


gual а cero, ya que si los valores de Tap Y Atp 898 correctos, С, y ¥ 


o G nó deben variar con tienpo y la pendiente debe ser igual a cero» 


En el análisis de intrusión de agua, el valor de m cs сего, 


А 4 
se asume o se conoce сол razonable exactitud, Havlena y Odeh |, diacu- 


о 


teu əl caso сп que ol valor de m soa una nueva incógnita y describen 


la forma de resolver el problema. Al mismo tienpo hacen énfasis en que 
este caso roquiere datos de extraordinaria exactitud, principalmente 
datos de presión. Además, el modo requiere conocer derivadas respacto 


a presión de determinados factores, lo que contribuyo aún mas a dispo- 


iva- 


nor de datos de extrema exactitud. Por ejemplo, 598 pequiera la der 


da respecto a presión de (5с-9:1)» o hace que en la 


práctica este mátodo tenga muy poco uso. 


Meto 10. £1 mitodo de van Everdingen,fimmetaan y McMahon + 


método VTM = discutido anteriormente, no tvabaja cuando existo errors 
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en los datos de presión, aún en el caso en que tal error sea pequeño. 
Habuen modifica el mávtodo de VTM, proporcionando resultados satisfacto 
vios, incluyendo el caso en que los datos de presión presentan dudas 
E 
respecto a su exactitud. 
En este método la ecuación de balance de materiales se 


arregla en tal forma de obtener la ecuación de una linea recta que par 


ta del origen, Y = ах. 


Si la Ec.(182), para yacimientos de petróleo saturado , ве 
multiplica рог D/N, cl término Wa se pasa al lado izquierdo, se роет - 


plaza por la Ге.(131) y además se define, 


үзээ (152) 
Xx a, [bt + (арка Ву] + Wp + Wi ~ GiBig (253) 
а = i (154) 
| N 2 
puede escribirse, 
{ = ax (155) 


Para yacimiento de potróleo no-satueados, las Ecs.(137) y (138) pus - 
den arroslerse en Хоста Similar y expresarse en la forma de ia Ec. (1353). 


Eu este caso, 


эссэ ep, (156) 


e 
ne pi-p) fuss Ep) (157) 


Bo + Wp- Yi (158) 


а = -— (igual a Ес.(15%)). 
N 


Finalmente, en forma similar, para yacimientos de gas, la Ec.(141) tam 


bién puede expresarse en la forma de la Ec.(155) En este caso, 


ү = ا‎ ~ | + Евр, (159) 
Bai G 
Y = OgBg + Wy - Yi (160) 


"EM (161) 


Como puede observarse de las ecuaciones anteriores el valor 


de Y contiene los términos con posibilidad de error mayor, especialuen 


to si el error ocuere on la presión, ya que todos 5 términos, excep- 
to C/N о С/С y m (сп caso de que exista capa de раз inicial) son Fun- 
ciones de este valor. Adewîs contiene térmiuos que son diferencias en 
factores voluméiricos, muy senciblos a presión. Por otro lado, X con- 


tiene términos de producción (e inyección оп caso de que axista), gene 


ralmente conocidos con mayor exactitud, 


El valor do C/U para el caso de flujo no-continuo, puede cai 


cularse aproximadamoate si зе conocen los valores de  conpresibilidad, 


saturación de agua connata y factor volumétrico inic del petróleo. 
Reemplazando a Cy de su definición, Ecs(123) y (124) y а N de la Ec. 


(96): 


Acuíferos Radi 


Зее: 09 
аа €— (162) 
200 (x-Sy) 


Acuíferos Lineales, 
с сеВоз (153) 
— = e хөжчөлэлэ 5 
u (1-8,) 


Pava yacimientos de gus pueden deducivss ecuaciones similasea paca C/G. 


Basta cambiar Во; por Виз, 


El procedimicato de cíleulo puedo resumir en log siguinates 


pasos: 


(1)-5е calcula ul valor de C/il o С/С de una de las осзас 


anteriores, Si se dispone de datos suficientes, se calcula аррохімаца - 


mente un valor de la consinte атр poe la Es.(128) о (127); се la поа- 


tario debe asumirse. Pl valor de гүр ве сепосе о пе asume, Coa los м 


lores de C/ o C/G у бебу, рага un determinado Pops 59 calcsliaa asa / 227 


miento (un valos 


ріс de valores de Y y X de la historia dol yas 


da tiompo), usando las Nes.(152) y (153), Есз.(150) y (i187) © 
y (160), de acusrdo al tipo de yacimientos. Aplicando mínisos cusdpa 


dos, so calcula el valor de a, 


n 
EX Yi.) 
1 1 153 
а = e Q man = т (154) 
N G n 
rou 


1 


12)-Con el valor de a obtenido del paso autueior y para cad 


valor de X se careto. "n valor de Ус, de la ecuación Y = ах. Con es 


tos valores y los valores du Y oves “Хав en ol paso sateriov, 59 cal- 
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cula la "Suwa Normalizada" del cuadrado de las diferencias de los valo 


res de Y, dada poo, ‹ 


1 
5508 по (165) 


donde Ymax es el valor máximo de Y calculado en e) paso anterior. 


(3)-Se asumen otros valores de Аср (valores mayores y menor 
res del valor calculado aproximadamente) y se repiten los Pasos 1) y 
2). Se grafica SSPp y hm (Но G) como función de Ар. Donde ocurra 
ol mínimo en el valor de iSDa se locen los mejores valores de Atp y N o 


С para el rap considerado. 


("n)-2i r es desconocido, se asumen diferentes valores de 
eD 
гор y sẹ vepitoa los Pasos 1) a 3). Па curva que presenta el $501 ni- 
nino en este punto so toma el mejor valor de 165, Но C HEP corres 
P 3 р» ср 


pondiente. 


calcula el mejor calor де C, а partir del mejor valor 


ёс 
o 


N o 6 obtenido y dol valor de C/N calculado on el Paso de 


Como puede observarse el cambio primordial en este método ces 
pecto al método VTM consiste єп que limita el análisis рог minimos cua 
drados a la determinación do N о G, Atp y Pap. La constante C se toma 
lucgo proporcional a N o G ya que C/N o C/G se calcula independiente - 
mente. Pov tanto, la exactitud de los resultados, dependerá del сопо- 
de 


cimiento esprto de los factores volumétricos iniciales, B о B 


о: ві? 


la compresibilidad a № -vna y de la saturación de agua сп la zona do p2 


tróleo. Además, en cl caso u. . ^imiento radiales, entra сд +, ‘ом de 
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6, o .ngulo de la civeunferencia sobre la cual actúa el acuífero, 


укедо __ da _ Stanley"? . - Este método, como los anteriores, permi- 


to calcular el petróleo o gas originalmente en el yacimiento y la cons 
tanta da intrusión, a través de un proceso de ensayo y error (tanteo). 
Utiliza las mismas ecuaciones de VTM, pero en este caso el término de 
se expresa por la ecuación, 


intrusión, Wa, 


Boj (n+1-3)” (166) 


ado pov (u:1-3)", 


EL valor Q[(ta-tj-12)5] de la Ec.(122) se ha recompila 
donde j es un número entero que varía спите 1 y n (níimero de periodos 
considerados en la historia del yacimiento); m es un exponente que va 


ría eutee 0,5 y 0,8. TL valor de C, es igual a la do la Ec.(122), mul 


m 


tipiicada por la erestente (Atp At) donde Atp es la constante de con 
vovsión de tiempo timensional a tiempo adimensional y At es ùl valor 
del iatervalo on que c2 ha dividido el tiempo де la histoc^i del yaci- 
miento, El mótodo aauue que el tiempo ве divide en periodos de igual 


ducación. 


Como se discutió anteriormente (tj) es función de ty. Stan 
, р D 


ley consideró, 


00650 = tp)” (167) 


Pero, tp = (Atg)(t) y a su vez, := (i)(at), luego, 


су) = (atp 4t)" (D) (1674) 


Q(tp) > (constante) Ci)" (167b) 
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donde 1 = número de períodos. 


junto con el valor de С de la Ec.(122 


lor de С de la Ec(166) debe tenerse 


La historia del yacimiento 


ración 1, 2, sse j.. n. Al aplicar 


va el cálculo де Че, se toma el Ару y 
que sigue en orden descendonte, etc. 
Por 


(037, el Ару con el (n-1)", otc. 


се el término (041-1) en lugar de i d 


En forma similar a los caso 


t(p,t) = 


Y 


CURE 


Stanley mismo arreglo de las e 


de VTM en línea recta deducido anteri 


cuerdo al tipo de yacimiento considor 


los y saturados y yacimiento de 


У = CX 

Я = CK 

Los valores de Y tienen la misma defi 
(192) y X la misma definición de las 

nícndo en cuenta que pava f(p,t) зе ц 


Пі procedimiento de cálculo 


pasos: 


La cosntante 


(St, At)? se ha incluído 


р 
). 


Por tanto, al determinar el va 


en cuanta esta constante, 


se divide en períodos de igual du 
el principio de superposición pa 


el Q(tp) máximo, el Ару y el Ер) 


Por tanto, el Ард se tomará con 
esta razón, en la Ec.(166) apare 
e la Ec.(167b). 
s anteriores, 

(168) 
tp,t) (134) 


enaciones ĝe belance Ја materiales 


ormente, Ess.(130) y (14%) y de a 
ado; yacimientos de petróleo no - 
gas, 

н (136) 
+ 6 (145) 
nición de las Ecs.(133), (139) y 
Ecs.(130), (180) y (193) pero te~ 
sa la nueva definición, Ec.(1££). 


puede resumirse en Jos siguientes 
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(1)-Se divide la historia del yacimiento en un nümero, n, ade 
acuado de iguales incrementos de tiempo y se numeran ascendente, 1, 2, 


3, ... j “ce Ne 


(2)-Se calculan las caídas de presión de acuerdo а la 


(125). 


(3)-Se calcula f(p,t) de acuerdo a la Ec.(188) para m = 0,5 


(h)-Se Calcula Y y X pira сада poríodo. 


(5) Se asunen valores de N o б. Para cada N o G asumido se 


calenlan una serios do valores de C (vno par cada poríodo), tanto pa 


ra m = 0,5 como m = 0,8. 


(0)-38 gráfica C como función de tiempos una vneva pava cada 
N о G asumido y para cada м = 0,5 y 0,8. De estos gráficos puede dedu 
cirse si el valoc do N o G verdadoro está couprondido зайге. los valor 
ros asunidos, o^;ervando si ocurre cambio de inflexión de las curvas, 
es decir, caubio de curvatura hacia arriba a curvatura hacia abajo, a 


medida que N о G aumentan, En cazo de que se tenga previamente un es- 


timado de М o Су además se conozca que el acuífero va in 


paso puede climinarse, pue canos seguros que en el intervalo de v 


lores asumidos de N о С está correcto. 


(7)-Рага cada N o G acunido se calcula el promudio arimótico 


de С, denominado Cp» у se са}сч1а la desviación normal corroaupondiernte 


E HANG. _с.)2 ; 
g = | 2 n (6, 21 (102) 
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Se gráfica 0 como función de N o G tanto para m = 0,5 como m = 0,8 y 


donde ocurra el mínimo de desviación normal alli se toma cl mejor valor 
de Н o б y m. El mejor valor de C, será el Cy calculado para los mejo- 


ros valores de N o G y m obtenidos, 


Debe tenerse presente за el Paso 5) que los valores de C cal 
culados al comienzo de la historia del yacimiento pueden ser muy dife- 
rentes de los calculados para el resto de la historia. Este comporta- 
miento puede duberse a que al comienzo de la historia, la intrusión uo 


siga exctamente la Ec.(166) o porque la presión по sca exactamente la 


que ocurre en el contacto agua-petróleo, etc. En tales casos зе го 


mienda no tomar estos valores de С para el cálculo de Ср. 


Una voz obt:nidos los mejores valores de Мо G y за pue- 


deu изар en la predicción en la forma conocida. 


El procedimiento anterior se aplica а acuíferos infinitos. 
Рига acuíferos finitos se usa la misma técnica pero al graficar ea el 
Paso 6), зо observará, que cuando el límite del acuífero llega а ser 
efectivo, C conenzará а disminuir a una rata constante, no importa el 
valor 40 Мо С asumido. Ги este caso, el valor де N o G зе cbteadri 
del Paso 5), analizando las diferentes curvas y seleccionando а de me 
nor variación de С con tiempo. Además, para la predicción es necesa- 
rio obtener una ecuación de С con tiempo. Ln vista de lo anterior, el 
mátodo es más recomendado, cuando se trate de acuíferos infinitos, 


Tiene la ventaja de ser wn mátodo sencillo y puede trabajarse fácil - 


mente con caleuladora de escritorio. 


Los métodos de análisis de la historia del yacimiento discu 


а = 


tidos en las secciones anteriores, son métodos relativamente simplifi- 
cados, purs no toman en consideración la reducción (disminución) de la 
zona de petróleo, debido 2 la intrusión (invasión) de agua y а ја ex- 
pansión de la capa de gas, en caso de que éste exista. Un análisis más 


completo debe considerar catos factores. 


Eficiencia de desplazamiento es una medida de la cantidad de 


petróleo que el agua да intrusión desplaza de la zona de petrólen y en 


caso de existir capa de gas inicial y cuyo comportamiento sea de expan 


sión, es vna uadida de la cantidad до petróleo despla 


Ло al expandir- 


se tal capa de gas dontvo do la zona de petróleo. 7 


definen, 


1) Ificiencias ds desplazenisnto on la zona da de 
ума, 
Intvusión neta de agua 
EDig = 7 “ТТ A иены (170) 
* Volumen de potvól:o de Ja zona de реьгбїео 
invadida por agua 


2) Eficioacia de desplazaniento en la zona de expansión de 
la сара de gas, 


Expansión total de la сара de Zas 
EDia AAA م س ت ےم‎ A O С س‎ ея (171) 
а Volumen de vetrólso де la zona de petróleo 
invadida por раз 


La intrusión ueta de agua es simplemente la intrusión total 
menos el agua producida, CAEN ni tiomino Че debe determinarse con 
tiempo según 32 ha discutido y Wp es dato de producción. һа expansión 
total де la сара бе gas astá dada par (Bi - бре) ба - 6, о sen que de 


be conocerse el gas originalmente en la copa de ваз y n1 gas producidd 


Gpe? de esta capa. Ambos términos asumen que el petrólco que queda el 
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las zonas invadidas uo sufren cambios (merma), ni ocurre liberación de 
gas de tal petróleo. Si esto último ocurre, tal gas liberado debo au- 
mentarse a la intrusión neta o a la expansión de la capa de gas, y los 
términos anteriores habría que modificarlos, dividdiéndolos por el fac 
tor [1 - SopíRsi-Rg) Ва/Во(1-Зог)] › donde бор es la saturación residual 
de petróleo en tales zonas de intrusión o expansión. Además se consi- 
dera que todo el putróleo ha quedado atrapado en estas zonas a la pre~ 


sión inicial, de lo contrario es necesario modificar el valor de Ryi. 


Los volúmenes de peiróleo de las zonas de poirólco invadidas 
por agua y gas, pueden a menudo obtenerse del comportamiento del cumpo, 
si зе logra conocer ol volumen del yacimiento invadido рог agua о gas 


sa 


de la сара de gas. Conociendo 13 saturación inicial de petrólen, se d 


3 


termina el volumen de petvódeo en tales zonas, 


Los volumones de yacimiento (formación productora) invadidos 
por agua o gas, pueden determinarse por el movimiento de los contactos 
agua-peteóleo y pas-petróleo. La posiciones de estos coutactos pueden 
estimarse de pozos o por registros. Si los contactos no se mueven uni 
fornente, ol voiumean invadido debe determinacse en cada poze y luego 
por mapas isopacos, determinar los volumenes invadidos. Las eficion - 
cias de desplazamionto de la historia, pueden usarse para predicciones, 


siempre y cuando las condiciones futuras de operación permanezcan cela 


tivamente iguales. 


i las eficiencias del comportamiento no pneden determinarse 


e 


dcl comportamiento del campo, existen otras técnicas quo pueden usarse, 

entro las cuales la más común es la 46 welge 5, El método de análisis 
OTT aM PENAT 2 raaf 

y prèdicciön de Woody y Moscrip"? también incluye el cálculo де eficien- 


cias de desplazamiento. 
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Weigo 26 basa su mátodo en las ecuaciones de desplazamiento 
de Buckley y Leverett"?, para lo cual,entre otras cosas, 39 requiere co- 
nocer las curvas de permeabilidad relativa. La ecuación de flujo fraccio 
nal, para una fase desplazante d y petróleo, despreciando el término de 


capilaridad, puedo escribirse, 


(0,8883(1077) А К КроАр sena 


Ра = سسس‎ 1 ат (172) 
хора Но At 
1 + س‎ 
Kd Yo 


dorde d se refiere a la fase desplazante que puode aer agua о Баз. El 
término Ap es la diferencia cn densidades: para petróleo y agus, áp = 


(ру-од) y para petróleo y gas, do = (P07Pg)» La Pig.59 presenta una 


curva tipica de fq como función de Sd (5, о Sg) y la curva de distri- 


bución de saturación correspondiente, El valor de бај corresponde en 


el caso da intrusión de agua, a la saturación de agua existente en la 


¿ona de petróleo antes de сопеплаг la intrusión y en el caso de expa 
sión de la capa de gas, a la saturación de gas existente en la zona de 


pesvóleo antes de comenzar la expansión. Este valor es genecalmente ce 


ro o muy cerca de cero y en tal caso la curva de flujo fraccional parte 


de cero. La tangente а la curva de fas a partir de ба}, determina la sa 


за (invasión) o expansión y la se 


turación бад en од frente de intrus 
turación promedia en dichas zonas, Sap> tal como indica la Fig.,59. De 
estos datos, puede determinarse la eficiencia de desplazamiento, Аде ~ 


más se determina el valor de saturación máxima de la fase desplazante, 


Sámaxo У POF consiguiento la saturación irreducible de petróleo, 35, 


El tórmino que incluye Ар, denominado término de gravodad, 


muestra oteos factores, además de la razón de permsabilidades relati- 


"253A ёр орош? 203 OJUANOZO]A 
041 cinj} эр $206 


-ѕәр әр 50:200:21:18 10018439 


127420 Я ШЕ 
6219499 - (506 ©) cejpsjad - (228 о) 
0050 6:201:02 сабо 06122100 
AA 
1 0:3421510 pj ‘24002010590 0 
: о 


волі) — 55 f! 


012099:01 
өр спа 


ЖЭ 


NOI2YuUn1VS 36 №01999 84 00074 30 VASNO 


174834) мО!20131:15:0 33 VASNI 
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$ 
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уаз, Ко/Ка y de la razón де viscosidades, нфд/Мо» que influyen en la cur 
va de fq contra $a y por consiguiente en la eficiencia de desplazamion 
то. Tales factores son: А, К, Но» APOG y Ч. Cada uno de estos facto- 
res, puede, teoricamente, estudiarse por separado y verse como influ- 
yen en la eficiencia de desplazamiento. Por ejemplo, un aumento en la 
rata de producción, aumentará las fracciones de flujo de agua o gas ра 
ra una saturación dada. Esto, por supuesto, disminuye la eficiencia 


де desplazmiento, 4 


Sin embargo, уа se discutió respecto а los, fucores que deben 
tenerse en cuenta en la rata 12 producción para una recuperación máximas. 
Entre otros, la permeabilidad da la fovmación, К, y la viscosidad del 


petróleo, son factores muy influyentes. 


Woody y Houscpip'?, combinan Лаз ecuaciones de balance d^ ma- 
terialos, de intrusión de agua y de desplazamiento (Buckley y Leverett) 
y las roSuolven simultáneamente. El petróleo originalmente en el yaci- 
miento lo consideran conocido. Determinan a diferentes tiempos de la 
historia del yacimiento, la constante de iutrusión de agua y las posi - 


ciones de los contactos agua-petróleo у gas-petróleo, lo mismo que las 


distribuciones de saturación en las zonas de intrusión y expansión. Con 


estos datos, obtienen las oficiencias de desplazamiento а estos 
D 
tiompos. Además cl mótodo permite estudiar la variación de la rata de 


producción con el fin de detevainapr la rata de producción para máxima е 


сісіспсја. Los cálculos Бол bastante laborioses e incluyen un proceso 


de ensayo y error (tanteo), por tanto, deben efecutarse а través de una 


computadora. 


El абтодо de Woody y Moserip también puede usarse para prede- 


cir el yacimiento. En este sentido, los cálculo exigen conocer la pre- 
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sión y la rata de producción para cada período futuro en se desee la 
predicción. Esto introduce en ol método una serie de cálculos de ensa 


yo y error: se asuma la presión y la rata de producción,se hacen los cal 
culos y se compara con los valores estimados. Si los dos no verifican 
satisfactoriamente, deben revisarse los estimados y repetir los cálcu- 


los hasta obtener una solución satisfactoria. Así se continua para ca 


da tiempo futuro. 


1 2 "m А 
Rachford y otros"9, presentan un método numérico que permi- 


to también calcular el movimiento de los contactos agua-petróleo у 635 


-petróleo con ticmpo а partir de la historia de presión y producgión. 
El método tanbién pernite predecir, asumiendo valores de las varones 


паз-реїгд1ео y agua-petvóleo. 


рыш: PORTANTENTO DU 


ACINILNT МРОТЕ ÜIDR 


JON "BL CO 


El tármino, predicción de yacimientos con expuje hidráulico 
puedo tener varios significados de lo que ue desee predecir: presión, 
rata de intrusión, intrusión total. etc. Por tanto, es más correcto, 
on forma similar al caso de yacimientos con empuje por gas on solución, 
hablar de prodicción del comportamiento de yacimientos con empuje hi~ 


dráulico, lo que combina predicción de pres ón,intrusión de agua, varón 


gas-putróleo y agua-petróleo, etc. 


El aátodo y la exactitud de una prodicción del comportan 


to de un yacimiento con empuje hidráulico, depende dci método con que 


se haya analizado la historia del yacimiento y de las suposicicnes he- 


chas para la predicción. Гог ejeuplo, las mismas ecuaciones que Wo 


оду у Moscrip usan para el análisis de la historia, se usan para la pre 


€ 


dicción tal como се moncionó anteriormente.: 
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En esta sección se presenta un método de predicción del comportamiento 
de yacimientos con empuje hidráulico, en base a que la intrusión de a- 
gua, pueda expresarse por una de ¡as fórmulas vistas anteriormente. bs 
te mótodo de predicción es un тё сойо de iteración (repetición de calcu 
los) o ensayo y error (tanteo) y en este sentido se asemeja a algunos 
de los métodos de predicción de yacimientos de empuje por gas en solu- 
ción. Además, en un empuje hidráulico, el tiempo es Un factor en са- 
da punto de la predicción. El método asume conocida la rata de produc 
ción de petróleo o permite determinarse de la historia del yacimiento. 
Igualmente, asume conocida la producción de agua, expredada generalmen 


te por la razón agua -potróleo, RAP. 


La rata de producción de petróleo puede estimarse a partir 


dal Índice de productivifad, j, 


Qo = 0(ре-р,ү? (173) 


donde J puede ser el Índice de productividad para cada pozo o del cam- 


po en general. e último caso, Оо será la rata de producción to- 
tal dol campo y en el primer caso, la rata total de produccióu, será ја 


suma de las ratas de producción individuales de los pozos. Además on 


esto caso, ol valor de pa es aprozimadamente igual à la presión del ya 


cimiento, valoe asenido en la pred El valor de рур es la pre- 


sión fluyonte de cada pozo, EL valor de J puede expresarse por, 


0,90708 k kpo В 
FE E EE (174) 


но Во lu(eQ/rg) 


El gvupo de tá evaluarse de la prueba inicial 


e 


kh/la(r,/r pue 
kn/la(r,/r,) pues 


de productividad del pozo. En eae momento Кро ёл muy coreano а la uni 
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dad y ug Y Во se evaluan а la presión inicial. En esta forma podría ob 
tenerse en índice de produtividad inicial, Jj. El índice de producti- 
vidad permanece prácticamente constante para presiones de producción 
por encima del punto de burbujeo. А medida que la presión disminuye se 
forma saturación de gas en la zona de petróleo, Кро У Во disminuyen у 
ug aumenta. Muchas veces, el cambio principal en el J puede ocurrir 48 
bido a la reducción de ha medida que el agua, gas,o ambos invaden el 


intervalo productor. 


La predicción de ratas de producción de pozos parcialmente in 
vadidos por agua es difícil, ya que Ку no se conoce. En la práctica, 
se produce muy poco petróleo una vez que ocurra la ruptura, Además es 
recomendado reparar los pozos раса prevenier producción excesiva de а- 
gua. А veces para predecir las ratas de producción de petróleo en pozos 


de alta producción de gua, зе hace uso de Jas cuvvas de dezlinación. 


Los siguientes pasos resumen el procedimientoa de predicción 


para el caso dc un empuje hidrfalico donde no existe сара de gas ini - 


ones порт 


spon- 


cial. Ea caso que exista, bastará modificar las ecuac 


dientes, por el facotr m. Adumás, tampoco existe ningún tipo de inyec- 


ción: 


(1) - Se selecciona un incremento (intervalo) de tiempo АЕ, 
generalmente un mes, un trimestre, un semestre o un año. En esta for 
ma se determina од бістро а la que se van à efeciuer los calculos, 

суо; tuo. + АЕ (17%) 


donde j es el nümovo de la etapa. 


(2) - Sa fija una rata de producción promedia del yacimiento, 
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Одо, para el período considerado, de donde puede escribirse, 
ANS = 004 (175) 
ЦЭ = ШЭРН + ANp (177) 


(3) - Se asume una presión Pj. Esta será la variable de en- 
sayo y error (tanteo) en los cálculos para la predicción. 

(+) - Usando el tiempo ty y la presión asunida p4, se deter- 
mina la intrusión de agua de acuerdo al tipo de intrusión considerado, 


Ses. (205), (112) о (122), 
мер = | вэ.) (178) 


(5) Se calcula la saturación ёв р 


¿tróleo en la zona de petró 


leo por la formula, 


“Np j - 


Soj E (179) 


“Boi 
(1 + 331) Sus - Swi) 


Esta ccuación es la rümoro (93) de la roferencia 1. Asume quo no exis- 


te capa de доо inicial, o si existe, ocurre difusión 021 gas en el pe- 


tpróleo y el contacto aguarnernóloo по so uae. En caso de existir ca- 


isma, se usa la te.(101) d2 la Refacen- 


cia 1. Еп la cencción anterior, ну бүр “on las saturaciones de ас 
gua (incluye Syj) y saturación de gas за la sona invadida por el aguas 


o ssa saturación ¿ccolucible de petró- 


además, (1 - Syp ~ бре) = 


leo en la sona invadida por agua. Estos vcloros do saturación daben 
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conocerea, bien sea de pruebas de laboratorio, de campo, método Helge, 


сте. 


El valor de Ур} debe determinarse del comportamiento pasado 
(historia) y generalmente s2 expresa como la razón agua-petróleo, RAP, 
с sea una fracción de ANy. El valor de Wej es el calculado en el Paso 
(4). Los demás datos son función de presión (Boi; Boj? о constantes 


(54193. 


(6) - Con el valor de Soj Y las curvas de permeabilidad rola- 
tiva en la zona de petróleo se determina la razón de permeabilidados rela 
tivas, (kg/ko)3» de donde se determina la vazón gas-petvólco instanta- 

А 


nea а la presión pj» 


Мо Во 
во в (| سل‎ прав ا‎ (189) 
1 jl ng jli 

(7) - Se calcula el valor de là галоп gas-potróleo acumulada, 


Rp, а ла presión Pj» 


(181) 


donde, Ар) = Npj - ЦЭНЭ 


(8) - ба calcula la intr ёл de agua, Меј» usando ia 202: 


ción de balance de materiales, 


Е (Brg Pri) + Ур] (182) 


(9) - Se compara el valor de Wa del Puso (8) con el valor de 
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Wa del Газо (+). Si coinciden dentro de una aproximación normal, la 
presión asumida р} es la correcta y se continua la predicción para el 
tiempo jri. Si son muy diferentes es necesario asumir una nueva pre- 
sión y repetir los cálculos hasta obtener valores satisfactorios. 0b- 
sérvese que Npj y tj permanecen constantes y solamente varía la pre- 


sión. Un aumento en presión,aunenta la intrusión total de agua por la 


‚а 
Fórmula, W¿=Cf(p,t) у la disminuye al calculavse рог la ecuación de ba 


lance de materiales. Por tanto: 


a) Si Ч, (Paso 8) > Hg (Paso !), se disminuye la presión 


b) Si Не (Paso 2) < Му (paso №), se aumenta la presión 


La predicción hasta este paso, permitiria obtener p, Че» Ny, como fun- 


ción de tiompo. ha predicción puede hacerse más completa, si сп el a- 


de la historia so ha determinado el avance de agua-potróleo den 
tro de la zona de petróleo y la eficiencia de desplazamiento, Además, 
también puede verificarse si la rata de producción usada en сада in- 
tepvalo es la correcta, sieupre y cuando de ја historia se haya detar- 
minado una relación entre producción y volumen de la zona de pereóleo 


Tavadido. En csta forma puede agregarse los siguientes pasos en la 


predicción: 


(10) - Se determina el volumen de petróleo de la zona de pe- 


tróleo invadido por agua о gas, conociendo la eficiencia de despla 
miento en estas zonas y la iatcusión neta de agua y expansión total de 


la capa de gas, Ecs.(170) y (171). De este valor зе deteeminan el vo- 


lumen poroso invadido y pueden localizarse los contactos agua-pot 


y gas-potróleo. 


(11) - Si se ha obtenido una relación entre la producción y el 


volumen de la formación invadida se verifica si el volumen invadido 
calculado en el Paso (10), está dentro de una aproximación normal 
(t 10 por ciento, puede considerarse satisfactorio). Si ocurre una 
diferencia muy grande, se toma otra rata de producción y se repiten 


los cálculos hasta obtener verificación aceptable. 


También podría considerarse el caso de un mantenimiento de la 
vresión en un valor p. En este caso, odría estimarse la cantidad de 
р 


fluído que es necesavio inyectar en la siguiente forma: 


(1) - Igual ul Paso 1) del procedimiento anterior. 

(2) - Igual al Paso 2) del procedimiento anterior. 

(3) - Se calcula Ус = C f(p,t) y luego el incremento en intru 
sión АЧе = (Wej - Woj-1- : 


Че por balance de miteriales y el incremanto 


de inter 


la inyección nota al yacimiento Читайте el in 


фогуајо considerado, 
Inyección neta = AWe (Paso 5) - Ma (Pazo 3) 


га el tiempo desca- 


En esta forma se repiten los válculos ha 
do, de acuerdo al agotamiento, invasión, etc. Ређа tenerse on cuenta 
que si la intrusión sigue flujo no-continuo y el acuífero es limitado, 
el valor de АМе del Paso 3) llega un momento en que es coro, de acuer 


do al valor de tp. 


PREDICCIÓN EN PASEA. 


ВАЗЕ $0 DE RATA TERMINAL CONSTANTE: 


SOLUCT 


ON. BET, 


Las soluciones de la ecuación de difusividad para el caso de 


rata terminal constante fueron dadas anteriormente, El teorema de su 
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perposición permite usar tates soluciones cuando la vata de intrusión 
cambia con tiompo en la forma indicada en la Fig.55. En forma general 


punde escribirse, Ecs.(79) y (89), 
on 
Ар = рүтрд FA E Сезе 4p [taj 125] (183) 


donde А es una constante que depende del tipo de acuífero: 


Acuífecros Lineales, 


887,50yL 


А = س‎ (188) 
kbh 
Acuiíferos Radiales, 
191,2 y 
А = o———— 105) 
kh 
Además, eg = 0 y to = 0. Para acuiferos lineales, Ар = Ру. Esta воћа 


ción рис4е cup 'varse ва dos formas para la predicción: 


1.- Pvedicción de presiones en el contacto agua-patróleo, pa- 


га una reta de intrusión dada. 


2.- Cálculo de ratas de iatrusión, para una presión dada en 


al contacto agua-petróleo. 


Eu cada uno de estos casos, de la historia (o predicción an- 
terior) se conoce La presión como función de rata de intrusión, hasta 


el tienpo tp, Гс.(183). Esta ecueción, para nn tiempo n+l, se escribe, 


n 
Pi Pnl = А T meu Ap lena tj-13, + (опат ба) Ар Сава tao! а 
j= 
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Esta ecuación puede usarse para predecir los dos casos mecionados. Si 
se conoce n+l al tiempo tg41, se despeja рд+і У si se conoce este va- 
lor se despeja ёрүүе En el siguiente paso, Pnt Y £n+1 serán conoci- 


dos, y pasan а ser рр y ёп. Se hace de nuevo la predicción para el 


tn+1» Que en realidad será tp+2. Se procede así sucesivamente. Los de 


са 
3 


más términos de la ecuación son conocidos. Ambos casos suponen conoci 


da la constante А. 


Si se conoce el valor del petróleo o gas orginalmente en el 
yacimiento, de la ecuación de balance de materiales puede calcularse 
Wo a diferentes tiempos y de aquí cálcular los valores de e. Luego, 
aplicando 12 ecuación anterior se calculan los valores ¿e p y se com- 
paran con los de la historia. Puede determinarse un valor de А que 


расе un chequeo aproximado de presiones. De la variación de 0 con 


tiewpo, puedo extrapolarse e y predocir р de la ecuación anterior, А 
veros, y particularmente en yacimientos que producen рог encima del 
punto de burbujeo, la rata de producción de petróleo puede consideran 
se igual а la rata de intrusión. fn este caso aún es más fácil heeor 


Ја predicción. 


El mejor valov de А puede obteners: fácilemente de la histo- 


ria del yacimiento cn la siguiente forma: 


1.- De la historia del yacimiento se calcula Ар = (pi-p3) 
denominados Ар} obse: vados,(Apjdobs ° 

2.- Conocido las ratas de intvusión, se calcula la sumato- 
ria де la £c.(183), una sumatoria para cada tiempo. 

3.- Se calcula un primer valor de A, dividiendo la suma de 
todos los Ap observados por la sumatoria de todes los términos (suma- 


torias) calculados un el Paso 2), 
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V.- Se calculan valores АВ; denominados бр] calculados, 
(Ар) сат multipicando el valor de А obtenido en el Paso 3) por cada 


una de las sumatorias del Paso 2) 


5. - Se obtiene la desviación normal correspindiente, 
n 
1 2 
= — r | - Р 
o p^ do САР сал 7 (883 obs] (187) 
qa 


6.- Se asumen otros valores de Ау se repiten los Pasos 4) y 
5 


5). Finalmente, se gráfica o coutra Ау donde ocurre el mínimo valor 


de о alli че toma el mejor valor de A. 


Además de los métodos estudiados anteriormente para el anflí 
sis y predicción de yacimientos con empuje hidraílico, ot»as técnicos 
nan sido usadas para hacer este tipo de estudios. Se discutivá en ез 
ta sección los métodos de Hutchinson y Sikora"? y de Funciones de Iu- 
о550 а 53 


Румепојај3»3% y brevemente las técnicas de Modelos Analógic 


(analizadores o computadores analógicos eléctricos) у Matemázicos. En 


particular, está última técnica está desplazando casi tutalomente el 
estudio en modelos analógicos, debido al desarrollo de computadores de 


mayor capacidad. 


odo 


Hure 


Este método permite analizar y predecir la historia de un ya 


lico de un modo general, sin necesidad de co 


cimiento coa empuje hidrá 
посер la geometría y tamaño del acuifero. Es un método sencillo de a 


nálisis y permite hacer los cálculos manualmente por medio de una cal 


culadora de озогдтогјо, 
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El método exige el conocimient» uel petróleo originalmente en 


el yacimiento y la historia normal del mismo: presión, datos de produ- 
ссібп, РУТ, etc. Пе la ecuación de balance de materiales, se calcula la 
intrusión de agua y luego la rata de intrusión. Con datos de presión 
y rata de intrusión se desarrolla una ecuación que relaciona estos 
términos a través del factor, denominado Función o Factor de Resisten- 
cia Z, el que incluye la geometría del acuífero, forma y distribución 


de flujo. La variación de este factor puede extrapolarse para prede- 


cir el yacimiento. 


Las Ecs.(26a), (28a) y (29a) expresan la caída total de pre - 
sión en el contacto agua-putróleo para una rata constante de intrusión 
en acuíferos lineales Y las Ees.(57) y (63) expresan igual сова para 


спев puedan : .uribirse, 


acuíferos radiales. En gencral, tales 
Ар = (pi-p) = А e Z(tp) (188) 


donde A es una constante que depende de la geometríe del acuífero, е, 
es la rata de intrusión y Z(tg) es una caída de presión adimensional, 
que puede hallarse en tablas en función de tp, según se discutió ante 
riormente. Si se yráficara Ар como función de t, рага una rata de in- 
trusión constante, e, se obtendrían las formas de las curvas 1 y 2 de 
la Fig.60, para acuíferos infinitos y finitos respectivamente, Al vo- 
mienzo, la presión caerá rápidamente а medida que el agua comienza a 
moverse en la vegión próxima al contusto agua-petróleo. А medida, que 
mayores volumens de agua entran en movimiento, la caida total de pre- 
sión con tiempo disminuye. Si el acuífero es uniforme en propiedades 
e infinito, Curva 1, la caída total de presión siempre aumneta cor Чет 


po pero cada vez а una rata menor а medida que al tiempo avanza. Un 


Presión, 5р = (pi-p) 


س 
CON COMUNICACION‏ 


RO 
IÓN SUPERFICIAL 


da Total de 


Ca 


Tiempo, t 
Fig.60 - Caída total de presión en el contacto agua- 
petróleo como función de tiempo para una ra 


ta de intrusión, e, constante y de acuerdo 
al tipo de acuífero 


gráfico de Ap contra log(t) para este tipo de acuífero presenta una 
pendiente que disminuye constantemente con tiempo hasta que se hace 
constante, Fig.19. Esta propiedad es básica para la predicción, со- 


mo se verá más adelante. 


Si el acuífero es limitado (finito) y cerrado, el área de dis 
turbio alcanzará el límite exterior y de allí en adelante la rata а 
caída 


la cual la presión disminuye es constante, Curva 2, ya que la 


de presión en el acuífero, se debe al agotamiento del volumen limita 
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do de fluído del acuífero. En este caso puede escribirse, 


dp e 
быш کاو‎ (189) 
dt сҮр 


donde с es la compresibilidad del agua y Ур es el volumen poroso del 


acuífero. 


En la Fig.60 también se indican otras curvas: Curva 3, сазо 
de un acuífero que comunica a la superficie y por consiguiente después 
de determinado tiempo, la caída de presión con tiempo es cero. Final 


mente, la Curva +, donde ocurre inhomogencidades en el acuífero. 


Las curvas de esta figura presentan una característica impor 
tante: la pendiente siempre disminuye con tiempo o se hace constante 
pero nunca aumenta, Esta propiedad se expresa diciendo que la pendien 
te disminuye monotdnicamente . Este principio y la forma de las curvas 


es la base del método de Huchinson y Sikora. B 


En la discusión correspondiente el principio de superposición, 
se dedujeron ecuaciones para calcular la caída total de presión en ol 
contacto agua-petróleo, cuando la rata de intrusión varía con tiom- 
ро еп la forma indicada en la Fig.54. Tales dcuaciónes (Есз.(81) y 
(863), pueden escribirse en forma ganeral, 


срна | (ео) EU Entar- o] (130) 


> 
о 
u 
> 
"ms 


donde А es una constante que depende del tipo de acuífero y F(tp) оз 
una presión adimensional que puede ser Fk(tp) o Ap(tg) de acuerdo al 


acuífero, según indican las Ecs.(81) y (86). 
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La presión adimensional (єр) depende del tiempo adimensional 
tp y êste depende del tiempo dimensional, t, a través de la constante 
Atp. Рог tanto F(tp) depende de t y puede escribirse F(t). Multipli- 
cando éste valor por A y definiendo el producto como una nueva función 
de resistencia, Z(t), o simplemente función 2, la ecuación anterior pue 


de escribirse, 


ваза i [2Cn7tn-j) Сатен 1 -40) (191) 


En esta ecuación, la función Z se obtiene a los tiempos dimensionales 


(tnTtn-3) У lta 0+1) 


Si lar intervalos de tiempo, (61,1753) = АХусоп iguales, lue 
go, tg = ПА y tj = jàt y la ecuación anterior puede escribirse, 
n 
Ар = E enp1-3[2(38t) - 2[(5-1)0t)) (192 
151 
п 
Ар = E en+1-j 154 (193 
3=1 
donde, 
42) = Z(jat) - 2(j-1)At (194 


Si se conoce una curva de 2 como función de tienpo, puede calcularse 
la caída total de presión en el contacto agua-petróleo por medio de la 
гс. (193). De la historia del yacimiento puede obtenerse tal curva de 
Z como función de tiempo у 5U extrapolación, de ucuerdo a la forma, pue 


de usarse para la predicción, usando también la Ec.(193). 


Los siguientes pasos, resumen la obtención de la mejor función 
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Z contra tiempo! 


1. - Se divide la historia del yacimiento en, iguales interva- 
los de tiempo, At, y usando la ecuación de balance de materiales, 56 
calcula la intrusión total de agua Wg, а cada tiempo de la historia. 


Por supuesto,esto considera que se conoce el petróleo originalmente en 


el yacimiento, N 


2 - $e calcula la rata promedia de intrusión de agua a cada 
tienpo, dividiendo por At, el promedio de la diferencia entre el $, al 


tiempo posterior menos el Не al tiempo anterior, 


Wej+l Wej-t 


)195 اتس انس = ер‏ 

3 ( 
2At 

Como antes se mencionó, algunos yacimientos que producen por encima del 

punto de burbujeo, la rata de producción puede considevavse igual a la 


pata de intrusión. En este caso, los Pasos 1) y 2) podrían eliminarse. 


3 - De la historia del yacimiento, se calcula la caída de pre 


sión total, 5p; * (рі-ру), denominada caída de presión observada, (4P Jobs 


» ц ~ Se asume una primera aproximación de Z como función de tiem 


po. Si cada periódo se numera ascendentemente, 1, 2, 3, œ J 


donde j, varía de 1 а n; se toma como primera aproximación de 


cuadrada del nûmero del intervalo, 


ай (71,2, 3, 0... п) (196) 


5 - Se determina 47) y рот medio de Ја Ec.(193) se calcula u- 


na caída de presión total, denominada caída de presión calculada. De- 


be tenerse cuidado al usar esta ecuación ya que los t rminos de la su 
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matoria estan arreglados en orden inverso, es decir, la rata de intru- 
sión para el intervalo n se multiplica por el AZ4 para el intervalo 1; 
la rata de intrusión para el intervalo (п-1), estará multiplicada por 
el AZ, del intervalo 2, etc. Por tanto, es conveniente escribir una 
columaa iavirtiendo el orden de la rata de intrusión, y así la multi - 
plicación de esta columna con la columna de los АЕ; dará directamente 


los Ap calculados, (Apj)oa] * 


6 - Se divide cada Ар observado (obtenido en el Paso 3) por 
el correspondiente Ap calculado еп el Paso 5). El cuocionte se denomi 


па, fej» factor de corveción j. 


7 - Se multiplica cada Ai per el correspondiente factor de со 
rrección, fejs del paso anterior: 2; (nuevo) = fej 2; (viejo). бе grå- 
fica Zj como función de j (equivalente 24 como función de tiempo) y se 
traza una curva continua a traves de los puntos, teniendo en cuenta que 


la pendiente siempre debe ser positiva y disminuir o permanecer constan, 


te con tiempo, 


8 - A cada intervalo se lee un nuevo valor де Zj del gráfico 
del paso anterior y se repiten los Pasos del 5) al 8) hasta que los va 
loras de Zj en dos iteraciones sucesivas sean iguales o cambien dentro 
de una aproximación aceptable. Generalmente es suficiente entre tres 
y cinco iteraciones. Este paso tambión puede formularse diciendo que 
las iteraciones se repetin hasta cuando los factores de correción, Fajs 


lleguen a ser relativamente constantes @ iguales o muy cercanos a la u 


nidad. 


9 - Con los últimos valores de Ap calculados y los Ap observa 


dos, se calcula un factor de correción promedio, Ғор, рот la Ёбооч1а, 


n 
E (Ару) оьз 
PEE! 
fop = —————————- (197) 
n 
E real 
151 
10 - Se cálculan los mejores valores de zm multiplicando los 


últimos valores leidos del gráfico en la última iteración (Paso 7), por 


el factor de correción promedio, Ғор» del paso anterior, 


11 - Se determinan los mejores valores de Ap calculados multi 
plicando los últimos valores de Ap obtenidos, por el factor de correc 


ción promedio, f Estos nuevos valores de Ар calculados deben che - 


ор», 
quear aproximadamente con los Ap observados, 

Este último paso no es absolutamente necesario уа que los mejo 
res valores de 7j se obtuvieron en el Paso 10), pero sirve como un che- 


queo más y muchas veces se desea conocer la diferencia exacta en los 


Ap. 


? 


Los mejores valores de 2j se gráfican como fun“ ¿ión de 
j y зе traza una curva continua a través do ellos. La predicción d2- 
pende de la forma en que csta curva se puede extrapolar de accordo а lo 


discutido respecto a la Fig.60. Pueden ocurrir varios casos: 


1) - Si la última parte de la curva de 2; como función de j es 
una línea recta definida, es lógico considerar un acuífero limitado y 
la exploración se hará linealmente, es decir 42; (pendiente) es cons 


tante, en papel normal. 


2) - Si la curva de Zj contra j en papel normal no presenta 


linearidad, puede que la presente al graficarsa como función de log(j). 
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Este caso indica que el acuífero es infinito. Para su extrapolación, 


se consideran dos casos: 


Caso_a. Que el acuífero se siga comportando como infinito 
сп el futuro y зе extrapolará el 2j como ura linea recta al gráficar- 
se como función de log(j). En este caso la extrapolación de la pen - 


diente puede hacerse usando la ecuación, 


j+1 
Тод 
i 
А = А72 ———— — (198) 
1 
log 
j-1 


donde j es un período cualesquiera en el futuro, ј> п. 


Ca 


o b. Que el acuífero sea finito, pero que se haya com- 
portado como infinito hasta el tiempo n. De aquí en adelante 50  com- 
porta como finito y se extrapolará Zj contra j como una linea recta en 


papel nornal, tomando la última pendiente conocida. 


3 - Puede ocurrir que tanto Zj como función de 1, 0 Zj como 


función de log(j) no presentan linearidad. En este caso nada puede de 
cirse respecto а qus el acuífero soa infinito, y lo más recomendado es 
usar el promedio de la pendiente para los 3 o Y últimos intervalos y 
extrapolar como una linea recta en papel normal, es decir, asumir un 


acuífevo limitado. 


4 ~ Ya que la intrusión de agua por caida de presión unitaria 
es menor para un acuífero finito que para uno infinito de caracterís 
ticas iguales, en algunos casos se toma un promedio de los valores de 


7, considerando acuíferos infinitos y finitos, 
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Una vez definida la extrapolación de la curva Función 4e Re- 
sistencia Z, por medio de la Ec.(193) puede predecirse la presión pa- 
ra una rata de intrusión dada, o lo contrario, puede predecirse rata 
de intrusión para una presión dada. Esto se hace paso por paso para 
cada intervalo en el futuro. Por ejemplo, la Ec.(193) escrita para 
el intervalo n+1, a la que se conoce 71,1» será, 

n+1 
(pi - Рафа) = €n+1 ۵71 + 152 *n«2-j ۵2j (199 
De esta ecuación, puede calcularse рд,» Si se conoce еду}, о lo con- 


trario, puede calcularse eg,1 sí se conoce руф: 


Punciones do Influencia. 


El método de Hutchinson y Sikora discutido en la sección ante 
rior es un caso particular de Funciones de Influencia, sólo que tal ca 
so, se denominó Función de Resistencia Z(t). En la presente discusión 
se denomina F(t) con el fin de usar la misma nomenclatura empleada еп 


la iitcvatura. 


En valor de presión adimensional Ap(tp) de la Ec.(83) dependa 
de tp y éste de t, por tanto puede escribirse como un valor Ap(t), 
Este valor multiplicado por A se define como Función de Influencia F(t) 


o simplemente Е y puede escribirse, 
Ар = (pi-pp) = E (епе1-3 - en-j Fj (200) 
Nótese que en esta ecuación se ha invertido la sumatoria comenzando 


por el último caabio en la rata de intrusión que actúa por el tiempo 


menor, luego el penúltimo cambio en la rata de intrusión que actúa 
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рог un tiempo un poco mayor y asi sucesivamento. En esta forma las f 


ciones de influencia aparecen en la sumatoria en orden ascendente. 


La ecuación anterior también puede expresarse en la forma de 


la Ec.(193), 


n 
Ap (ра-ра) = Е ёда1-4 АР (201! 
T 


donde АР; = Fj-Fj 4 


La forma de las curvas de F(t) contra tiempo es del tipo de las curvas 
de la Fig.60. Además, F(t) en unidades prácticas, genoralmente se ех- 
presa en l1pc/bis/día, tal como puede verse de las Ecs.(200) y (92017. 
También de estas ecuaciones puede verse que la presión y rata de intru 


sión están relacionadas a través de la función de influencia, la que 


puede definirse como la respuesta on presión en el contacto a 


leo a una rata de intrusión unitaria causada por una caída de presión 
unitaria en el mismo contacto agua-petráleo. La función de influencia 
refleja la heterogeneidad y geometría del acuífero. Por tanto, existe 


una función particular para cada sistema yacimiento-acuífero. 


Si se conoce la F(t) contra tiempo, pucde predecirse la pres 
sión a partir de un programa de rata de intrusión o viceversa, prede -. 
cir rata de intrusión, 51 se conoce la presión. En este caso, a tra- 
vés de balance de materiales puede ргедесјезе producción acumulada y 


rata de producción de petróleo. 


Como antas se mencionó, el método de Hutchinson у Sikora ез 
un caso particular de obtención de la función de influencia а partir de 


la historia del yacimiento. Las soluciones de Habor?, van Everdingen 
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y Hurst? у Hueller?1 para los sistemas lineal y radial ideales vistos 
anteriormente, permite obtener fácilmente las funciones de influencia 
para cada caso particular, conociondo las constantes А y АХр, en la 


forma indicada al discutir la Ec.(200). 


Hicks y otros3?, discuten еј uso de analizadores elécticos 
para determinar la función de influencia que describe la historia de 
un yacimiento y permite efectuar la predicción del mismo. Coats y o- 
tros? desarrollan un método para obtener la función de influencia 
por programación lineal a través de computadores y a partir de la his 
toria de presión y rata de intrusión. El programa, además de utili- 
zar la Ec.(200), debe cumplir otras condiciones, ya que la sola progra 
mación de la Ec.(200) presenta variaciones en la pendiente de Fj y tal 


como se discutió en el método de Hutchinson y Sikora, la pendiente de 


la curva de P(t) contra t debe disminuir monotónicamente. 


Modelos 


Más comunmente se conocen con el nombre de analizadores elec- 
tricos9* o simplemante mallas RC, donde R significa Resistencia y С 
significa Condensador. Cada modelo se compone de una red de condensa- 
dores y resistencias en base a una división del ávea del yacimiento y 
del acuífero. Cada una de estas áreas encierpa un volumen de forma - 
ción dado (bloque) у la capacitancia del fluido, dafinida como el pro 
ducto del volumen de fluido dentro de los poros y 12 compresibilidai 
del fluido, se representan por un condensador, Los valores relativos 


de la capacidad de los condensadores se fijan por la ecuación, 


М анус 


(202) 
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donde, Ср = capacidad eléctrica del condesador, en mircofaradios; M, = 
constante del escalamiento entre el fluido y las unidades eléctircas, 
en microfaradios/bls./lpc; a = área de cada división de la red, en pies 
cuadrados (ya que el área se divide generalmente en cuadrados, a veces 
se usa 1? en lugar de a, donde L es el lado del cuadrado); h = espesor 


en el área a, en ples; ф 5 porosidad promedia del volumen ah, fracción; 


e = compresibilidad del fluído en los poros, 1рс-1, 


La resistencia al flujo entre los centros de los bloques de la 
red se representa por una resistencia variable, conectada a los conden 
sadores correspondientes. Los valores relativos de las resistencias se 


vepresentan por la ecuación, 


No y 
R 
Rg = ни (203) 
e 0,001127 kh 
donde Вр = resistencia eléctrica, em mega-ohmios; Na = constante de es 


calamiento del analizador, en mega-ohmios/lpc/bls./dia; k * permeabili 
dad, en md. ; и viscosidad del fluido, en cp.; h = espesor promedio del 


bloque, on pies. 


El voltaje en los condensadores representa la presión del ya- 


cimiento, 
AE = РАр (298) 
donde, AE = voltaje a través de los condensadores, en voltios; P = 


constante de escalamiento, еп voltios; Ap = caída de presión entre 


bloques, еп Тре. В 


Finalmente, el flujo de corriente de cada condensador repre 


senta la intrusión de agua o producción de petróleo, 


1=0 (205) 


donde, 1 = flujo de corriente de cada condensador, en micro-amperios; 
: : : " 
Q = constante de escalamiento, en micro-amperios/bls./dia; e = pata 


de intrusión (o producción) en bls/dia. 


La Fig.61 presenta un ejemplo del Campo East Texas (Estados - 
Unidos) rodeado por un acuífero, la división eu áreas (bloques) y el a 
nalizador eléctrico correspondiente, con el que se hizo un estudio del 
campo y predicción del mismoS1, Cada punto representa el centro de un 
bloque donde se asume que ocurra la expansión del fluido. Además, а са 
da punto va unido un condensador. Los bloques están conectados entre 
51 por resistencias. El número de bloques depende de la aproximación 


deseada ya que se supone una presión promedia para cada bloque. 


Para similar la producción del campo East Texas, se retira co 
rriente de la malla (analizador) en el punto situado a la derecha, que 
representa el campo. Por otro lado, para simular inyección de agua a 
lo largo del extremo (límite) oeste del campo, se introduce  corvointe 
proporcional al volumen de agua inyectada en el punto donde se unen 


las resistencias del acuífero y del campo, tal como indica la Fig.61. 


En general, para el estudio del comportamiento de un yacim 
to, los condensadores se cargan inicialmente al mismo voltaje, lo que 
representa igual presión inicial, y sus capacidades son proporcionales 
al producto del volumen comprimido en cada bloque y a la compeesibili- 
dad efectiva. Si cl bloque está localizado en un acuífero, el fluído 
será agua y la compresibilidad es la del agua. Las resistencias repre 
sentan la resistencia al flujo entre las unidades (áreas en que se ha 


diviáido el sistema), En esta analogía, la caída en voltaje es propor- 
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cional a la declinación de presión. 


Con el fin de analizar la historia de un yacimiento y prede- 
cir su comportamiento, las resistencias y condensadores del acuífero se 
regulan por ensayo y error (tanteo) hasta que el instrumento reproduz- 
ca aproximadamente la historia pasada de presión y producción del уа 
cimiento. Para regular (calibrar) inicialmente el instrumento en for- 
ma adecuada, se requiere por lo menos dos o tres años de historia del 
yacimiento. En la práctica, se escoge una escala adecuada de tiempo, 
por ejemplo, un minuto = un año y se varía paso a paso la resistencia 


R, lo que representa el flujo del yacimiento a intervalos exactos. Las 


corrientes producidas, leídas en el amperímetro А, son proporcionales 
а las ratas de intrusión de agua durante cl intervalo. Estas ratas se 
calculan рог los métodos vistos anteriormente, Cuando la curva entre 
voltaje y tiempo, generalmente reproducida en un voltímetro con regis- 
trador, coincide aproximadamente con la curva tiempo-presión de la his 
toria pasada del yacimiento, se asume que los condensadores y resisien 
cias del instrumento están debidamente regulados para representar el 
acuífero. Con el instrumento regulado en esta forma, puede simularse 


el comportamiento del campo cuantas veces зе quiera, bajo las ratas de 


producción e inyección deseadas. 


Los analizadores eléctricos más conocidos son el Analizador 
Carter? y el Analizador Sun, Este último realmente es un Computa ~ 
dor Analógico, pues además de componerse de una serie de resistencias y 
condensadores tiene incorporado un sistema que permite resolver automa 
tica y simultáncamente la ecuación de balance de materiales y la ecua- 


ción de intrusión de agua de flujo no-continuo. Са esta forma se 


elimina la solución por ensayo y error (tanteo) de estas ecuaciones. 
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Tanto los analizadores eléctricos como los computadores analó 
gicos tienen ya muy poco a casi ningún uso en la industria petrolera y 
en su gran mayoría hoy hacen uso de modelos numéricos para la simula- 


ción de yacimientos. 


El desarrollo de computadores digitales, trajo como consecuen 


cia la posibilidad de resolución de problemas cada vez más complejos, 
especialmente en la solución de ecuaciones diferenciales, imposibles de 
resolver analíticamente (ecuaciones no-lineales). En la industria pe- 
trolera se habla comunmente do Modelos Matemáticos o simplemente Mode- 


los Numéricos de Yacimientos??, 


En geneval,los modelos numéricos de yacimientos consisten еп 
la solución numérica de las ecuaciones fundamentales de flujo ® para 
cada una de las fases considerados en el sistema y considerando iayec- 
ción, producción, intrusión de agua, etc. además de la información 
misma sobre propiedades de la formación y fluidos del yacimiento. De 
acuerdo a las ecnaciones fundamentales de flujo usadas, se conocen va- 
rios tipos de modelos: (a) de una fase y una, dos o tres direccionas, 
(b) dos fasos y una, dos o tres direcciones у (c) tros fasos y una, dos 


o tres direcciones. 


Las ecuacioens fundamentales de flujo se escriben en forma le 


2 


diferencias finitas, tanto en espacio (direcciones) como en tiempo. Se 
coloca una malla de diferencias sobre el campo de petróleo о yacimien- 
to como ilustra la Fig.62- А cada punto de la malla (nodo) se le asig 
na un valor de permeabilidad, espesor, porosidad, contenido de fluidos 


(saturaciones iniciales), elevación, y presión. A los fluídos (petró- 
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leo, gas y agud, en caso de un modelo trifásico) se le asignan valo- 
res de factores volumétricos, solubilidad y densidad, También debe co 
посерѕе la compresibilidad de la formación. Las fuerzas eatre roca y 
fluídos se conocen de las siguientes funciones, dependientes de la sa 
turación: permeabilidad relativa para cada fase, presión capilar entre 
petróleo y agua, presión capilar entre agua y petróleo. Otra inferma- 
ción necesaria proviene de los pozos e incluye: intervalo productor ra 
ta de producción de petróleo, gas y agua como función de tiempo y pre 


sión observada como función de tiempo. 


En un proceso de simulación se asignan ratas de producción рог 


pozo para cada una de las fases. Luego se considera una diferencia Fi 


Limite del Campo 


Fíg.62 - Malla superpuesta en un campo de petróleo. 
о pozos, 
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nita de tiempo y se calculan nuevas presiones y ratas de producción. 
Estas se calculan de las saturaciones existentes en cada punto (nodo) 
de la malla. El proceso se repite para determinado. número de interva 
los de tiempo y en esta forma, se calcula tanto la historia de produc 


ción como de presión por pozo. 


Una de las aplicaciones más importantos de modelos numéricos 
de yacimientos es el estudio de la historia y predicción del comporta 
miento de un yacimiento. En cl estudio de la historia del yacimiento, 
se trata de reproducir la misma a través del modelo para lo cual se 
calculan presiones y ratas de producción como función de tiempo y se 


comparan con los valores observados. Si no coinciden dentro de una a 


proximación aceptable pucde ser nocosario modificar «algunos datos has 
ta que la historia observada y calculada coinciden con la exactitud de 


seada. 


Para la predicción, se asignan ratas de producción, а 105 po- 
zos о al campo en total lo mismo que ratas de inyección on puntos de- 
terminados. Соп este esquema de producción e inyccción se estudia el 
comportamiento del ссоро empleando la información necesaria usada ра- 
ra el análisis de la historia. En esta forma, pueden evaluarse va - 
rios esquemas (programas) de producción e inyección. Introduciendo 


consideraciones económicas, puede wnleccionarse la forma óptima de ex 


plotación (producción) de un campo o yacimiento. 


INTERFERENCIA БИТВЕ С 


MPOS LOCALTYADOS ПМ UN MIS'O АСОТРЕКО. 


А veces se encuentran zonas donde varios campos de petróleo 


B 


estan localizados dentro de un mismo acuífero, es decir, comparten o 


estan afectados por un acuifero común. Tales campos se encuentran 
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en comunicación hidrodinámica; por consiguiente, la producción de flui 


do de uno de estos campos, trae como consecuencia una declinación de 


А 


presión que se transmite a otro и otros campos а través del acuífero, 
ocurriendo un proceso de interferencias de presión entre campos. Este 
comportamiento, puede conducir a una reducción en presión en puntos lo 


calizados a varios kilómetros de un campo en producción. 


“Uno de los ejemplos más típicos de interferencia de campos en 
un acuífero común зе refiere a la formación Woodbine en el oriente del 
Estado de Texas, E.U91, En esta formación se encuentran mas de doce 
campos importantes de petróleo con empuje hidráulico activo, compartien 
do el mismo acuífero. Entre estos campos se encuentra el famoso campo 


Kast Texas. 


La producción de varios de estos campos (principalmente cl cam 


Toxas, descubierto en 1930) trajo como consecuencia que la pre- 


sión inicial del campo Hawkins (ver Fig.61), descubierto en 1940, fuera 
280 1рса por debajo de la presión inicial del аси своб. ¿sto indica 
que el efecto de interferencia en cuanto a presión se refiere, puedo se 
y debe tenerse en cuenta en el estudio del comportamiento de campos si- 
tuados en un acuífevo común. Otros ejemplos tipicos de campos que com- 


parten un acuífero comúa pueden сисопегарзе en la literatura5?,58,59. 


Uno de los mátodos más comunes para evaluar problemas de inte: 


26, 


ferencia en un асі Сего común, es 91 descrito por Mortada te mt, 
do tiene las siguientes suposiciones: (1) cada campo dentro del acuífe: 
го común, tiene forma circular y la caída de presión determinada en e. 
contacto agua-petróleo de cada campo representa la caída de presión рг 


modia en el campo, (2) acuífero infinito,(3) la compresibilidad efectiv 


(eutep?) es constante e iudependiente de presión, (^) formacióu uniform 
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es decic, valores constantes de $, k, h y también constantes las propia 


dades de los fluídos. 


Considérese un acuífero donde existen tros campos de petróleo 
А, В, С, tal como lo indica la Fig.93. Se conoce la distancia entre 
les campos, AB, АС, ВС y los radios de los mismos, Ria, Rip,Ric- De- 
bido a la producción de cada uno de estos campos, ocurre interferencia 
entre sí a través del acuífero. а producción del campo Á trae como 
consecuencia una caída (disminución) de presión cn los campos B y [un 
Similarmente la producción del campo B lleva consigo una pequeña cai- 
da de presión en el mismo campo B y una pequeña caída de presión en 


los campos A y C. Igual cosa sucede con el campo C. 


ACUÍFERO 


Campo 


c 
Rie 


Pig. 63 - Esquema de un acuífero común que contiene tres 
campos de petróleo A, D y C. 
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Sea, (рї-р)дд, la caida total de presión (a partic de la con- 
dición inicial) en el Campo A debido a su propia producción.  Igualmen 
te, (pi-p)gp Y (pi-p)gc, las caídas totales de presión en los campos B 
y С debido a sus propias producciones. Además, (pi-p)a4g es la caída to 
tal (a partir de la condición inicial) de presión вп el Campo A debido 
a la producción del Campo B. Igualmente, (рї-р)дс es la caída total de 
presión en el campo Á debído a la producción del campo C y así sucesi- 
vamente en caso de que existieran más campos. (pi-p)ya4 Y (рі-р)вс 59 
rán las caídas totales de presión en el campo B debido a la producción 
de los campos A y C y finalmente, (рі-р)сд Y (pi-p)gg son las caídas 
totales de presión en el campo € debido a la producción de los campos 
A y B. De lo anterior puede escribirse la caída total de pvesión en 


el campo A, definido por (pi-p),, 
(рі-р)д = (рі-р)дд + (рі-р) адв + (vi-pdac (206) 


Similarmeuto, las caídas totales de presión en los campos 


вус, 
(pi-p)p = (рї-р?вв + (pi-p)na t (рі-р) зс (207) 
(рі-р)с = (Pi-PJec + (рі-р) сд + (рі-р)св (208) 
Los términos que llevan el subíndice mozclado (AB, AC, BA, 


etc) se denominan "cambios (caídas) de presión debido a intecfsren- 
cia" para distinguirlos de los términos con subindice de la misma 
letra (AA, BB, CC) que representan el "cambio (caída) de presión de 


bido a la misma producción del campo". 


Para el cálculo del cambio de presión debido а la misma pro 


ducción, (pi-p)aa, (pi-p)pgp, cte., se hace uso do la solución de van 
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Everdingen y Hurst para sistemas infinitos discutida anteviormente, 
Ec.(37), cuando ryzi. 105 cambios de presión debido a interferencia, 
(рі-р) дв» (рі-р) рд» etc. se obtienen a partir de la solución de Mor- 
tada26 para sistemas infinitos, también discutida anteriormente, Figs. 
20, 21, 22, 23 y Tabla 111, para el valor apropiado de гр. Este valor 
se obtiene dividiendo el valor de la distancia entre los campos consi- 
derados por cl valor del radio del campo que causa la interferencias 
Con refevencia a la Fig.61, por ejemplo, si se está estudiando el efec 


to de la producción del Campo B sobre el Campo А, para calcular (pi-P) 8 
, 


se usará рух ÁD/R;g; en cambio, para evaluar el efecto de là producción 
D iB 
del Campo А sobre el Campo B, (pi-p)sj» se usará гр = AD/Rja, У asi su- 
Pi-P/BA р ЗА 


cesivamente. 


Recudedese que las soluciones mencionadas en el párrafo ante- 
rior se refieren al caso de rata terminal constante. En este сазо debe 
usarse la rata de intrusión y aplicar el teorema se superposición para 
considerar variación en la rata de intrusión, Ec.(84), 

141,2 q 


n 
Ар = (pi-p) = —— - Y 
kh j= 


, (ој - ера) ap Со Е; 1p] (эз) 
donde to y eg боп iguales a cero. Para el caso de calcular los valo - 
ros de (рї-р)дд» (рї-р) вр» etc., las presiones «dimensionales Ap(ty) 

se obtienen de la Tabla II, Pig.19 o Есз.(81) y (42). Para el caso de 
calcular los valores (pi-p)an, (pi-p)ga, ete., las presiones adimensio 
nales se obtienen de las Figs. 20, 21, 22, 23, Tabla III o Ecs.(u9) y 


(50), para el valor correspondiente de rp. 


Debe tenerse en cuenta, que para yacimientos con empuje hi- 


draúlico, la rata de intrusión de agua más la rata de expansión de hi 
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drocarburos es igual a la rata de producción de petróleo. Еп general 


puede escribirsc9?, 


du ан dN ан 
р р 
EP, a e а во —— + (R-Rs)Bg — + By ——— (209) 
dt ах dt dt 
ах ам ан 
р 
e = By = (R-R53)Dg р Ву (230) 
dt dt dt 
e = ао [Bt + (R-Rgi)Bg)] + Bu Qu (211) 


En caso de que el yacimiento se encuentre por encima del punto de bur- 
bujeo y si no existe producción de agua, Bt-Bo, y la ecuación anterior 


puede escvibirse, 
e = doBo (212) 


En base a јаз ecuaciones anteriores puede convertirse fácil- 


mente чо a e o viceversa. 


Las Ecs.(208),(207) y (208) permiten calcular las caídas tota 
les de presión en los campos A, B, C, cuando se conocen las ratas de 
producción y por consiguiente intrusión. Muchos otros problemas pueden 
resolverse, de acuerdo a los datos disponibles. Por ejemplo: se ticnen 
dos campos en un mismo acuífeeo y se conoce la historia de producción 
de A y la historia de presión de B. Para hacer la predicción en este 


caso, deben hacerse los siguientes pasos: 


1)- Con los datos de historia de producción de A se determina 


(pi-p)aag У (bi-p)naj- 
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2) - Dıl dato (р1-РЭв4 У el valor (pi-p)paj 3® calcula. 
(рї-рЭвв» 


(pi-p)pnj = (ri-p)ej 7 (pi-p)paj 


3) - Utilizando la Ec.(95), se calcula la intrusión a сада 


tiempo (período) para el Campo B, 


21 
маву 7 1,119 Феб Rig 14 Ap; o[cep-t5-400) (95) 
j= 
pi valor de Ар = Gi-punjti 7 (рі-р)ввј. Соп estos valores se calcu 


la la historia de intevsión del Campo В, nj 


њу - Su datormina, (pi-P)aBj» aplicando là Ее. (81) y con los 


valores de Cp; detarminados en el Paso 3) ara rpaAB/RiB. 
Bj y P D i 


m 
1 


Zinalmonte, 


(pi-pdaj = (рі-р) АА; + (рі-р) АВ (213) 


Robinsont, presenta un programa de computadores para corre- 
gir la ecuación de balance de materiales, escrita en la forma de de la 
Ec.(135), debido a la interferencia de 01108 campos localizados en un 
also acuífero. En caso de interferencia, a la Ec.(136) debe agregár 
«sole otro térmivo que tome en cuenta los efectos de interferencia. 
El trabajo de Robinson presenta un ejemplo, donde considera tres cam- 


pos en un misno aculfero. 


11,- 


12.- 


13.- 
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APENDICE В 


Determinación del Cambio de Presión Adimensional, Ap(tp), en Termi- 


и 


nos de tp, para el Caso 3 de Sistemas Radiales, usando la Función 


128, 


F(tpa) de Matthews, Brons y Hazebroek, МВ! 


La ecuación para calcular la presión a un tiempo t, de un 
pozo situado еп un sistema radial limitado y que fluye a una rata 


constante puede expresarse por?8, 


2 
pnt) = pi = хн -Eil- PUE ) + Yita? (1B) 
umfkh tkt 
donde tp, está definido por, 
kt kt 
"^ auri ^ uen x 


У (Ерл) es solo función de tp, para valores relativamente grandes de 
тер» tal que 1/12, Zo, y representa, al multiplicarse por див/чикћ 
la caida extra de presión comparada con la solución para yacimientos 


infinitos, debido al límite del yacimiento. (Recuérdase que en el cas 


de intrusión de agua a un yacimiento, г, = Ri). 
МВН?8 desarrollaron una relación para calcular Y(tp4) da- 
da por, 
ЕК(крд) = Зебра ~ ҮС) (3B) 


La función P(tpa¿) se representa grálicamente en la Fig.1B en térmi- 


nos de рд, 


Introduciendo la Ec.(3B) en la Ec.(1B), 
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{ 


с 0i3u3J9j3Ej ¿0/69/10 обима 
po cdwei о 2 


I) y par ар 9012373 “Bi б 


to 


= © 


VG, IVNOISNAWICY CONSUL 


INAJ 


NOI 


ја 


} 


qua 


дес? 
22) + ~ F(tpa) (4B) 
DA DA 
њткћ tkt 


~ ——— | -Eil- 


p(ry,t) = p 


H- 


van Everdingen y Hurst? definenAp(ty), Ec.(35) del texto, 


2ffkh 


Ap(tp) = g - Pot] (58) 


quB 


recmplazando la Ec.(*B) en la Ec.(5B), 


1 gucr? 
Apltp) = ——-Ei(- -5-) + atga = Elpa) (8B) 
2 цк 
уа que, 
kt 
tp = —— (тв) 
paer 
luego, 
tp 
tpa = -- 2 (ов) 
Цан? 
y la Ec.(6B) puede escribirse, 
1 1 
Ap(ty) = — | Г (- ———) + ТЕЛА - FUtpa) (98) 
2 "to 


La ecuación anterior se voa en el cálculo do Арс), para diferente 


valores de rep» cn función de tp» 


51 tpa 20,1, puede escribirse??, 


СЕ А) = 10(0ү5рд) (108) 
donde 9. 
уе 
CA mo 5 (113) 


PU 
donde, d = constante de Euler . Por tanto,la Ec.(9B) pue 


de escribirse, 


1 1 
ЗЕ ===) + "үд = 1n(31,6tpa) (128) 
Ht р 


Apítp) 


Esta expresión es la Ec.(61) del texto. 
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APENDICE C 


1 - Cálculo del Integral f(p,t) de la Ecuación de Schilthuis, 


Ec.(101). El integral de esta ecuación se define por, 


t 
f(p,t) = (pi-p)dt (1c) 

o 
La Fig.55 presenta la historia de presión de un yacimiento 
(contacto agua-petróleo). De la historia se conocen los valores de 
presión pi, Ру» Од» P3 <+.. Pn-2» Pn-1» Pn A los tiempos 0, tq, 12, 
ta .... 1-2, fn-1, tys El integral f(p,t) representa el área por 


uncima de la curva, es decir, el área comprendida entre la curva con 


tinua y la linea punteada horizontal supevior. 


а E | 
El método del cálculo del integral f(p,t) puede hacerse де 

i ета : ү os 
varias maneras de acuerdo a la definición del diferencial de presion 


Ap. Si el Ap se define, 


де, $ ТЕ (20) 
donde 4 = 1, 2, dy v. N ‚ luego el área por encima де la curva 


uede calcula or la sumatoria de una serie de áreas de trapecios 
р , 


excosto para ol primer At = (t,-0) que será el área de un triángulo. 


Si además se define, +. 


Ati = (t3-t3-1) (30) 
donde j = 1, 2, 3, +... n (donde to = 0), el área encima de la curva 


Ar 
f(p,t) = CAL) Ata + (— 


2 


Ура МАР 


2 


+ (- 


f(p,t) = 


) Ata 


Apt др, 


2 


1 n-i 
ES Ata + Ap 
2 јЕ1 


Apot Ap 
A АЕ 


2 


Reemplazando las Ес5.(20) у (30) en la Ec.(4C), 


1 


1' At; + Аура 


(nc) 


f(p,t) = [eere + (ру-р1+рї-Р2)(%2-%1) + (pi-p2*P4-P3) 


(03-12) To. 


donde Ро = Pie 


decir, los intervalos de tiempo son iguales, las Ecs.(5U), (7с) 


+ Cpi-Pn-itfPi7Pp)(tn-tn-1) 


(2pi-bj70j4-1) (t3-t3-1) 


(8C), pueden escribirse, 


С(р,%) 


2 í 


n 


> (рү"Рү.-1051-5ү-19 


81 


351 


s 1 өч 
| gc "з 4 руку.) 
15 


| (2n-1)pi ~ 


(60) 


(76) 


(8с) 


Si el At; = (tj-tjg.1) es constante e igual а At, es 


y 


(106. 


(tic 


(120 
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Si en el análisis de la historia del yacimiento о en la 
predicción es necesario utilizar la Ec.(119), el diferencial de pre- 
sión usado en esta ecuación se define en la forma indicada en la Pig. 


55, Ec.(122) del texto, 


1 
Ару = —(Ф;-2-Р) (130) 
2 


donde j = 1, 2, 3, .... n y pag = P4 = Pj. Si además, el Ar; se de 
fine por la Ec.(3C), en este caso el área por encima de la сийуа es~- 
tará dada por la suma de las áreas dc los diferentes rectángulos que 


aparecen on la Pig.55, 


f(p,t) = Ор, At4 + CAp1+Ap2) Ata + ( Ар1+ДАр2+ Ара) Аа + 0... 


+ (Apat Ард? Араз + o Ари) Ар (140) 


n n 
#(р,ї) => ( 2 Pj) АЕ; (150) 
j=1 Isi y 


Ruomplazando el An; de la Ec.(13C), se obtienen las Ecs.(7C) y (86) 


deducidas anteriormente. 


Si los int:rvalos de tiempo son iguales, АЕ: es constante e 


igual а At, la Ec.(15C) puede escribirse, 


f(p,t) = (160) 


и 


~ ~ 
At 2 Of) Др; (176) 
151 


Reemplazando los valores de Ар de la Ec.(130), se obtiunen las Fes. 
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(8C) a (12C) deducidas anteriormente. 


Finalmente, si el diferencial de presión se define por, 


donde j = 1, 2, 3, .... n. Para j = 1, Po = Ра. Como antes 61 va 
lov de Ах, se define por la Ec.(3C) y el área por encima de la cur- 
va estará dado por la suma de las áreas de rectángulos y triángulos 


en la siguiente forma, 


Арі AP2 A 
f(p,t) = ---г- At, + CApy* шалны At, + САру+ Apart — ta 
2 
Ap 
но... + (Apat Apot Араг +... o (190) 
1 
хишиг (200) 


Reemplazaado el valor de Ар де la Ес. (180), se obtienen las Гез. 


(те) y (86). 


Si los intevvalos de tiempo son iguales, £M es const 


rual а At, la ecuación anterior puede чзлрірірзе, 


" л Aem 
f(p,t) = === Ару +2 ри (210) 


Finalmente reemplazando Ар; por la Ее.(180), se obtienen Jas Eco. 


(100) a (120) 


2- Cálculo dol Integral f(p,t) de la Есллс 6) de Pur 


dificada, E (109). El integral de esta ecuación se define por, 


t 
p 
аё ғ A d E (220) 


log(at) я log(at) 


f(p,t) = 


Como antes, se conoce là historia de presión, o sea una se 
rie de presiones ру, Ру» Pg '''' Pn а los tiempos t4, t9, 53 s». ta 
respectivamente. Para obtener gráficamente el valor del integral de 
ја Ec.(22C), basta calcular los valores de Ap; de acuerdo a la Ec. 
(20), dividir cada Др; por el loglaty) correspondiente y graficar 
los valores de Apz/loglaty¿) como función de t. Zl área bajo la 
curva resultante, será ol calor del integral. De nuovo, el tiempo 
estará dividido en intervalos de acuerdo а la Ec.(3C). El Дева ba- 
jo la curva será la suma de las áreas de una serie de tveapecios, ex- 


cepto en el primer intervalo que será el área de un triángulo, 


1 A 1 Ap АР 1 
f(p,t) = -— MES деу bmi is oce) Ata + — 
2 log(at4) 2 loglat1) log(at2) 2 
AP ДР: 1 Ара - 
( 2+ A At ен S TEE 
log(at5) log(at3) 2 log(atgy-1) 
APn 
в —_ ) Ata (230) 
loglatp) 
n-1 
1 АР, Орд 
= — | س‎ — Atp + e ——( ДЕ; + АЕ; +1) (240) 
2 log(atp) e loglatj) 
1=1 


Reemplazando los valores de Ap y At de las Есв.(20) y (30), 


n-1 
рі-р —  Pi-Pj 
f(p,t) = X 2» ЛЭХ а ES A о) 025€ 
и сг )t1 7j-1 
2 loglatp) log(atj) p 


ja 
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Si los intervalos de tiempo son iguales, es decir, At; = 


constante = At, las Ecs.(24C) y (250) pueden escribirse, 


nzi 
f(p,t) = ¿A Хра +2 > „Арі (260) 
2 loglatp) Sar Іов (аё) 
п-1 


и 
Ї 
! 
| 
i 
| 
+ 
~ 


At Pi-Pn pi-pj 
хэс - —— (270) 
2 loglaty) loglaty) 


¿oria f(p,t) de la Ecnación de van ~ 
LM GI ы ыг ГЕЙ ы ДИА. a па. ЫА ду 


Everdingen Hurst, Ec.(119). La sumatoria de esta ecuación 
Ever Y 2 PLE 


definida por, 


APj У КИТА (280) 


donde tq es el tienpo a la cual se doses calcular Ja sumatovia £(p,0. 


El valor de Ар; se obtiene en la forma indicada en la Fig.55 у se¬ 


gûn se discutió en el texto, Ec.(122). 


Doterminación del Petróleo o Gas Originalmente en el Yaci- 


miento, М ó Gy la Constante de Intrusión de Agua, C, por Mínimos Cva 


drados?9, 


Las Ecs.(136) y (144) del texto pueden os 

Y =CX+ A (1D) 

donde А, constante; para yacimientos de petróleo es Ч y para yaci- 
mientos de gas es 6. Los valores де C y ^ (constantes) pueden ge- 
terminarse por mínimos cuadrados, siempre y cuando se tenga una se- 
rie de valores de Y y X. Cuando la ecuación за escribe en Ja forma 


anterior so habla de LINEA DE REGRESION DE Y SOBRE X. 


Considerando la ecuación znteior para yacimientos de ре- 
tróleoo, А=Н, escribiéndola para заа serie de valores de Y y X, y Su- 
nando las ecuaciones resultantes puede escribirse la denesinada pri- 


mera seuación normal, 


n 215. 
Ух = ср Хү + он (20) 
1 1 
Si cada una de las ecuaciones anteriores se multiplica por Х, y se 


suman las ecuaciones resultantes, 50 obtiene la segunda ecuación nor 


mal. 
e 2 n 
SX E Хр 82 Xi (3D) 
Y 1 


despejando C y N de estas dos ecuaciones, 


n. لگ‎ 
a2 XQ Yi - (QN 10 
1 1 1 
с = - = 
Exon 
n Жу = х 
1 Los 
л. 2 n n, n 
" СЭН. Yi) (223068 Үү 
n n 
мр ка -($ хр" 
1 1 


o también si el valor de С se reemplaza en la Ec.(2D), 
1,2 a 
N= — (> - can 
n 1 


La Ec.(1D) puede arreglarse en forma diferente. 


dieado por C y transponiendo términos, 


3: 1 
X EN 8 
c c 
Si se define, 
1 
б® = — 
с 
1 
N! = - —— N = C'N 
e 
la Ec.(7D) puede escribirse, 


X = C'Y + N” 


Cuando la ecuación de la linea recta за escribe 


(14D) 


(5D) 


(6D) 


Divi- 


(8D) 


(9D) 


(10D) 


en esta 


forma, se habla de LINDA DE REGRESION DE X SOBRE Y. Los valores de 


C' y N' son constantes, 


En forma similar al caso anterior pueden escribirse las e 


cuaciones normales, 


(11D) 
n цээр n 
Уха тр = oh Y] +N Y (12D) 
1 1 1 
de donde puede obtenerse C' y N' 
p. n. 
Xi Ys = ( х)о Y 
1 1 
C! = — ЈЕ. SE (13D) 
n 2 n 
i Yi С Y)? 
1 1 
eL qn cm. qu. 
ФРУ к) - ОИ 54 
E i 1 1 
y? = == AAA (14D) 
п о n 
n2 Yi - (Ty? 
1 1 
o también se el valor de 0! se reemplaza еп la Ec.(11D), 
n n 
"Ww. _ Ын 
N's -ži бар (15D) 


Una vez obtenidos C' у N' , puede determinarse Су Ма travós de las 


Ecs.(8D) y (9D). 


Cuando se trata de yacimientos de gas, se cumplen exacta - 
mente las mismas ecuaciones, solo que en lugar de caleular Ну N' el 
valor correspondiente será б y G'. Además debe tenerse cuidado al 
calcular Y y X, usar los definiciones correspondientes de acuerdo al 


tipo de yacimiento y tipo de intrusión. 


YACIMIENTOS CON EMPUJE HIDFAULICO 


RAMIRO PEREZ PALACIO 


ERRORES 
Página Línea Dice Debe Decir 
1 5 concecnencia consecuencia 
2 24 variado variada 
3 17: son resultado es resultado 
5 14 lates tales 
5 24-25 energis energía 
16 23 desplacen desplacen 
17 t coverdenadas coorisnadas 
17 13 pororso poroso 
17 22 (13/1312 = L/T) (3/7/17? = 1/1) 
22 Ec.(11) 2р. 1 р LC UM 
2 р др © ср 24 
25 13 pude puede 
ЭР Эр 
32 12 (——) = 0 ——), = 0 
эх ( эх e 
eo ° 
33 Ec.(28) bi 5 
n-1 n=1 
35 2 imitario unitario 3 
37 15. . presentandas presentadas 
48 1 aporzimarse aproximarse 
50 22 luego Chatas?, luego por Chatas?, 
57 11 Ayúndiece Apéndice 
58 Ec.(58): El numerador dentro del signo integral debe ser: 


(1 -се-ч?®р) (л, (u)Yelurp) - үз (а) баер) Ја 


' 

58 última página siguiente. página 65. 

65 El valor de л que entra en las definiciones da - 
ierf(x), erfc(x) y erf(x) debe sor: фт 


65 Los límites de integración de la definición de 
23 о 
erfc(x) debe ser: | y no | 
И х о 
77 18 con las raices son las raices 


97 7 tp л 5,5(0, 5)? cp > 0,5(гср3 


143 
145 
145 


15 


147 
148 


. 143 


157 
158 
159 
161 
162 
164 
165 
166 


167 
167 


169 
169 


Ec.(118): 


o 


19 


20 


Ec.(130): 


Ес.(160): 


13 
13 
última 


uw ore 

Dice Debe Decir 
dividada dividida 
histroia historia 
obtenrse obtenerse 


Debe ser: Ме = NBoi ~ (N-Np)Bo + Мр - Wi 


eucaciones ecuaciones 
eliminana eliminará 
concer conocer 
calcularse hacerse 
conocen conoce 
series serie 


días, meses 
ya que Cg se 


días meses 
ya que se Cg se 


La última sumatoria a la derecha del nume 


n 
тайор debe ser: Е Cgj tj log(tj) у по 
n 1 

E Сеј log (tj) 

1 
represetna representa 
disctutido discutido 
conocen conoce 


La parte derecha de la ecuación para Ард 


1 
debe ser: -- (pj - рз) y no EL Bg =p 
2 2 


convertir 
conocimiento 


converitr 
conociemiento 


Debe ser: а = (Cup - суу? 


calculado cálculo 
obtenrse obtenerse 
razoanes razones 
infoermación información 
fraficados graficados 
que para que para cual 
inconaitivo insensitivo 
Mabwen McEwen 


Debe ser: X = GpB, + Wp ~ Wi 


g 
sencible sensible 
calor valor 
:imicuto yacimientos 


Línea 


Ec.(167): 


~ 
o o £m 


> 


юк 
ona 


хл 


~ 
~ 


Q(tg) = Atp)” 


cosntante 
cuanta 
yacimiento 
taels 
dividdiéndolos 


ño (ро - Pg) 


facores 

en se desee 
expredada 
productivifad, j 
en 

produtividad 
prevenier 
hidrüalico 
facotr 

a la que 
conocerce 

Un aumento 
mecionado 
Fácilemente 
correspidiente 
aumneta 
constante, Figs 


Debe Decir 


21433) з (Ае р)" 


constante 
cuenta 
yacimientos 
tales 
dividiéndolos 


Ар = - (po - Pg) 


factores 

en que se desee 
expresada 
productividad, 4 
el 
productividad 
prevenir 
hidráulico 
factor 

à la cual 
conocerse 

Una disminución 
mencionado 
facilmente 
correspondiente 
гмпопфђа 
constante. 


Debe ser: 523 = Z(JAt) - 2001-1345) 


repetin 
correción 
a la que 
En 


ador 
mircotaradios 
elúctircas 
cocveinte 
ecuacioens 
896; 

anreior 


repiten 
corrección 
en la que 
El 
eléctricos 
condensador 
microtaradios 
eléctricas 
corriente 
ecuaciones 
896, 1963; 
anterior 


